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NOMENCLATURE 


B  streamsheet  thickness  in  normal 
direction 

E  total  energy 

-* 

F  force  vector 

Ff  distributed  friction  force  vector 
H  total  enthalpy 

I  rothalpy 

J  Jacobian 

M  Mach  number 

Q  heat  sources 

Qg  additional  heat  source 

S  stream  surface 

T  temperature 

V  average  absolute  velocity  vector 

W  average  relative  velocity  vector 

a  speed  of  sound 

b  streamsheet  thickness  in 

3 

£  -direction 
b  channel  width 

d  blade  thickness 

ei  unit  base  vector,  i  =  1,2,3 
fe  external  force 

fg  body  force 

h  static  enthalpy 

hi  metric  coefficients,  i  =  1,2,3 
k  thermal  conductivity 
k  averaged  kinectic  energy 
m  meridional  distance 
it  normal  to  streamsur f ace 
p  pressure 


p*  rotary  stagnation  pressure 
r  radius 

r  gas  constant 

s  entropy 

t  time 

u  x-velocity  vector 

v  y-velocity  vector 

w  relative  velocity  vector 

(J  relative  flow  angle 

8*  blade  angle 

Y  specific  heat  ratio 

Y  station  angle 

£ 1  blade  lean  angle 

contravar iant  base  vector 
covariant  base  vector 
C  vorticity  vector 

0  angular  coordinate 

£ a  curvilinear  coordinate,  u  *  1,2,3 
P  density 

o  meridional  flow  angle 
o  stress  tensor 

t  shear  stress  tensor 

=S 

T  secondary  stress  tensor 
^  stream  function 

il  volume 

oj  angular  velocity 

to  pressure  loss  coefficient 

V  streamwise  gradient  operator 

V  gradient  operator 
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Superscripts 


averaged  quantity  A 
density  averaged  quantity  A 
A!  fluctuating  quantity  with  respect  to  A 
A"  fluctuating  quantity  with  respect  to  A 

Tq  relation  stagnation  quantity 

Subscripts 

r  cylindrical  coordinate  components 

0  cylindrical  coordinate  components 

z  cylindrical  coordinate  components 

m  stagnation  quantities 
p  pressure  side 

s  suction  side 
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INTRODUCTION 


The  computation  of  turbo.nachinery  flows,  particularly  in  the  multistage 
v  :nf igur  <tion ,  is  based  on  a  complex  set  of  approximations  and 
interaction  with  empirical  data.  It  combines  an  approximate  space 
description  with  a  representat ion  of  the  shear  stresses,  which  reduces 
to  their  effect  on  the  stagnation  pressure  losses,  or  entropy 
production,  provided  by  empirical  information. 

In  this  report  a  rigorous  derivation  of  these  components  is 
pre3ent€.'d . 

Fart  1  contains  a  systematic  derivation  of  the  conservation  laws  for  the 
two  quasi  three  dimensional  space  approximations,  namely  the 
st reamsurf ace  approach,  initially  introduced  by  C.H.Wu,  and  the 
a veraging  procedure . 

Part  2  describes  the  formulation  of  the  Through-flow  models  and  their 
coupling  with  the  blade  to  blade  flows,  including  a  discussion  of  their 
numerical  discretization.  In  particular,  the  derivation  o.f  the 
streamline  curvature  method  is  presented  in  details  and  its 
implementation  can  be  compared  with  the  streamf unction  approach.  In 
both  cases  the  limitations  of  the  methods,  in  presence  of  supersonic 
flows,  are  discussed  and  their  origin  are  pointed  out. 

This  work  has  been  performed  while  the  author  was  holding  the  NAVAIR 
Research  Chair  at  the  Naval  Postgraduate  School.  The  stimulating 
atmosphere,  and  discussions  with  Professor  R.  Shreeve  have  greatly 
contributed  to  the  realisation  of  this  work. 
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PART  1  :  THE  SPATIAL  LEVEL  OF  APPROXIMATION 


The  most  general  representation  of  a  fluid  flow  is  obviously  the  fully 
three-dimensional  one,  since  it  corresponds  to  the  dimensionality  of  the 
real  world.  But  needless  to  say,  this  representation  puts  in  general 
very  strong  requirements  on  computer  storage  and  computational  times 
and,  although  three-dimensional  computations  are  to  be  applied  in 
increasing  measure  in  the  coming  years,  many  physical  situations  can 
still  be  described  with  reasonable  accuracy  by  an  approach  with  a 
reduced  number  of  space  variables. 

Many  configurations  are  indeed  developed  which,  although  not  strictly 
two-dimensional  in  nature,  can  be  considered  as  not  "strongly" 
three-dimensional  in  the  sense  that  one  can  expect  that  the  variation, 
or  the  magnitude,  of  one  velocity  component  is  less  important  than  the 
variations  of  the  other  two.  A  typical  example  of  such  a  situation  is 
the  wing  of  a  modern  airplane  with  varying  chord  length  and  spanwise 
varying  profile  shape.  Within  a  certain  approximation  which  excludes  in 
particular  the  wing-body  junction  and  the  tip  regions,  the  flow  along 
the  wing  could  be  approached  two-dimensionally .  Alternatively,  a 
cylindrical  airfoil  under  a  sweep  angle  with  respect  to  the  upstream 
flow  will  present  limited  three-dimensional  effects  for  relatively  small 
values  of  the  sweep  angle. 

Axisymmetric  configurations  are  similarly  suitable  for  descriptions  with 
a  reduced  number  of  space  variables. 

It  is  clear  that,  neither  for  the  spatial  level  of  approximation  nor  for 
the  other  levels,  is  it  possible  to  put  forward  general  conditions  of 
validity  for  a  given  set  of  assumptions.  Depending  on  the  initial  flow 
conditions,  the  geometrical  configurations  and  the  required  level  of 
accuracy  for  any  system  of  approximations,  one  will  have  to  check  the 
calculated  flow  behaviour  against  either  experimental  data  or  results 
from  computations  at  a  higher  level  of  approximation. 

Hence,  a  fully  three-dimensional  computation  will  allow  to  check  by 
comparison,  the  limits  of  validity  Of  a  two-dimensional  approach,  in  the 
same  way  as  experimental  data  or  a  reliable  viscous  calculation  will 
allow  to  ascertain  the  limits  of  validity  of  an  inviscid  approach. 

For  example,  it  is  well  known  that  for  flow  incidences  along  an  airfoil 
lower  than  certain  values,  the  boundary  layers  will  not  separate  and  an 
inviscid  two-dimensional  approach  might  be  of  sufficient  validity  to 
predict  the  lift  coefficients  with  acceptable  accuracy;  while  for 
higher  incidences  the  occurence  of  large  separated  regions  might  require 
a  three-dimensional  viscous  computation  in  order  to  maintain  the  same 
level  of  accuracy. 

In  the  field  of  internal  flows,  such  as  channels  and  turbomachinery 
passages,  three-dimensional  effects  are  induced  by  the  presence  of 
material  walls  and  the  geometrical  configuration  will  lead  to 
three-dimensional  flow  components,  even  with  a  two-dimensional  inlet 
flow. 
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In  order  to  treat  these  situations  an  intermediate  description  between 
the  fully  three-dimensional  and  the  two-dimensional  can  be  introduced. 
This  intermediate  level  of  approximation  which  can  be  called 
Quasi -Three-Dimensional ,  approximates  the  flow  as  a  succession  of 
interacting  families  of  two-dimensional  flows  along  intersecting 
surfaces.  Appropriate  families  of  surfaces  are  defined,  which  can  be 
treated  as  streamsurfaces ,  along  which  a  two-dimensional  velocity  field 
projection  is  determined.  Obviously  a  fully  three-dimensional 
time-dependent  flow  field  will  require  three  families  of  intersecting 
surfaces  in  order  to  be  determined  completely.  This  corresponds  to 
requiring  three  inter-dependent  scalar  functions  of  three  coordinates  in 
order  to  describe  the  three  component  velocity  field  restricted  by  the 
continuity  equation.  But  only  two  families  of  surfaces  are  required  for 
steady  flows,  see  Section  2-8.  Such  a  representat ion  is  equivalent  to  a 
complete  three-dimensional  description  and  no  spatial  approximations  are 
involved.  The  Quasi-Three-Dimensional  approximation  comes  in  when  the 
number  of  surfaces  is  reduced;  that  is,  when  it  is  considered  that  a 
valid  description  of  the  flow  behaviour  is  obtained  by  neglecting  the 
flow  contributions  along  either  certain  families  of  surfaces  or  along 
certain  members  of  a  given  family  of  surfaces. 

Actually,  one  could  visualize  these  families  of  surfaces  as  generating  a 
system  of  coordinate  surfaces,  such  as  the  *  £,  £2  -  n,  £3  ^  ; 
families  of  figure  1.1.1.  The  projections  of  the  velocity  vectors  in 
these  surfaces  are  computed  on  a  two-dimensional  basis,  and  the  Quasi-3D 
approximation  is  contained  in  the  assumption  that  either  the  projection 
of  the  velocity  field  in  one  set  of  surfaces  is  independent  of  the  third 
coordinate,  with  the  consequence  that  the  flow  in  one  of  the  members  of 
the  family  is  representat ive  of  all  of  them,  (for  instance,  the  flow  in 
a  (£,;)  surface  being  independent  of  n),  or  eventually  that  the  velocity 
field  in  one  family  is  of  negligible  magnitude  with  respect  to  the 
others,  (for  instance  the  velocity  projection  in  the(;,n)  surfaces). 

An  alternative  to  the  stream-  or  pseudo  stream-surface  method  for 
obtaining  a  quasi-three-dimensional  approximation  is  offered  by  the 
Averaging  Method.  This  approach  consists  in  averaging  out  the 
conservation  equations  with  respect  to  a  chosen  coordinate,  for  instance 
the  £3  direction,  obtaining  equations  with  the  remaining  coordinates  as 
independent  variables . 

Two-dimensional  equations  are  obtained  in  this  way,  representat i ve  of 
the  average  flow  with  respect  to  the  £3  coordinate,  but  containing 
geometrical  terms  function  of  the  averaged  space  direction.  Clearly, 
limits  of  integration  have  to  be  defined  and  hence  this  procedure  is 
best  suited  for  internal  flows  with  streamwise  varying  cross-sections. 

This  procedure  can  also  be  extended  to  an  averaging  over  a 
two-dimensional  region  such  as  the  cross-section  of  the  channel,  leading 
to  a  quasi-two  dimensional  flow  description  in  the  direction 
perpendicular  to  this  section  involving  terms  which  are  functions  of  the 
varying  cross-section. 

Although  both  approaches  lead  to  equivalent  representat ions  in  Quasi-3D 
approximations ,  with  differences  only  in  the  way  the  interaction  between 
the  various  families  of  surfaces  are  defined,  the  streamsurface  approach 
can  be  used  to  define  a  description  of  a  fully  three-dimensional  flow 
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whereas  the  averaging  procedure  can  not. 

The  two  models  are  presented  and  discussed  in  the  following  sections. 
It  is  also  clear  from  the  preceding  considerations  that  the  interest  of 
this  quasi  three-dimensional  approach  lies  mainly  in  the  field  of 
bounded  or  internal  flows.  Therefore  we  will  present  this  approach  in 
view  of  its  original  development  in  the  field  oi  turbomachinery  and 
channel  flows. 
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1.1.  THE  QUASI -THREE-DIMENSIONAL  REPRESENTATION  (Q3D)  FOR  INTERNAL  FLOWS  : 
C.H.  WU’S  STREAMSURFACE  METhQD  Wu  (1952),  Wu  (1976) 


The  streamsurface  method  described  in  this  section  has  been  introduced 
by  C.H.  Wu  (1952)  in  order  to  approximate  the  internal  flow  in  a 

turbomachinery  blade  passage. 

We  will  consider  an  arbitrary  channel  bounded  by  an  inlet  and  an  outlet 
surface,  for  instance  a  duct  or  the  passage  between  two  blades  in  a 
turbomachine,  or  the  part  of  a  river  between  two  cross  sections.  A 
typical  configuration  is  represented  schematically  in  figure  1.1.1.  An 
arbitrary  curvilinear  coordinate  system  £l,  C2»  £3  (or  £,  n,  0  is  shown 
with  taken  approximately  in  the  mainstream  direction  and  £2,  f;3 

defining  roughly  transverse  directions.  The  inlet  and  outlet  sections 
may  be  defined  by  constant  surfaces.  A  typical  representat ion  for  a 
turbomachine  passage  is  shown  in  figure  1.1.2. 

In  this  figure  the  direction  of  the  axis  of  rotation  of  the  machine  z 
may  be  identified  with  £ 1 ,  while  the  cylindrical  coordinates  Q  and  r  are 
respectively  to  be  considered  as  £2  and  53  coordinates. 

Representative  streamsurf aces ,  such  as  the  Sx  or  S2  surfaces  in  the 

above  figures  are  obtained  by  following  the  paths  of  the  fluid  particles 
lying  along  selected  lines  in  the  inlet  section.  For  instance,  the 
fluid  particles  along  the  line  AiBx  of  figure  1.1.1  will  generate  a  E, 3 

type  of  surface,  or  S!  surface,  and  the  particles  along  A2B2  will 

generate  a  £2  type  of  surface  (S2  surface) 

More  generally,  the  considered  surface  S  is  described  by  a  time 
dependent  equation 


SU\  C,  t)  =  0 


(1.1.1) 


obtained  by  following  the  paths  of  the  fluid  particles  emanating  from  an 
initial  line.  In  a  general  time  dependent  flow,  S  will  not  be  a  strict 
streamsurface ,  in  the  sense  of  the  velocity  vector  being  normal  to  the 
surface.  This  will  however  be  the  case  when  the  flow  is  steady. 

The  surface  (1.1.1)  can  be  written  explicitly  in  terms  of  the  reference 
coordinate  £3  as  follows, 


c3  =  S3U\  S2,  t) 


(1.1.2) 


where  the  direction  £3  can  be  visualized  as  being  roughly  in  the 
direction  of  the  normal  to  S,  and  where  £l  and  £2  are  coordinates 
defining  the  surfaces.  In  the  following,  we  will  restrict  ourselves  to 
time  independent  surfaces  ;  since  S  is  then  a  streamsurf ace ,  one  has  at 


each  point  of  the  surface  S,  with  n  being  defined  as  normal  to  S 


(1.1.3) 
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for  n  defined  by 


(1.1. H) 


Two  possibilities  are  open  in  order  to  describe  the  family  of  surfaces 
with  respect  to  the  arbitrary  curvilinear  system  of  coordinates,  figure 
1.1.3. 

1)  The  surface  S  does  not  contain  the  and  £2  lines  and  is  situated 
arbitrarily  with  respect  to  the  local  coordinate  system  5°,  (a-1,2,3). 

In  this  case,  the  surface  can  be  represented  by  equation  (1.1.2)  and  the 

-► 

normal  to  the  surface  is  proportional  to  the  vector  n,  defined  by  its 
covariant  components  nQ(-3i^3,  -32£3,  1).  Indeed,  from  equation  (1.1.11) 

with  S  =  C  -  £5U\  52). 

n  =  -  (|^r).e‘  -  £||-).e2  +  l.e*  (1.1.5) 


where  e  ,  e  ,  e  are  the  contravariant  base  vectors  of  the  coordinate 

system.  Note  that  the  components  of  n  as  defined  by  equation  (1.1.5) 
are  not  projected  along  unit  vectors,  since  the  base  vectors  are  not  of 
unit  length.  For  instance,  in  an  orthogonal  curvilinear  system  with 
metric  coefficients  hi,  h2,  h3,  the  length  of  the  contravariant  basis 

vector  e  is  equal  to  1/h  . 

M  a 

“*Ql 

Hence,  one  would  have,  with  e  denoting  the  unit  base  vectors 


+  -1  +1  1  dC  + 2  1 

n  3  hi  W  e  “  h2  dC  e  *  h3  e 


■►u 

=  n  e 
—  a 


(1  . 1  .6) 
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defining  hereby  components  n  as  the  projections  of  the  normal  n  along 

■►ot  a 

the  unit  vectors  e  . 


It  is  to  be  remembered  that  the  covariant  base  vectors  e  are  tangent  to 
ct 

the  £  -coordinate  lines,  while  the  contravar iant  base  vectors  e  are 

normal  to  the  surface  £a*=  constant.  These  two  sets  of  base  vectors  are 
aligned  in  the  same  directions  when  the  coordinate  system  is  orthogonal. 

ii)  The  surface  S  is  defined  by  the  coordinate  lines  V  and  £2 ;  that 
is,  S  is  a  £3  *  constant  surface,  and  the  surface  S  contains  the  lines 
£l  »  constant  and  £2  «  constant.  In  this  case,  it  is  seen  from  equation 

(1.1.5)  that  n«e3  and  if  the  coordinate  system  is  orthogonal,  the  unit 

4  4 

vector  along  the  normal  direction  is  given  by  1  -  h3.n  . 
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1.1.1.  Streamsurface  variations 


With  the  general  representation  of  S  by  equation  (1.1.2),  the 
(covariant)  components  of  the  normal  vector  n,  are  given  by 

||t  =  3,£3  -  "  n,  (1.1.7) 

and 

3  r  3 

g|r  =  a2S3  -  ■  n,  (1.1.8) 


with  n 3=  1  . 

This  allows  us  to  define  the  variations  of  any  scalar  flow  quantity 
along  the  streamsurfaces  S,  as  functions  of  and  £2 ,  by  the  following 
relations,  where'the  overbars  indicate  variations  along  the  surface, 

_L  JL  +  ( 9  (1  •  1.9a) 

at,1  -  at;1  at7 

or  in  condensed  notation 

^ i  =  ^ i C  ) • ^ 3  =  ~~  n j  •  <)  3  (1.1.9b) 

Similarly 

_L  a  wiki) 

as2  -  as2  as2 

or  in  condensed  notation 

a2  »  a2  +  o2s3).a3  =  a2  -  n2  .  a3  (i.i.iob) 

Obviously,  the  variations  a3  with  respect  to  £3  along  the  surface  are 
identically  equal  to  zero,  for  arbitrary  functions  g 


3  (1.1.10a) 

w 


(1.1.11) 


Equations  (1.1.9)  to  (1.1.11)  can  be  grouped  in  the  form 


3aB  =  9a8  - 


3jg 
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a  = 


1  ,2,3 


(1.1.12) 


or  in  vector  notation 


n  . 


9  jg 


(1.1.13) 


for  any  scalar  quantity  g. 

In  the  particular  choice  of  coordinates  where  the  surface  S  is  the  *(5l, 
52)  surface,  that  is  a  surface  53  =  constant,  case  ii  above,  the 

variations  3,  and  32  are  identical  to  the  ordinary  partial  derivatives 

3,  and  32,  since  the  components  nj  and  n2  of  the  normal  vector  are  zero 

according  to  equation  (1.1.5). 


The  gradient  with  overbar  Vg  is  the  projection  of  the  three-dimensional 

operator  Vg  on  the  surface  S.  This  can  be  seen  from  the  definition  of 
the  gradient  of  a  scalar  in  curvilinear  coordinates,  when  the 
coordinates  5*  and  £2  are  in  the  surface  S  ;  case  ii.  With  the  index  Y 
ranging  from  one  to  two,  one  has 


Z  +  0  -  +Y  ^  3  IK 

Vg  =  e  .S^g  =  e  .3yg  +  e  .3sg 


a  -  1 ,2,3 
Y  =  1  ,2 


(1.1.14) 


-¥  -4 

Since  in  this  case  e3  =  n,  this  equation  is  identical  to  equation 
(1.1.13)  with  the  definition 


Vg 


=  e 


9yg 


Y  =  1  ,2 


(1.1.15) 


when  the  f;1  and  52  lines  are  in  the  surface  S. 

Note  that  equation  (1.1.13)  remains  valid  for  any  system  of  coordinates. 

This  approach  can  be  applied  to  obtain  the  streamsurface  variations  for 

-► 

a  vect°r  quantity  F.  from  the  definition  of  the  divergence  of  a  vector 
in  curvilinear  coordinates,  one  has 


-£  -►  -►  a  -►  -+Y  -►  _  -> 

V.F  =  t .3  F  =  e  . 3  F  +  e3.33F 
a  Y 


+y  +  +  + 

t  .3^F  +  n.33F 


(1.1.16) 


If  and  are  orthogonal  coordinates  in  the  streamsurface  S,  the 
first  term  can  be  worked  out  explicitely  as  follows, 


8  - 


(1.1.17) 


-  •+  -fY  ■+ 

V.F  -  e  .3  F 


cY.3  (F°c  ) 
T  a 


; [  3v(FYh  jh2 )  ♦  F33,(h,h2)] 


3v(FYh,h2)  ♦  3j(h,h2) 


h,h2  t 


h,h2 


The  first  term  is  the  two  dimensional  form  of  the  divergence  operator, 
while  the  second  term  represents  a  contribution  to  the  ffstreamsurface,! 

-¥ 

“ 

divergence  operator  V.F,  arising  from  the  deformation  of  the  surface  in 
the  third  direction. 


*►  + 

This  contribution  vanishes  when  the  vector  F  lies  in  the  surface  (F.n-0) 
or  when  the  streamsurface  is  independent  of  the  third  coordinate  £3. 

In  both  these  cases  the  flow  can  be  considered  as  two-dimensional. 

With  the  definition  (1.1.17),  equation  (1.1.13)  can  be  generalized  to 
vector  quantities  as 


-  -►  +  +  +  + 
V.F  =  V.F  -  n. 3 jF 


(1.1.18) 


There  are  various  ways  to  obtain  the  flow  equations  considered  along  the 
streamsurface  S.  One  way  can  be  termed  as  an  "algebraic  method"  and  was 
originally  followed  by  Wu  (1952)  with  cylindrical  coordinates,  and 
further  extended  to  formulations  in  arbitrary,  curvilinear  coordinate 
systems  by  Wu  (1976). 

This  approach  consists  in  introducing  directly  the  streamsurface 

derivatives  3i  and  32  in  all  the  equations  of  motion  written  in  the  E? 

coordinate  system  in  place  of  the  ordinary  derivatives  3i  and  32, 

obtaining  in  this  way  two  dimensional  equations  where  all  flow  variables 
are  considered  as  functions  of  £*  and  £2  along  the  selected 
streamsurface. 

The  derivatives  of  the  flow  variables  with  respect  to  53  appearing  in 
the  obtained  "streamsurface"  equations,  are  considered  as  known 
quantities.  This  is  also  the  case  for  the  geometry  of  the  streamsurface 

S,  which  has  to  be  given  through  the  knowledge  of  the  normals  n  in  each 
point.  Since  the  flow  surface  will  only  be  known  once  the  whole  flow 
field  has  been  obtained,  the  present  decomposition  of  the  three 
dimensional  flow  field  is  clearly  an  iterative  procedure. 
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In  order  to  obtain  the  necessary  information  about  the  derivatives  of 
the  flow  variables  in  the  third  direction  53,  the  two  dimensional  flows 
have  to  be  solved  along  all  the  surfaces  of  the  same  family  53  - 
£3(£l,£2).  This  corresponds  to  the  different  surfaces  of  the  type 
AxB^jD,  of  figure  1.1.1  or  to  the  different  surfaces  Sj  on  figure 

1.1.2. 


On  the  other  hand,  the  information  which  is  necessary  in  order  to  define 
the  shape  of  these  streamsurfaces  can  only  be  obtained  through  the 
determination  of  the  flow  components  along  the  streamsurfaces  of  another 
family,  such  as  the  surfaces  of  the  type  A2B2C2D2  on  figure  1.1.1,  or 

the  S2  family  of  figure  1.1.2,  when  the  first  family  is  taken  to  be  S, 

surfaces . 


Therefore,  for  time  independent  flows,  the  complete  three  dimensional 
flow  field  will  be  reconstructed  through  the  iterative  solution  of  the 
two  dimensional  flows  on  two  families  of  intersecting  streamsurfaces. 

In  the  following  we  will  follow  an  approach  which  differs  from  the 
"algebraic"  method,  in  order  to  derive  the  two  dimensional 
"streamsurface"  equations.  We  will  express  directly  the  conservation 
laws  on  a  finite  volume  delimited  by  two  neighbouring  streamsurfaces, 
forming  a  streamsheet  of  thickness  b.  This  will  lead  to  a  formulation 
which  is  written  in  vector  form,  with  gradient  operators  along  the 

surface,  such  as  the  V  operator  of  equation  (1.1.17).  The  obtained 
formulation  will  therefore  be  valid  in  any  system  of  coordinates  and  the 
explicit,  algebraic  form  of  the  flow  equations  can  be  obtained  by 
expressing  the  gradient  operators  in  function  of  the  selected  coordinate 
system. 


1.1.2.  Streamsurface  Formulation  of  the  Mass  Conservation  Equation 


Referring  to  figure  1.1. A,  we  consider  the  domain  enclosed  between  two 
streamsurfaces,  forcing  a  streamsheet  of  thickness  B.  A  reference 
streamsurface  is  considered  as  formed  by  the  mid-points  of  the 
streamsheet  and  the  obtained  differential  equations  are  referred  to  this 
surface.  In  order  to  define  the  streamsurfaces,  a  curvilinear 
coordinate  system  is  introduced  with  the  direction  53  representing  the 
coordinate  which  is  to  be  eliminated  from  the  flow  equations.  The 
streamsurface  will  be  assumed  to  contain  the  and  £2  coordinates  and 
is  therefore  defined  as  a  (’  =  constant  surface.  We  will  select  5*  and 
S2  to  form  an  orthogonal  system  in  the  streamsurface  with  metric 
coefficients  hi  and  h2. 

a 

The  unit  vectors  along  the  £  axes  are  defined  by  e^  related  to  the  ea, 
if  the  coordinate  system  is  orthogonal,  by 


-► 

e 

a 


h 


e 

a 


h 

a 


-►a 

e 


-►a 

e 


no  summation  on  a 
a  -  1,2,3 
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The  normal  vector  n  is  equal  to  e*  according  to  equation  (1.1.5).  The 
physical  streamsheet  thickness  in  the  direction  of  the  coordinate  £*,  is 

measured  by  B. 

The  thickness  B  is  also  written  as  B  -  bhl(  where  hj  is  the  metric 
coefficient  associated  to  £*,  that  is  the  elementary  length  dlj  along 
this  axis  is  dlj-hjdS’.  The  quantity  b  is  then  a  thickness  measured  in 
units  of  the  variable  £*.  For  instance,  if  £’-0  the  angular  variable  in 
cylindrical  coordinates,  then  b  is  an  angular  thickness  and  B-br  where  r 
is  the  local  radial  coordinate. 

The  volume  di  formed  by  variations  d£ 1 ,  d£2  and  d£*  along  the  three 

£a-axes ,  is  measured  by  the  determinant  of  the  inverse  Jacobian  matrix 
of  the  transformation  from  a  cartesian  coordinate  system  to  the 

(X  1 

£  -system,  1/J.  That  is  dT  -  j  d£ldS2d£*  and  when  the  ^-system  is 
orthogonal,  1/J  ■  h,h2h,. 

Therefore,  the  elementary  volume  dft  formed  with  the  length  elements  dl,» 
hid5!»  dl2-  h2d£2  on  the  surface,  and  b  in  the  direction  £*,  see  figure 
1 . 1 .  JJ,  is  defined  by 

dfi  -  j  dC‘dt2  (1.1.19) 

If  the  metric  coefficients  ht  and  h2  are  introduced,  we  can  define  the 

streamsheet  thickness  B,  measured  in  the  direction  normal  to  the 
surface,  by 

^  -  Bh,ha  (1.1.20) 


When  the  coordinate  system  is  orthogonal,  B-B-bhj,  but  in  general  the 

ratio  B/b  is  equal  to  the  ratio  of  the  elementary  distance  dn  along  the 
normal  to  the  surface,  and  the  variation  d£*  along  the  direction. 
That  is 


B  dn 

b  ”  dTr 


(1.1.21) 


The  general  form  of  a  scalar  conservation  law  for  a  quantity  u,  is  given 
by 
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(1.1.22) 


a 

at 


pudfl  ♦  d» 


'fl 


?.d3  «  0 


For  instance  for  the  mass  conservation  law,  u*1  and  F*pv. 

This  Integral  formulation  is  applied  to  the  infinitesimal  volume  dfl  of 
figure  1.1, U,  The  side  faces  1  and  2,  formed  by  the  c*  and  C* 
directions  have  the  surfaoe  elements 


dS(,)  -  5  «*.!> 


(1 .1.23) 


If  the  coordinate  system  is  orthogonal  dS^  ^-Bhid^.e'-Bl^hjd^.e1,  The 
balanoe  of  fluxes  for  faces  1  and  2  is  given  by 

f  ?.ds(l) 


(Flb/J  -  (Flb/J  de2)^ 


1+2 


w 


(bF‘/J)  dCldC* 


(1.1.2H) 


(Bh^jFMdC'dC* 


where  Fl  is  the  projection  of  the  flux  vector  on  the  base  vector  cl, 


,  *r  -r  , 

Fl  -  F. c 1 


(1.1.25) 


and  the  overbar  on  the  partial  derivative  Indicates  a  derivative  along 
the  streamsurface .  This  implies  that  with  another  system  of 
coordinates,  whereby  the  £l,  axes  are  not  in  the  surface,  the 
derivatives  will  have  to  be  evaluated  by  relations  (1.1.12).  A  similar 
contribution  is  obtained  for  the  balance  of  fluxes  through  the  faces  3 
and  *1 ,  with 


dS( 2 ) 


(1.1.26) 


leading  to 
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|  F.dS(2) 

3+J4 


( F2b/J  dS1) 


S2+d£2 


(F2b/J  d£*)^2 


( F2b/J)d£*d£2 


(Bh,h2F2)d^1d^2 


with 


F2  =  F.  e  2 


■+( 2 )  -*• 

If  the  coordinate  system  is  orthogonal  dS  -  Bhid£l.e2- 

The  balance  of  fluxes  for  the  remaining  faces  5  and  6 
element 


->( ») 
dS^  ' 


hih2d^ld^2.es  «  —  d5*d£2 


z  h1h2d^1dC2.£3 


is  obtained  as  follows, 

j  F.dS(  )  -  (f3h1h2dc‘dt2)^,+b  -  (f’hihjdS'd^2)^ 
5+6 

=  b  |^r  (h,h2f3)  d5‘d£2 

■  b  W  (B  h»h«?-") 


since,  from  equation  (1.1.29) 


f3 


-  F.  e 3 


F. 

n 


■+  -»  n 

-  <F-">  5 


F3 


B 

b 


Grouping  all  the  terms,  one  obtains  after  dividing  by 
conservation  law  for  the  streamsheet  of  thickness  B 


9_ 

at 


( puBh  jh2 ) 


(Bh ,h2F* ) 


(Bh,h2F2) 


<F.,n  h, 


In  the  left  hand  side  one  recognizes  the  expression 


(1.1.27) 

(1.1.28) 

Bh ih2d£l .  e2 . 

,  with  surface 

(1.1.29) 

(1 .1 .30) 

(1 .1 .31 ) 

d£‘d£2 ,  the 

H  ,  (1.1.32) 

h2) 

of  the  two 
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dimensional  divergence  operator  in  the  coordinates  £*  and  £2  along  the 
surface.  The  conservation  equation  becomes,  for  time  independent 
surfaces,  applying  the  definition  (1.1.17) 


at  (puB)  + 


V  (BF)  =  -b  a, 


(F.1  ) 
n 


(1.1.33) 


where  the  overbar  on  the  gradient  operator  indicates  a  "streamsurface  " 
derivative. 

The  right  hand  side  term  is  a  reminder  of  the  fact  that  the  present 
description  is  not  strictly  two  dimensional,  but  describes  a  selective 
sheet  out  of  a  general  three  dimensional  flow. 

The  condition  for  the  considered  surface  to  be  a  streamsurface  is,  from 

■+  -¥ 

equation  (1.1.3),  v.n=0.  Consequently,  the  right  hand  side  of  equation 

■+  ■+ 

(1.1.33)  vanishes  for  the  mass  conservation  equation  (since  F=»pv)  which 
reduces  to  the  following  form 

|^(pB)  ♦  V  (pB  v)  -  0  0  -I-3*4) 


Considering  now  equation  (1.1. 3*0  as  a  starting  point,  one  can  apply 
this  equation  to  any  system  of  coordinates,  with  the  definitions 
(1.1.1 3)  and  (1.1.17)  of  the  streamsurface  derivatives . 


1.1.3.  Definition  of  the  Streamsheet  Thickness 

The  above  formulation  of  the  continuity  equation  along  the  streamsurface 
has  to  be  identical  locally  to  the  full  three  dimensional  form  of  the 
mass  conservation  equation.  This  condition  provides  a  relation  between 
the  streamsheet  thickness  B  and  the  three  dimensional  properties  of  the 
flow. 

Indeed,  applying  the  relations  (1.1.13)  and  (1.1.17)  to  the  general 
three  dimensional  form  of  the  continuity  equation, 

If  ♦  $  (pv)  =  0  (1.1.35) 


gives 


i£  ♦ 

3t 


?( pv) 


-pn.33v 


(1.1.36) 


Equation  (1.1.36)  should  be  identical  to  the  streamsheet  form  (1.1.34) 
of  the  mass  conservation  law.  Identifying  these  two  equations,  leads  to 
a  relation  defining  the  streamsheet  thickness  as  a  function  of  the  three 
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dimensionality  of  the  flow,  namely 


II  + 

3t 


(v.tf)B 


B  n. d3v 


(1.1.37) 


For  time  independent  flows,  the  time  derivative  vanishes  and  the  above 
relation  becomes 


5  (v.V)B 


n.3jv  -  -v.3sn 


(1.1.38) 


If  the  three  dimensional  flow  field  is  known  and  the  streamsurface 

normals  n  defined,  equation  (1.1.38)  will  determine  the  evolution  of  B 
along  the  selected  streamsurface.  Note  that  the  operator  on  the  left 
hand  side  is  the  convective  derivative  in  the  streamsurface.  This 
equation  expresses  that  the  variation  of  the  streamsheet  thickness  is 
defined  by  the  variation  of  the  velocity  field  in  the  third  direction. 
Therefore  if  the  flow  is  independent  of  this  coordinate,  the  right  hand 
side  of  equation  (1.1.38)  will  vanish  and  B  may  be  taken  as  constant  : 
the  flow  is  then  two  dimensional. 


Geometrical  interpretation 

The  condition  (1.1.38)  on  B  can  also  be  written  as  a  purely  geometrical 
relation  between  the  variations  of  B  and  the  normal  to  the 
streamsurface,  expressing  that  a  divergence  or  contraction  of  the 
streamsheet  necessarily  implies  that  the  direction  of  the  normals  to  the 
surface  change  during  a  displacement  accross  the  streamsheet,  as 
illustrated  on  figure  1.1.5. 

From  figure  1.1.5,  considering  a  displacement  dl  in  the  streamsheet,  it 
can  be  seen  that  the  divergence  angle  a  of  the  streamsheet  is  given  by 

tan  a  ’  \  HI  (1.1.39) 


and  that  the 

normal  unit 

vectors  rotate  by 

the  same  angle  in  the 

displacement 

from  P  to  Q. 

Hence 

1 nQ  ’  1  nP 

B  3r„  . 

'  2  Sn  '  ‘  l'“ 

(1  . 1 .140) 

where  d/dn  represents  the  derivative  in  the  normal  direction.  With  the 
approximation  tan  a  -  a  one  obtains 

1  dB  ♦  ^_n  (1  .1  .141  ) 

B  dl*  1  3n 
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Since  1^  is  an  arbitrary  direction  within  the  streamsurface ,  the  left 
hand  side  operator  is  equal  to  VB/B,  and  equation  (1.1.41)  becomes  with 
(1.1.21) 

1  in  _  b  9ln  (1.1.42) 

B  V  3n  B  3^ 

When  C  is  along  the  normal  dn=h3df;  ,  1  -  h3n  -  h3e  and 

|  VB  =  -  3,  (h3n)  -  -  3,n  -  jj  33h3  (1.1.43) 

Introducing  B»B»bh3,  in  this  relation,  leads  to  the  following  equation 
for  b 


b  5b  ‘  -3*"  -  h,  5h> 


(1.1.44) 


Equation  (1.1.44)  can  also  be  written  in  the  following  form,  which  is 
more  appropriate  for  explicit  computations, 


(v.V)b 

b 


a 


-v 


.  3  3n 

a 


n 

a 


.  3  3v 


(1.1.45) 


where  the  orthogonality  property  (1.1.4)  has  been  introduced  under  the 
form 


a 

v  .n 


(1.1.46) 


Example  E. 1 . 1 . 1 .  Cartesian  coordinates 

In  cartesian  coordinates  and  a  streamsurface  defined  by  z-z(x,y)  one  can 
eliminate  the  z-variations  along  this  streamsurface  in  the  continuity 
equation  (1.1.34).  One  has 

|j(PB)  .  |j  (pB  u)  *  |y  (PB  V)  -  0  (E.l.1.1) 


The  streamsurface  derivatives  are  defined  by 
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3_  _  3_ 
3x  3x 


(E.1 .1 .2) 


3_  a  3_  _  3_ 

3y  3y  ny  3z 


(E.1 .1 .3) 


with  the  normal  vector  n  given  by 


-  <!*>  r  -  <jj*>  t  *  r 

3x  x  3y  y  z 


(E.1  .1  .11) 


Equation  (1.1.38)  becomes 


1 

B 


(u.|-  B  +  v.|—  B) 
3x  3y 


n. 


-► 

3v 

3z 


(E. 1 . 1 .5) 


Observe  that  B~b,  since  all  metric  coefficients  are  equal  to  one. 

Example  E, 1 . 1 ,2.  Cylindrical  coordinates  -  SI  family 

If  the  streamsurf ace  SI  of  figures  1.1.1  or  1.1.2  is  considered  as  a 
surface  defined  by  r-r(0#z),  in  cylindrical  coordinates,  one  can 
identify  £*=0,  52=z  and  £3«r.  The  metric  coefficients  are  hx  -  r  , 

h2=h  j  =  1  . 


The  continuity  equation  (1.1. 3*0  becomes 


ft  (pB)  *  r  fe  (pBve)  *  r  b  <pBV> 


(E.1 .1 .6) 


where  the  velocity  components  vr*v2'v0  denote  the  physical  components. 
The  streamsurface  derivatives  are  defined  by 


*£  .  *£  -  n 

30  30  0  3r 


-  n 

3z  3z  z  3r 


(E.1 .1 .7) 


where  the  components  of  the  normal  vector  n  -  V  r(0,z)  are  given  by 


17 


n 


0 


8r(0  f z) 

86 


r.n 


0 


n 

z 


-8r(6,z ) 

8z  =  n  z 


n  -  1  -  n 
r  r 


(E. 1.1 .8) 


where  n  represents  the  projections  of  n  along  the  corresponding  unit 
vectors . 

Equation  (1.1.38)  for  the  streamsheet  thickness  B=b  becomes 
-  5  _  Ve  9B  A  §B  _  r  9  ,V  9  9  ,  (E.l.1.9) 

v.V  B  »  —  +  v  x—  =  B  [n.  —  ( — )  ♦  n  r—  v  +  n  tt~vJ 

—  r  96  z  9z  6  9r  r  r  9r  r  z  9r  z 


Example  E. 1 . 1 . 1 .  Cylindrical  coordinates  -  S2  family 

If  the  streamsurface  S2  of  figures  1 . 1 . 1  or  1.1.2  is  considered  as  a 
surface  defined  by  0=6(r,z)  in  cylindrical  coordinates,  one  can  identify 
El-z,  C2=r  and  C3=0.  The  metric  coefficients  are  h,=  1  ,  h2  =  1  ,  hs  « 

r.  The  continuity  equation  (1.1.34)  becomes,  with  B"b.r 


It  *  h  (pBv2>  *  I?  (cBvr>  ■ 0 


(E. 1 . 1 .10) 


The  streamsurface  derivatives  are  defined  by 


is  ,  1ft  -  n  1ft 

9z  9z  z  90 


(E.1 .1 .11) 


1ft  _  1ft  -  n  is 

9r  9r  r  90 


where  the  components  of  the  normal  vector  n 

8e(zfr) 

n  =  -  — r -  «*  n 

z  8z  z 

80(z  fr ) 

n  *  -  — 1 —  -  n 
r  8r  r 


V  0(z,r)  are  given  by 


(E. 1 . 1 . 12) 


ne  =  1 


r .  n. 


Equation  (1.1.45)  for  the  streamsheet  thickness  b  becomes, 


.2  9v  9v  n  9v 

nt.  8 b  8b  .  r  z  r  0  0-i 

v.Vb-v  t—  +  V  r —  •  b  [n  -r—  +  n  — r—  + - r— ] 

—  z  9z  r  9r  z  96  r  96  r  90 


(E.1. 1.13) 
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1.1.^.  Streamsurface  Formulation  for  the  Momentum  Conservation  Equation 


The  momentum  equation  is  written  in  integral  form  in  an  absolute  frame 
of  reference,  in  presence  of  external  forces  f 

=  -4  (  -► 

»  F.dS  »  pf  dfl  (1  .1 .147) 

s  e 


at  L 


pvdfi 


where 


F  *  pv  ^  v  +  pi 


T 


(1 .1 .H8) 


is  the  momentum  flux  tensor. 


Following  the  procedure  which  led  to  equation  (1.1.32)  in  the  previous 
section,  one  obtains, 

ft  +  hV  ^9l  ^'•^1  Bhih2)  +  a2  (F.e2  Bh  ,h2 )  ] 

(1.1. 149) 

=  p^-B  "  <f.?  h.h*) 

e  h  jh  2  n 


x  -4 

The  product  (F.1  )  reduces  to  the  contributions  from 

4  4 

shear  stresses  since  v.n=0  expressing  the  fact 
streamsurf ace .  Hence , 


the  pressure  and 
that  one  follows  a 


-F.  1 

n 


(1.1.50) 


The  derivative  terms  on  the  left  hand  side,  with  an  overbar  to  indicate 
that  they  are  to  be  considered  along  the  streamsurface ,  have  the  form  of 
the  divergence  of  a  tensor  in  the  two-dimensional  space  Hence 
the  quasi-three-dimensional  momentum  equation  becomes,  applying  equation 
(1.1.17) 

B  +  +  g  V[(pl  -  t )B]  -  pfe  +  |  3j(o.Tn)  (1-1-5D 


written  completely  in  the  two-dimensional  space  ^-S2. 

The  additional  term  in  the  right  hand  side  appears  as  an  extra  force 

4 

term  f  acting  on  the  "streamsurface"  flow  and  originating  from  the  three 
dimensionality  of  the  flow. 
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(1.1.52) 


pBf  »  -b3,  [ (pi  -  t).T  ]  »  b3j(o.t  ) 

n  n 

An  alternative  expression  for  the  additional  force  term  can  be  obtained, 
if  the  relation  (1.1. 42)  between  B  and  the  derivatives  of  n  is 
introduced  in  equation  (1.1.51). 

One  has 


V  (o  B)  ■  B(V.o)  ♦  (o.V)B  =  B(V.o)  -  o  b3(?  ) 

n 


(1.1.53) 


and  transferring  the  second  term  to  the  right  hand  side  of  equation 
(1.1.51),  leads  to  the  following  form  of  the  quasi-three  dimensional 

momentum  equation  applying  b1n=  B.n  as  derived  from  equation  (1.1.29) 


1  r3  , 


[■rr(pvB)  +  V  (pBv(x)v)  -  V.o]  =  +  n.3so  ♦  pf 


(1.1.54) 


B  ‘‘St 


The  first  term  on  the  right  hand  side  is  the  expression  of  an  additional 
force  fg,  obtained  from  the  relations 


pfB  =  pf  ♦  (°gV)B  =  n.3,3 


(1.1.55a) 


or  applying  the  relation  (1.1.21) 


n?  t  90 
pfB  =  'n  '  3^ 


(1.1.55b) 


Equation  (1.1.54)  takes  the  following  form 


B  *"ft^p^B^  +  +  ^(P1  “  t)  -  p?B  ♦  p f 


(1.1.56) 


Non-conservative  formulation 

Substracting  the  quasi-three  dimensional  continuity  equation  (1.1.34) 

multiplied  by  v,  from  the  momentum  equation  (1.1.51),  one  obtains  the 
non  conservative  formulation  for  the  quasi-three-dimensional  momentum 
equation, 


-► 

3v 

dt 


+ 


( v . V)  v 


"V(° B) 


f  +  f 


(1.1.57a) 
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or,  from  equation  (1.1.56) 


-»  -*• 


-► 


(1.1.57b) 


The  last  equation  has  the  usual  form  of  a  two-dimensional  momentum 

equation  in  the  subspace  of  the  streamsurface ,  with  however  the  addition 

-► 

of  the  distributed  force  fg  acting  on  the  flow.  Hence,  the  whole  effect 

of  the  three-dimensionality  on  the  non  conservative,  two  dimensional 

flow  description  is  contained  in  the  force  term  £3.  This  extra  force 
arises  from  the  projection,  on  the  normal  direction,  of  the  pressure  and 
shear  stress  gradients  in  the  third  coordinate  direction,  as  shown  by 
equation  (1.1.55).  Note  also,  that  the  conservative  formulation 
(1.1.56)  contains  also  the  influence  of  the  variation  of  the  streamtube 
thickness  B. 


In  absence  of  shear  stresses  and  in  a  steady  flow  configuration, 
equations  (1.1.55)  and  (1.1.3)  show  that  the  additional  force  fD  is 

D 

orthogonal  to  the  velocity  field  along  the  streamsurface  S.  One  has  in 
this  case 


^T 


l  iR 

n  3n 


(1.1.58) 


and 


P  fB.v  ^  0 


(1.1.59) 


This  is  an  interesting  relation  which  will  be  used  in  conjunction  with 
the  distributed  loss  model  in  the  field  of  turbomachinery  flow  models. 


The  results  obtained  in  this  section  remain  valid  for  other  forms  of  the 
momentum  equation,  in  particular  for  Crocco’s  form,  or  for  the  equations 
in  a  relative  system. 


Crocco's  form  of  the  momemtum  equation  becomes,  when  written  along  the 
streamsurface  S 


-► 

9v 

3t 


v  x  c 


TVs 


i  -  =  -►  -► 

VH+-V.x+fD+f 
p  Be 


(1.1.60) 


whit  the  absolute  vorticity  ;  defined  along  the  streamsurf ace  by 


V  x  v 


(1 .1 .61 ) 


where  all  the  derivatives  are  streamsurf ace  derivatives.  Note  that  by 


i 
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application  of  (1.1.13)  to  the  entropy  relation,  one  can  eliminate  the 
pressure  from  fDt  since 

D 


-  jj  3sp  =  n  (T3  3s  -  33h) 


(1.1.62) 


Example  E. 1 . 1 . 4  :  Cylindrical  coordinates 

Take  51-  0,  £2»  z,  r,  for  the  streamsurface  Si  defined  by  r  =  r(0,z) 
with  h i -r  ,  h2«  h3  =  1  . 

“►  -4 

The  additional  force  term  fg,  equation  (1.1.55),  becomes  with  n  defined 
by  equation  (E.1.1.8) 


Pf« 


[ 


_  !e 

3r 


3t 

3r 


]  n 


(E. 1 .1.14) 


For  the  S2  type  of  surface,  figure  1.1.2,  0=0(z,r),  one  obtains  the 

■+  4 

force  term  fg,  with  h3-r,  hi-h2-1  and  n  defined  by  equation  (E. 1.1. 12) 


[_  iE  +  ii]  ^ 

L  36  36J 


(E. 1 . 1 . 15) 


The  momentum  equations  are  obtained  in  cylindrical  coordinates  with  the 
addition  of  fg  and  with  the  derivatives  replaced  by  the  streamsurface 
derivatives  3,  with  3/30  set  to  zero. 

For  inviscid  flows,  one  obtains  for  instance 

pfg  -  -  n  (E.1 .1.16) 


and  the  radial  component  of  the  momentum  equation  becomes 


3v  3v  3v  ,  r 

_e:*v  _e  +  v  _n  ,  - 1  iE  +  f  +  f 

3t  r  3r  z  3z  p  3r  B,r  e,r 


(E. ) .1.17) 


1.1.5.  Streamsurface  Formulation  for  the  Energy  Conservation  Law 

The  energy  equation  can  be  treated  exactly  in  the  same  way  as  the  mass 
conservation  equation.  The  integral  form  of  the  energy  equation  is 
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(1.1.63) 


where  Q  is  a  source  term  and  the  energy  flux  is  defined  by 

F  =  pvH  -  i.v  -  kfo 


(1.1.64) 


Following  the  same  steps  as  in  section  1.1.1.,  equation  (1.1.63)  becomes 


&  (pBE> 


+  V  [ ( pvH  -  t  , 


kVT) B]  =  QB  -  b9  3 ( F. Tn ) 


(1.1.65) 


where 


F.  1 

n 


-(i.v  +  k^T) . T 

n 


(1 .1.66) 


since  v.n=0. 

Applying  the  continuity  equation  leads  to  the  alternative  form 
P  [!r  +  (^)H]  =  |f  +  g  V  [(i.v  ♦  kVT)B]  ♦  Q  -  |  3,(F.1n)  (1*1,67) 

or  introducing  equation  (1.1. ^2) 


-r  -r 

♦  ( v .  V)H]  =  I?  ♦  V  (i.v  +  k^T )  ♦  Q  -  n.  9  3F 


(1.1.68) 


9t 


9t 


Hence,  the  confinement  of  the  energy  balance  to  the  streamsurf ace  S, 
leads  to  a  conservation  equation  for  the  energy  fluxes  within  the 
streamsheet  of  thickness  B=bh3,  whereby  an  additional  heat  source  term 

appears,  due  to  the  three  dimensionality  of  the  flow,  given  by  equation 
(1.1.69). 


Qg  =  -  n. 83F  -  n.33(kVT  +  t.v) 


n  3n 


(1.1.69) 


These  additional  sources  result  from  the  normal  components  of  the 
gradients,  or  the  normal  direction  of  the  heat  conduction  fluxes  and  the 
work  of  the  viscous  shear  stresses,  in  the  third  direction. 

In  absence  of  shear  stresses  and  heat  conduction  effects,  the  quasi 
three  dimensional  form  of  the  energy  equation  reduces  to,  with  the 
source  term  Q~0, 
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where 


(1.1.71 ) 


is  the  total  material  derivative  along  the  streamsurface ,  and  is 
identical  to  its  two-dimensional  restriction. 


1.1.6.  Quasi-Three  Dimensional  Interaction  Procedure 

Summarizing  this  approach,  we  have  seen  that  the  flow  behaviour  can  be 
followed  along  a  streamsurface  £3=£ 3 ( £ 1 , )  and  the  conservation 
equations  for  mass,  momentum  and  energy  can  be  written  in  a 
two-dimensional  form  in  function  of  the  space  variables  El  and  E2  with 
space  derivatives  taken  along  the  streamsurface.  The  influence  of  the 
three-dimensionality  of  the  flow  appears  in  all  three  conservation 
equations  through 

i)  a  variable  streamtube  thickness  B 

ii)  through  an  additional  force  term  whose  non-viscous  contribution  is 
the  pressure  gradient  in  the  E3  direction  acting  in  the  direction 
normal  to  the  surface,  in  the  momentum  equation  and 

iii)  for  the  energy  conservation  equation,  through  an  energy  source 
term  whose  non-viscous  contribution  is  a  normal  flux  arising  from 
the  thermal  conductivity  in  the  E3  direction. 

In  order  to  solve  this  two-dimensional  system  of  equations,  the 
extra  terms  just  mentioned  have  to  be  known,  as  well  as  the  shape 
of  the  streamsurface.  These  informations,  known  only  if  the 
complete  three-dimensional  flow  description  is  available,  have 
therefore  to  be  obtained  in  an  iterative  procedure. 

In  practical  applications  to  internal  flows,  the  present  approach 
is  best  suited  for  steady  state  flow  situations,  since  in  this 
case,  one  can  completely  define  the  velocity  field  through  two 
families  of  streamf unctions  ip i  and  vz.  Tih  (1979). 

"►  -► 

The  steady  state  continuity  equation  V(pv)  =  0  is  satisfied  by  the 
existence  of  two  functions  tp t  and  such  that 


pv 


(1.1.72) 


By  definition,  one  has 


(v.$)Y,  =  0 


(1.1.73) 


(v.$H2  -  0 


(1.1.74) 


which  shows  that  the  velocity  vectors  are  tangent  to  the  surfaces 


( £* ,52 »53 )  *  constant  and  fZ2  ,Z3 )  -  constant  and  hence  that 

the  intersection  of  any  of  these  two  surfaces  is  a  streamline. 
Therefore,  these  surfaces  form  two  families  of  streamsurfaces 
generated  by  the  various  streamlines  issued  from  the  fluid 
particles  situated  in  the  inlet  surface  of  the  flow  domain.  In 
figure  1.1.1,  the  surfaces  AiBiCjDj  (or  Sx  in  figure  1.1.2)  and 

A2B2C2D2  (or  S2  in  figure  1.1.2)  are  typical  representations  of 

surfaces  of  constant  Vi  and  values. 


The  complete  three-dimensional  flow  description  is  therefore 
reconstructed  by  writing  the  two-dimensional  form  of  the  equations 
of  motion  along  two  families  of  intersecting  surfaces  such  as 

and  \p2.  The  information  necessary  on  one  family  of  surfaces  is  to 

be  obtained  from  the  flow  behaviour  on  the  other  family.  For 
instance,  if  the  surface  S  considered  above  is  a  member  of  Vi 

family,  the  shape  of  this  surface  as  well  as  the  streamtube 
thickness  B  and  the  additional  terms  which  are  all  determined  by 
the  variations  of  certain  flow  variables  in  the  direction,  will 
be  obtained  by  the  solution  of  the  same  two-dimensional  form  of  the 
equations  written  along  the  other  ^2  family  of  surfaces,  which 

contains,  by  definition  the  variations  in  the  £s  direction. 

For  this  latter  family  of  two-dimensional  flows,  the  knowledge  of 
the  additional  terms  can  only  be  obtained  from  the  family  of 

flows  and  this  sets  up  the  iterative  procedure. 

In  practical  applications,  the  computation  will  start  by  assuming  a 
reasonable  shape  for  the  first  family  of  surfaces  (vi)  and  by 

neglecting  the  contributions  of  the  additional  terms.  After  the 
solutions  along  this  first  family  are  obtained,  the  flow  is 
calculated  along  the  second  (\^2)  family  of  surfaces  with  an  input 

from  the  first  family  in  order  to  determine  the  three-dimensional 
contributions  to  the  equations  of  motions.  After  this  step,  a  new 
approximation  to  the  \j;x  flows  can  be  obtained  and  the  whole 

computation  is  continued  solving  alternatively  the  flow  along  each 
set  of  streamsurf aces  until  a  chosen  convergence  criterion  is 
satisfied . 

This  calculation  procedure,  in  a  simplified  form,  is  widely  used  in 
the  field  of  single  and  (or)  multistage  turbomachines  and  will  be 
described  more  in  details  with  regard  to  computational  details  and 
additional  approximations  in  part  2  of  this  report  dealing  with  the 
internal  turbomachinery  flows. 
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1.2.  THE  QUASI-THREE-DIMENSIONAL  REPRESENTATION  (Q3D)  FOR  INTERNAL  FLOWS  : 

THE  AVERAGING  PROCEDURE 


As  discussed  in  the  introduction,  an  alternative  to  the 
streamsurface  approach  in  reducing  the  complexity  of  the  three 
dimensionality  of  the  flow,  can  be  obtained  by  taking  the  average 
of  the  conservation  equations  over  one  of  the  spatial  directions. 
This  procedure  will  lead  to  exact  equations  for  the  flow 
conservation  laws  in  the  remaining  space  variables,  reducing  in 
this  way  the  flow  description  to  a  lower  number  of  space  variables. 

Since  this  averaging  technique  requires  limits  of  integrations  for 
the  corresponding  variable  the  method  takes  its  full  potentiality 
within  the  framework  of  internal  flows,  such  as  channel  or 
turbomachinery  flows. 

Let  us  consider  therefore  a  channel  limiting  the  flow  region  by 
upper  and  lower  surfaces,  as  well  as  side  surfaces,  these  surfaces 
being  either  material,  solid  walls  or  fluid  boundaries,  figure 
1.2.1.  Average  values  of  any  flow  variable  A  over  say,  the  5 3 
direction,  can  be  defined  through 


b, 

»  -  bu’ls’.t)  I  *  (c-.c’.c’.t) 

b  o 


(1.2.1) 


The  quantities  ^ 3 «b i(Cl »Ca.t)  and  S’-boCS1 ,52,t)  define  the  limits 
of  the  integration  region  corresponding  to  the  surfaces  A^CjC, 
and  B!B2D2Dj ,  with 


b  -  b,  -  b0  (1 .2.2) 


(b»A,B,  on  figure  1.2.1).  These  boundary  surfaces  are,  in  general, 

functions  of  the  other  coordinates  and  eventually  of  time.  Note 
that  b  is  the  width  of  the  integration  region  expressed  in  units  of 

Relations  for  the  averaged  derivatives  are  obtained  as  follows  : 
from  the  general  properties  of  the  derivative  of  an  integral  one 
has,  the  subscript  a  indicating  either  time  t  or  the  coordinates  51 
or  52, 
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9  A 
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l 

b 


bi 

■ 

9  A  df;3 
a 

b0 


b, 


J  9a 


A  d£ 


3  _ 


-  [A  (el»ea»blit)3ob1  -  A  (tl.5a»b0.t)  9ab0] 


(1  .2.3) 


1 

b 


9 

a 


(bA) 


1 

b 


(A  1 9^b ! 


A09ab0) 


where  A!  and  A0  are  the  values  taken  by  the  variable  A  on  the  upper 
and  lower  limits  of  the  integration  domain. 

Similarly,  for  a=3»  one  obtains 

97a  =  J  (A,  -  A0)  (1 .2.14) 


On  the  other  hand,  the  average  of  a  product  of  two  quantities 
requires  the  introduction  of  the  deviations  from  the  average  value 
following 


A  =  A  +  A' 


(1.2.5) 


with,  by  definition 


A'  =  0 


(1.2.6) 


Hence,  one  obtains  for  the  average  of  a  product  of  two  quantities  A 
and  B  and  for  their  derivative,  the  following  relations 

AB  »  A.B  +  A'B'  (1.2.7) 

9  AB  =  i  9  (b.A.B)  +  -■  9  (b.A'B' )  -  l  (A,B,9  b,  -  A0B09  b0)  (1.2.8) 

a  ba  ba  b  a  a 

where  a  represents  the  coordinates  £2  or  time. 

The  boundary  surfaces  of  the  integration  region  will  be  defined  by 
S,  -  C  -  b,  Ul  -  0  (1.2.9) 

and 

So  ;  (3  -  bo  Ul,52,  t)  =0  (1  .2.10) 

for  the  upper  and  lower  limits  respectively. 
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An  important  simplification  of  the  averaging  process  occurs  when 
the  boundary  surfaces  can  be  assumed  to  be  streamsurfaces ,  that  is 
surfaces  which  do  not  allow  any  cross-flow.  This  assumption  will 
be  satisfied  for  the  solid  surfaces  bounding  a  duct  or  channel 
flow,  but  is  also  valid  in  free  surface  problems  where  one  of  the 
boundaries,  represents  the  time  and  space  dependent  free  surface  of 
an  open  channel  fluid  flow,  as  occurs  in  river  and  sea  hydraulics 
or  tidal  flows. 

Expressing  that  the  velocity  vector  lies  at  each  instant  t  in  these 
surfaces,  that  is 


d_ 

dt 


S0 


3So  * 

T7 —  ♦  (v.V)So  =  0 


(1.2.11) 


along  the  surface  So,  one  obtains  with  (1.2.9)  and  (1.2.10)  the 
relation 


3b0 

3t” 


■4  -4 

♦  v.n0 


0 


where  n0 


=vs, 


is  the  normal  to  the  considered  surface. 


(1 .2.12) 


4 

The  components  of  n0  along  the  unit  vectors  of  an  orthogonal 
coordinate  system  are  defined  by,  see  equation  (1.1.6) 


(Z±  3b°  3b°  1_ 
n,  as*  ’  h2  aiF  ’  hj 


(1.2.13) 


Similarly,  one  has  for  the  second  boundary 


3b, 

at- 


■* 

*  v.n, 


0 


with 


-1  3b, 

(hT  W 


-1  3b, 

hi  W 


(1 .2.  Ill) 


(1.2.15) 


With  these  relations,  the  averaging  rules  can  be  written,  grouping 
(1.2.3)  and  (1.2.4)  as  the  average  of  the  gradient  of  a  quantity  A. 


$A  »  ■[  $  (bA)  +  j-  [A.n] 
b  b 


(1.2.16) 


where  the  square  brackets  indicate  the  difference  over  the  two 
boundaries 
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(1.2.17) 


[An]  -  A^,  -  A0n0 


The  gradient  operator  is  restricted  to  the  two-dimensional  space 
,  since 


(1.2.18) 


for  any  averaged  quantity  A. 

The  same  relations  apply  to  the  divergence  of  a  vector  F  or  a 

a 

tensor  quantity  xf  that  is, 

$.F  =  l  $  (Fb)  +  £  [F.n]  (1.2.19) 

b  b 

and 


n  -  j  t;.S] 

b  a  b 


(1.2.20) 


where  is  the  averaged  value  of  t.  These  averaging  rules  express 

that  the  averaged  gradients  of  a  flow  quantity  over  a  given 
direction  are  equal  to  the  gradients  in  the  (two)  remaining 
coordinates  of  the  averaged  quantity  weighted  by  the  width  of  the 
integration  region,  plus  a  contribution  expressing  the  balance  of 
the  normal  components  of  the  considered  quantities  over  the  two 
boundaries . 

It  is  essential  to  observe  that  the  explicit  form  of  the  divergence 
operator  acting  on  the  averaged  quantities,  has  to  be  written  in 
its  three-dimensional  form,  where  all  derivatives  with  respect  to 
5 3  are  set  to  zero  afterwards. 

These  two  relations,  written  in  operator  form,  are  valid  on  an 
arbitrary  system  of  coordinates  provided  the  metric  coefficients 
are  Independent  of  the  integration  variable.  When  this  is  not  the 
case,  the  divergence  operators  have  to  be  written  first  in  the 
system  of  coordinates  considered  and  then  integrated  over  the 
variable  5*.  Averaged  values  of  the  metric  coefficients  will  then 
appear  in  the  explicit  formulation.  For  instance,  if  the 
integration  variable  is  the  radial  coordinate  in  a  cylindrical  or 
spherical  coordinate  representation,  an  average  radius  will  appear 
in  the  equations. 

Hence,  in  orthogonal  curvilinear  coordinates  one  has 


$.(Fb) 


1 


hih*hj 


[9,  (hjhjhjb  F1)  +  92  (h,h2h,b  FJ ) ] 


(1 .2.21 ) 


This  has  to  be  compared  with  the  streamsurface  divergence  operator 
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V  as  defined  by  equation  (1.1.29). 

The  explicit  occurence  of  hs  in  the  above  equation,  explains  why 

the  quantity  b  appears  in  the  averaged  equations  while  the  quantity 
B-bhs  appears  in  the  streamsurface  equations  of  the  previous 

section. 

1.2.1.  The  averaged  form  of  the  continuity  equation 

9p  i  -► 

Applying  the  above  rules  to  the  continuity  equation  ^  +  vpv  »  0 

leads  to 

(pb)  +  $  (pvb)  +  [pv.n]  +  Otb0  -  3^,)  -  0  0-2.22) 


Expressing  that  the  two  bounding  surfaces  are  streamsurfaces , 
equations  (1.2.12)  and  (1.2.14),  the  averaged  continuity  equation 
becomes 


(pb)  +  $  (pvb)  =  0 


(1.2.23) 


Applying  equation  (1.2.7)  to  the  averaged  product  pv,  leads  to 

(pb)  +  $  (p  v  b)  -  -  $(p'v'b)  (1.2.24) 

If  an  integrating  factor  p  is  defined  such  that 


(p'v'b)  -  -  p(v.V)  p 


(1.2.25) 


equation  (1.2.24)  becomes 


+  $  (p  v  B)  -  0 


(1.2.26) 


where 


B  -  b. p 


(1.2.27) 


Comparing  this  with  equation  (1.1.34),  a  two-dimensional  mass 
conservation  equation  is  obtained  where  the  B-factor  represents  a 
varying  streamtube  thickness  which  is  equal  to  the  product  of  the 
geometrical  thickness  of  the  fluid  layer  over  which  the  averaging 
is  performed,  multiplied  by  a  correction  factor  which  represents 


-  30  - 


the  effect  of  the  mass  flux  fluctuations  over  the  flow  domain. 
Hence,  this  B  factor  is  dependent  on  the  three-dimensional  flow 
properties . 

It  is  to  be  noticed  that  for  incompressible  flow,  there  are  no 
density  fluctuations  and  hence  that  B»b. 

The  present  space  averaging  technique  has  strong  resemblances  with 
the  time  averaging  procedure  applied  in  the  study  of  turbulent 
flows.  An  elegant  way  of  accounting  for  the  effects  of  the 
turbulent  density  fluctuations  on  the  time  averaged  flow  is  to 
define  density  averaged  quantities  instead  of  the  usual  time 
averaged  variables.  A  similar  procedure  can  be  applied  in  the 
present  context  where  the  compressibility  effects  included  in  B  and 
described  by  the  terms  in  equation  (1.2.25)  can  be  absorbed  through 
the  definition  of  the  following  density-average  of  a  variable  A, 
see  for  instance  Hirsch  and  Warzee  (1979) 


fbl 

b0 


PA  d£3 


(1 .2.28) 


Denoting,  in  the  following,  the  ordinary  average  as  defined  by 
equation  (1.2.1)  by  an  overbar  and  the  density  averaged  values  by  a 
tilda,  one  has 


A 


p  pb  ;b0 


PA  dC3 


(1 .2.29) 


The  fluctuation  of  the  variable  A  with  respect  to  its  usual  average 
A  is  indicated  by  A'  while  A"  will  indicate  the  fluctuations  with 
respect  to  the  density  average  A.  Therefore,  the  following 
relations  hold, 


A  =  A  +  A '  »  A  +  A" 


(1 .2.30) 


with 


A'  =  0 
A"  #  0 


but 


pA"  =  0 

One  has  also,  from  the  above  definitions 

~  —  o  '  A  • 

(A  -  A)  -  -  A"  = 

P 


(1.2.3D 


(1.2.32) 
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and  for  the  product  of  two  quantities  A,  B 


pAB  -  p  A  B  +  pA"B"  '  JJ' 

while  the  above  relations  for  the  average  of  a  derivative  remain 
unchanged. 

Applying  these  various  properties  to  the  averaged  mass  conservation 
law,  equation  (1.2.23)  leads  to  the  equation 


3F  (pb) 


V  (p 


V  b) 


(1 .2.3*0 


writing  V  for  the  averaged  velocity  vector  v. 

This  exact  relation  has  the  advantage  that  the  flow  dependent 
streamtube  thickness  B  in  equation  (1.2.26)  reduces  here,  in  all 
cases,  to  the  geometrical  thickness  b  of  the  flow  region  in  the  5* 
direction,  since  the  influence  of  the  density  fluctuations  are 

-¥ 

included  in  the  definition  of  V.  Therefore,  if  the  surfaces  b0  and 

b,  are  not  flow  dependent,  the  streamtube  thickness  of  the 

considered  fluid  layer  b»bi-b0  is  a  known,  geometrical  quantity. 

This  represents  an  important  conceptual  simplification  of  the 
averaged  compressible  flow  equations  compared  to  the  previous  form 
of  equation  (1.2.26). 

Comparing  equations  (1.2. 3*1)  with  the  form  of  the  continuity 
equation  (1.1.3**)  obtained  following  the  streamsurface  approach, 
one  notices  that  not  with  standing  the  close  formal  resemblance, 
the  former  equation  describes  the  two-dimensional  evolution  of  the 
averaged  flow  considered  in  the  e;‘-e;4  coordinate  system  while  the 
latter,  equation  (1.1. 3*0,  describes  the  evolution  of  the  local 
flow  variables  followed  along  a  given  streamsurface  of  the 
three-dimensional  flow. 


1.2.2.  The  Averaged  Momentum  Conservation  Equation 

Considering  the  differential  form  of  the  momentum  equation,  in  an 
absolute  frame  of  reference,  eventually  as  a  Reynolds  averaged 
equation  for  turbulent  flows 

(pv)  +*V  (pv(x)v  +  pi  -~t )  -  p£  (1.2.35) 


the  application  of  the  averaged  rules  (1.2.16)  to  (1.2.20)  leads  to 
the  following  averaged  momentum  equation 


!^(pvb)+$(  pv  ©vb  +  pb  -  tab)  -  P?e*b  -  [pv  (v.n  -  |^)]-[pn]  +  [t  .n]  ^ 1  * 2‘  ^ 


writing  t  for  the  averaged  viscous  shear  stress  tensor,  plus 

3 
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eventually  the  Reynolds  stress  tensor. 


The  first  term  under  brackets  in  the  right  hand  side  vanishes  when 
the  streamsurface  assumptions  equations,  (1.2.12)  and  (1. 2.1*0,  are 
satisfied.  Otherwise  they  contribute  to  an  impulse  flux  through 
the  boundary  surfaces.  The  last  two  terms  represent  the  components 
of  an  additional  force  term  which  appears  as  an  outcome  of  the 
averaging  procedure  in  a  way  that  completely  parallels  the  way  the 
additional  force  appears  in  equation  (1.1.52)  derived  in  the 
streamsurface  approach.  The  force  term,  defined  by 

P  f  *  -  g  [o.n]  =  -g  [ (pi  -  t  )  .n]  (1.2.37) 


results  from  the  difference  in  the  normal  components  of  the 
pressure  and  shear  stresses  acting  on  both  limiting  surfaces  bj  and 

b0. 


The  close  resemblance  between  the  two  expressions  of  the  extra 
force  term  equation  (1.1.52)  and  (1.2.37)  is  obvious,  so  more  that 
they  become  identical  if  the  averaging  is  performed  over  an 
infinitesimal  distance  d^3»b,  with  metric  coefficients  independent 
of  £3.  Hence,  equation  (1.2.36)  becomes 


ft(pvb) 


$  ( pv  (x)  v  b) 


(o.b)  + 


pfeb 


pfb 


(1.2.38) 


writing  og  for  the  averaged  shear  stress  tensor  o. 


If  the  averaging  rule  of  products,  equation  (1.2.33)  is  introduced 
with  density  weighted  averages,  one  obtains  the  following  equation 

-► 

for  the  density  averaged  velocity  vector  V*v, 


f^pVb) 


$  (pbV(x)V)  -  +  $  (o  ,b) 


pfeb 


pfb  -  $  (  pv"  (x)  v"b) 


(1.2.39) 


or  alternatively,  taking  into  account  the  averaged  continuity 
equation  ( 1 .2. 3*0  » 


f-  $  *  (V.$)  V  -  -  —  $  (o  .b)  ♦  f  ♦  f  -  —  ^  ( pv"  fx)  v"b)  (1.2.110) 

at  .  •  a  e  - 

bp  bp 


This  last  equation  can  be  made  formally  identical  to  the  standard 

two-dimensional  momentum  equation,  by  introducing  the  force  f 
■* 

instead  of  f,  with 


pb  f  -  -(pi  -  t  ).7b  ♦  pbf 
t>  a 


(1 .2.111  ) 
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leading  to 


4  ■  ■ 

♦  (V.$)  V  =  $  (pi  -  Ta)  ♦  f  +  ?  -  $  (pv"0v"  b)  (1.2.42) 

p  bp 


With  the  exception  of  the  last  term,  this  equation  and  equation 
(1.1.58)  obtained  by  following  a  streamsurface ,  have  very  similar 
forms,  although  the  interpretations  of  each  term  differ  by  the  fact 
that  here  local  derivatives  of  averaged  quantities  appear  while  in 
the  former  approach  one  has  to  do  with  streamsurface  derivatives  of 
the  local  flow  quantities.  In  both  cases  however,  the 
quasi-three-dimensional  approximation  leads  to  the  introduction  of 

4 

variable  streamtube  thickness  (B  or  b)  and  to  additional  forces  f 
whose  non-viscous  components  act  in  the  direction  normal  to  the 
considered  surfaces. 

In  the  present  approach  an  additional  term  appears  on  the  right 
hand  side  of  the  momentum  equation.  This  term,  which  has  the  same 
structure  as  the  Reynolds  stress  tensor  in  the  turbulent  averaged 
Navier-Stokes  equations,  contains  the  influence  of  the 
three-dimensionality  of  the  flow  on  its  averaged  representation. 

It  is  to  be  kept  in  mind  that  the  v"  components  of  the  velocity 
field  represent  the  deviations  from  flow  uniformity  in  the  5 5 
direction  and  can  therefore  be  considered  as  large  scale  velocity 
fluctuations  around  the  averaged  flow. 


Hence,  the  products  pv"  x  v"  are  an  expression  of  the  large  scale 
momentum  transport  around  the  averaged  flow  and  its  gradients  act 
as  a  diffusion  effect  on  the  resulting  averaged  momentum.  If  the 
flow  is  mildy  non-uniform  in  the  averaging  direction  these  terms 
may  become  neglectable  with  respect  to  the  other  terms.  If  not, 
they  can  only  be  evaluated  through  a  knowledge  of  the  flow  in  a 
family  of  surfaces  containing  the  ^“direction,  in  the  following 

„S 

we  will  designate  these  terms  as  secondary  stresses  x 


-S 

T 


(1.2.43) 


1.2.3.  The  Averaged  Energy  Conservation  Equation 

The  energy  equation  can  be  derived  in  the  same  way  as  described 
above  using  the  local  form,  corresponding  to  equations  (1.1.62)  to 
(1.1.64) 


+  $  ( pvH  -  G)  =  Q 


(1.2.44) 


where  G  is  the  flux  of  the  shear  stresses  work  and  heat  conduction 


+  jt  rs  4 

G  -  kVT  +  t.v 


(1.2.45) 


The  averaging  procedure  leads  to  the  following  form  of  the  averaged 
energy  equation 


|^(pE  b)  ♦  $[(pvH  -  G)b]  -  Q.b  -  [p  |g]  ♦  [G.n]  -  [pH  (v.n  -  ||)]  (1-2‘J<6) 


The  last  term  vanishes  when  the  boundaries  are  streamsurfaces , 
following  equations  (1.2.12)  and  (1.2.14). 


Two  additional  contributions  to  the  averaged 
in  the  right  hand  side  is  source  terms.  The 

energy  balance  appear 
first  one 

r  3b,  3b"  Sb' 

[p  3tJ  *  Po.  3t  Pi-  3t 

(1.2.47) 

is  the  work  of  the  pressure  forces  against  the  moving 
while  the  second  term 

boundaries , 

[G.n]  -  [ (kVT  ♦  t.v)  .n] 

(1.2.48) 

is  the  transfer  of  heat  flux  and  work  of  the  viscous  stresses 
accross  the  boundaries  of  the  integration  region.  These  two  terms 
will  be  grouped  in  the  source  term 

b  Qb  -  [G.n]  -  [p  |£]  (1.2.49) 


When  the  averaged  products  on  the  left  hand  side  are  worked  out, 
one  obtains 


|-  (pEb)  +  $[(pVH  -  G)b]  =  Q.b  ♦  Qb.b  -  $  (bpv"H") 
or  by  application  of  the  averaged  continuity  equation 
p  [~  ♦  (V.$)  H]  =  g  Ig-(pb)  ♦  g  $  (Gb)  +  Q  +  Qb  “  l  $  (b  pv"H" ) 


This  form  of  the  averaged  energy  conservation  equation  parallels 
completely  equation  (1.1.68)  obtained  by  the  streamsurface  method 
with  the  difference  in  Interpretation  already  discussed  above  for 
the  momentum  equation.  An  additional  energy  flux  term  appears, 

pv"H",  which  represents  the  effect  of  the  total  enthalpy 
fluctuations  as  convected  by  the  large  scale  velocity  fluctuations. 

Making  use  of  the  equations  (1.2.12)  and  (1.2.14),  the  additional 

source  term  can  be  written  as  follows 
b 
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or 


-  [(t.v).n 


P  n.v]  +  g  [k  $T.n] 


g  [(T.n).v 


p  n.v]  +  -  [k  ^T.n] 


(1.2.52) 


In  general  circumstances,  the  first  term  of  equation  (1.2.52)  shows 

that  the  additional  heat  source  Q.  arises  from  the  difference  in 

b 

the  work  of  the  shear  stresses  and  pressure  forces  along  the  two 
limit  surfaces  b0  and  bi,  in  addition  to  the  balance  of  heat 

conduction  through  these  two  surfaces. 

An  interesting  particular  situation  occurs  when  both  end-surfaces 

b„  and  b,  have  the  same  velocity,  say  v^.  In  this  case,  the  first 

“►  ♦ 

term  of  represents  the  work  of  the  additional  force  f,  equation 
(1.2.37),  acting  on  the  moving  boundary  surfaces, 

Qw  *  w  Ct-n  -  p.n].v  ♦  -  [k  VT.n] 

D  D  S  D 


or 


pf.v  ♦  i  [k  VT.n] 
s  u 


(1.2.53) 


A  typical  situation  of  this  kind  occurs  in  the  passage  between  two 
blades  of  a  rotating  turbomachinery  blade  row. 

The  averaged  energy  equation  can  be  put  under  different  forms,  when 
the  stagnation  enthalpy  fluctuations  H"  are  separated  into  the 
static  enthalpy  fluctuations  h"  and  the  contributions  from  the 
kinetic  energy. 

The  total  enthalpy  fluctuations  can  be  written  as 

H"  =  H  -  H  =  h"  ♦  v".V  +  —  " 


(1 .2.5^) 


where  the  average  kinetic  energy  of  the  velocity  fluctuations  is 
introduced , 


P  k 


(1 .2.55) 


Referring  to  equation  (1.2.50),  the  last  term  of  the  right  hand 
side  becomes  after  introduction  of  equations  (1.2. 5*0  and  (1.2.H3) 
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pv"H" 


pv"h" 


(1.2.56) 


-S 

T 


+ 


where  t  designates  the  secondary  stresses  (1.2.43) 

The  energy  equation  (1.2.50)  becomes 

|^(pbE)  ♦  $  ( pVHb)  $  [Gb  -  pv"  (h"  +  ^-)  b]  ♦  $  (bViS)  +  b  (Q  +  Qfe) 

(1.2.57) 

-  -  JJQ  -  - 

=*  $(v.t  +  V.TS  -  pv"h"  +  £-2  h  -  pv"k"  )b  +  b(Q  +  Q  ) 

r 


Note  the  appearance  of  the  energy  due  to  the  work  of  the  secondary 

— S  ^  s . 

stresses  t  against  the  averaged  flow  field  V,  V(V.t  )>  in  the 

right  hand  side  as  well  as  the  transport  of  total  energy  of  the 
fluctuations  (hM  +  v,,2/2)  by  the  large  scale  fluctuating  flow  field 

V". 


If  the  deviations  from  the  average  velocity  are  small,  this  last 

contribution,  being  of  third  order  in  v"  might  be  neglectable 

=*S 

compared  to  the  work  of  the  secondary  stresses  t  against  the  mean 
flow  V.  In  addition,  if  the  work  of  the  fluctuating  viscous  (and 

turbulent)  shear  stresses  is  neglected  v"t  •  0,  then  the  averaged 
energy  equation  (1.2.57)  simplifies  to 


^(pbE) 


$(pVHb)  »  $[(V.t  -  pv"h" 


H!!e 

p 

r 


^T)b]  +  b(Q 


V 


(1.2.58) 


where  the  total  shear  stress  appears  as  the  sum  of  the  viscous 
(plus  turbulent)  and  secondary  stresses 

-  t  +  (1 .2.59) 


One  can  also  introduce  the  total  energy  of  the  averaged  flow 


H  =  h  + 


pv" .v"  “  - 

H  -  .-v  =  H  -  k 


(1 .2.60) 


2p 


into  equation  (1.2.57),  in  order  to  obtain  an  equation  for  H. 
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1.2.^.  Crocco's  form  for  the  averaged  momentum  equation 


If  the  entropy  relation  is  averaged  over  the  width  b  of  the  channel,  one 
obtains  from  equations  (1.2.7),  (1.2.16)  and  (1.2.33) 


)T$s  «  p^h  - 


(1 .2.61 ) 


or 


pT.^s  -  p.$h  -  $p  -  ( pT) " ( $s )"  +  p'($h)' 


(1.2.62) 


leading  to 


T  $(bs)  -  V(bh )  -  :  $(bp)  +  [ (h  -  £  -  Ts)n]  -  z  (pT)'(Vs)'  +  p ' ( Vh) ' ( 1 . 2. 63) 

P  p  p 


or,  grouping  the  three  last  terms  into  a  force  term  bfg , 


T  $(bs)  -  $(bh)  -  ~  $(bp )  +  b^c 

P 


(1 .2.6*0 


Eliminating  the  pressure  term  with  the  help  of  equation  (1.2. 6*0,  one 
obtains  Crocco's  form  of  the  momentum  equation  for  the  averaged  flow, 


SV  ■*  •+  ~+~  ■**  1  •+  =t  ~  ■* 

rr  “  V  X  5  «  TVs  ~  VH  +  ~  y.T  +  f  +  f' 
3t  -  e  s 

P 


(1.2.6b) 


with 


pbf '  «  (f  -  f  )pb  +  p(Ts  -  h)$b 
s  s 


(1.2.66) 


[i.n]+[ (Ts  -  h)n]p  +  b(pT)'($s)'  -  p'(^h)')  +  p(Ts  -  h)$b 


and 


;  =  V  x  V 


(1.2.67) 


Observe  that  the  pressure  term  is  eliminated  from  the  new  force  term  f ' , 
and  the  appearance  of  H,  the  total  energy  of  the  average  flow,  instead 
of  H. 
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STREAMSURFACE  APPROACH 

AVERAGING  APPROACH 

CONTINUITY  EQUATION 

(OB)  +  V(p  v3)  -  0 

-£■  Gb)  ♦  V(  pdb)  -  0 

MOMENTUM  CONSERVATION  EQUATION 

(pBv)  +  V  [  (pv  ■  v  ♦  pi  -  t)  B  ]  -  0B^e  ♦  p»? 

— ■  (pb\J)  +  V  KoV  ■  V  ♦  p  -  t^)  b  ]  ■  pbf^  *  “  *(pv"  ■  V”B) 

o?  -  -  |  33  l  (pi  -  t) .n  ] 

0?  ■  -g-  I  (pi  “  t)  .n  1 

ALTERNATIVE  FORMULATION 

-  NON-CONSERVATION  FORM 

♦  (v.V)v  -  -  -  V(p" 
jC  p  e  B 

TT  -  -  —  V(pT  -  T  )  *  i  *  t - —  9(pv"  ■  v"b) 

at  p  a  e  B  b- 

fB  -  ?  -  (pi  -  O.VB  -  -  n3j  (pi  -  ') 

r  -  f  — —  (pi  - "  > .  vb 

B  Pb  a 

CROCCO* 

S  FORM 

—  -  V  X  7  -  TVS  -  VH  ♦  -  V.T  ♦  t  ♦  ?_ 

Jt  s  0  e  B 

?  x  c  ■  W!r  -  vh  ♦  —  7."t  ♦  ^  ♦  ?’ 

3t  p  e  s 

C  ■  V  x  v 

C  ■  Vx  ? 

ENERGY  CONSERVATION  LAW 

(pEB)  ♦  v  (  (pvH  -  kVT  -  T.v)B  ]  -  Q.B  ♦  Q.  B 

It  D 

(pih>)  ♦  Ptptffl'  -  G)b  -  O.b  -  7(pv"H"b)  ♦  Qg-b 

Qb  -  !  33  1  (k^  *  1 

Ob  -  6.n  1  -  ^  (  p  atbl 

5  -  kTT  ♦  T.v 

NON-CONSERVATIVE  FORM 

o(~  ♦  V .  V)  H  -  ^  p  ♦  V(kVT  ♦  T.v)  ♦  Q  ♦  Qg 

P  [  ♦  (v.7)?ll  -  +  ^-(pb)  ♦  f  5(Cb)  ♦  Q  *  Qg  -  i  7(pv”H"b) 

Q  -  n.3,  (kvT  ♦  7-v) 

!  B  J 

NOTATIONS 

The  operators  7  -  ( 3 j ,  3^)  are  defined  by 

The  averaged  quantities  are  defined  by 

b. 

7  -  V  -  n^3 

-  l  !  PA  dC3 

'  b 

0 

where  n  is  the  normal  to  the  surface 

where  is  the  averaged  coordinate  direction 

b  -  (b,  -  b  ) 

1  0 

-¥ 

V  is  the  averaged  velocity 

(A]  -  A,  -  Ao 

Table 


I.  :  Quasi-Three  Dimensional  Formulation. 


Streamsurf ace  and  averaging  approach 


SUMMARY 


The  system  of  equations  obtained  by  following  the  flow  along  a 
streamsurface ,  or  the  system  of  equations  found  by  averaging  the 
conservation  laws  over  a  chosen  direction  lead  to  two  equivalent 
representat ions  of  the  two-dimensional  restriction  of  a 
three-dimensional  flow.  The  former  approach  defines  a  rigorous 
iterative  procedure  between  two  families  of  streamsurfaces  allowing  to 
reconstruct  the  complete  three-dimensional  flow.  On  the  other  hand  the 
averaged  equations  define  a  "mean"  flow  in  the  two  remaining  coordinates 
and  time,  which  is  not  associated  with  any  particular  streamsurf ace . 
The  solution  of  this  two-dimensional  mean  flow  requires  however 
informations  from  the  three-dimensionality  of  the  complete  flow,  as 
expressed  by  the  additional  force  and  energy  source  terms  but  also  by 
the  averaged  products  of  velocity  and  energy  fluctuations  around  the 
averaged  flow  quantities. 

The  quasi-three-dimensional  approximation  will  consist  in  introducing 
some  assumptions  about  the  unknown  contributions  from  the 
three-dimensionality  by  imposing  for  instance,  a  given  shape  for  the 
(unknown)  streamsurface  along  which  the  flow  is  followed  or  by 
neglecting  certain  terms  in  the  equations  related  to  the  additional 
fluxes  introduced  through  the  averaging  procedure.  In  both  cases,  the 
simplest  approximation  consists  in  keeping  only  the  influence  of  the 
varying  streamtube  thickness  B  or  b,  and  neglecting  all  other 
influences . 

Both  approaches  are  applied  in  the  field  of  internal  flows,  such  as 
turbomachinery  and  channel  flows,  while  the  averaging  technique  is 
widely  applied  in  the  field  of  hydraulics,  in  particular  in  the  analysis 
of  tidal  flows  of  rivers,  estuaries  and  seas.  Table  1.1  presents  a 
summary  of  the  flow  equations  as  derived  from  the  two  approaches. 
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Figure  1.1.1.  :  General  configuration  of  a  channel  flow 


Figure  1.1.2.  :  Turbomachinery  blade  passage 

* 

-  Hi  - 


F igure  1 , 1 . 3«  :  S  surface  generated  by  t1  and  £2  lines  and  defined  by 
£3  -  constant. 


-  U2  - 


F igure  1.1.4.  :  Elementary  volume  of  a  streamsheet  of  thickness  B. 
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Figure  1.1.5.  :  Divergence  of  streamsheet  and  relation  to  thickness  variation 


Figure  1.2.1.  :  General  configuration  of  a  channel  passage 
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PART  2  :  THE  COMPUTATION  OF  TURBOMACHINERY  FLOWS 


The  second  part  of  this  report  is  devoted  to  a  presentation  of  some  of 
the  computational  models  and  techniques  used  in  the  prediction  and 
analysis  of  multistage  turbomachinery  flows. 

Turbomachinery  internal  flows  constitute  one  of  the  most  complex  flow 
patterns  in  the  field  of  technology.  The  detailed  flow  behaviour  within 
a  blade  row  is  mostly  strongly  three-dimensional,  viscous  dominated,  and 
often  influenced  by  the  presence  of  shock  waves.  Figure  2.1.1  is  a 
sketch  of  the  typical  flow  phenomena  in  a  compressor  blade  passage, 
where  the  various  contributions  to  the  overall  flow  picture  are 
illustrated.  Due  to  the  turning  of  the  flow  through  the  blade  passages, 
an  initial  boundary  layer  along  one  of  the  end  walls  of  a  blade  row  will 
generate  streamwise  vorticity  within  the  passage,  through  the  secondary 
flow  mechanism,  and  give  rise  to  additional  three  dimensional  effects 
which  can  become  predominant  (as  in  turbine  blades  with  high  turning 
angles),  creating  an  interaction  between  the  viscous  regions  and  the 
three-dimensionality  of  the  flow.  Another  aspect  of  these  phenomena  can 
be  found  in  turbines,  where  the  thick  leading  edges  of  modern  turbine 
blade  sections  generate  horse-shoe  vortices  in  the  inlet  regions  which 
influence  strongly  the  internal  passage  flow.  This  is  illustrated  in 
figure  2.1.2,  showing  a  visualization  of  these  flow  structures.  In 
addition,  the  clearance  space  between  blade  rows  and  casing  (or  hub) 
walls  generates  leakage  flows  and  local  vorticity.  Furthermore,  in 
modern,  high  performance  machines  the  flow  reaches  supersonic  velocities 
with  the  consequence  of  the  appearance  of  shock  waves,  shock  boundary 
layer  interactions  and  choking  phenomena.  Figure  2.1.3  represents  a 
typical  flow  pattern  in  a  two-dimensional  tip  section  of  a  transonic 
compressor  fan  as  determinated  experimentally  and  figure  2.I.H  shows  a 
possible  model  for  these  complex  processes. 

Even  in  steady  rotating  machines,  the  wakes  of  an  upstream  blade  row 
create  an  unsteady  incident  flow  on  the  subsequent  blades  due  to  the 
relative  motion  between  rotors  and  stators.  This  adds  to  the  overall 
complexity  of  the  flow,  making  it  a  formidable  challenge  for  the 
developer  of  computational  prediction  tools. 

The  above  figures  represent  typical  flow  patterns  in  single  blade 
passages.  In  multistage  configurations,  figure  2.1.5,  the  accumulation 
of  these  effects  can  lead  to  severe  deterioration  of  initially  uniform 
velocity  profiles  over  the  height  of  the  blades,  even  when 
circumferential  averaged  profiles  are  considered.  Figure  2.1.6 
illustrates  some  aspects  of  the  flow  migration  as  visualized  by  the 
stagnation  pressure  distribution  in  a  multistage  compressor.  The 
isopressure  lines  show  the  distribution  and  concentration  of  losses  in 
the  exit  plane  of  this  blade  row. 

Notwithstanding  these  complexities,  the  design  process  requires  the 
ability  to  predict,  at  least  in  some  approximate  way,  the  flow  behaviour 
in  a  multistage  machine  in  order  to  optimize  the  performance  and  to  be 
able  to  localize  the  sources  of  dissipation  and  energy  losses.  In  view 
of  the  complexity  of  the  flow,  the  computation  of  a  flow  structure  such 
as  depicted  in  figure  2.1.1  requires  high  levels  of  approximation  of 
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flow  models,  such  as  is  the  full  Navier-Stokes  equations  with  adequate 
turbulence  models  for  the  Reynolds  averaged  formulation.  Methods  able 
to  predict  aspects  of  the  viscous  three-dimensional  nature  of  the  blade 
row  passage  flow  are  appearing,  and  will  be  used  at  increased  rates  in 
the  future.  However,  these  computations,  very  instructive  at  the  single 
passage  level,  will  not  be  able  to  be  used  currently  in  initial  design 
processes  of  a  multistage  machine,  at  least  in  the  foreseeable  future  of 
computer  resources.  Therefore  typical  multistage,  simplified  models 
have  been  developed  in  order  to  single  out  an  approximate  view  and 
prediction  of  the  main  aspects  of  the  turbomachinery  flow. 

For  instance,  one  of  the  aspects  of  the  simplified  models  deals  with  the 
properties  of  the  flow  averaged  over  the  blade  spacing  in  the 
circumf erential  direction,  that  is  in  the  direction  of  the  blade 
rotation.  This  pass age -aver aged  flow  can  be  considered  for  a  rotating 
blade  row,  as  the  flow  seen  by  a  fixed  observer,  attached  to  the  stator 
blade  rows,  looking  at  the  rotor  flow.  Referring  to  figure  2.1.1,  if 
the  flow  is  averaged  at  a  given  spanwise  position  over  the  blade 
spacing,  the  details  of  the  internal  three-dimensional  flow 
configuration  are  lost,  but  instead,  a  two-dimensional  flow  distribution 
in  function  of  mainstream  and  spanwise  coordinates  is  obtained.  Of 
course,  the  lost  information  will  have  to  be  supplied  by  some  form  of 
empirical  input  as  will  be  discussed  in  the  following.  This 
passage-averaged  flow  actually  removes  the  secondary  and 
three-dimensional  flow  details  as  well  as  any  flow  distribution  in  the 
blade-to-blade  direction.  It  represents  therefore  an  axisymmetrio 
approximation  of  the  real  flow  and  is  called  the  meridional 
through-flow,  or  radial  equilibrium  flow.  From  this  point  of  view,  the 
alternative  representation  of  the  through-flow  by  the  flow  projection 
along  one  representati ve  streamsurf ace  of  S2-type  shown  in  figure  1.1.2, 
is  an  equivalent  approximation  and  either  the  stream-surf ace  approach 
discussed  in  section  1.1  or  the  averaged  formulation  of  section  1.2  can 
be  used.  Hence,  the  flow  description  in  rotating  machines  will  rely  on 
the  quasi-three-dimensional  spatial  approximation  described  in  Part  1. 
Indeed,  if  the  through-flow  is  determined,  it  can  be  coupled  to  flow 
computations  in  the  blade-to-blade  surfaces,  of  the  SI -type  of  figure 
1.1.2,  in  order  to  reconstruct  an  approximate  picture  of  the 
three-dimensional  flow.  A  still  further  step  consists  in  determining 
the  through-flow  along  various  S2-type  of  surfaces  attempting  in  this 
way  to  reconstruct  the  secondary  flow  pattern  in  the  blade  passage. 
This  last  step  can  be  considered  for  simple  blade  passages  but  is  seldom 
performed  in  multistage  configurations.  A  more  frequent  approach 
towards  this  goal  consists  in  computing  separately  the  secondary  flow, 
see  for  instance  Atkins  and  Smith  (1983),  based  on  inviscid  secondary 
flow  approximations  as  developed  by  Hawthorne  (1955),  Smith  (1955),  and 
many  others.  The  interested  reader  will  find  an  excellent  and  rigorous 
presentation  of  secondary  flow  theories  in  Horlock  and  Laksminarayana 
(1971).  In  all  these  approaches  the  flow  is  decomposed  in  combinations 
of  two-dimensional  flow  components  which  can  be  reconstructed  in  a 
quasi-three-dimensional  way,  but  which  are  individually  easier  to  solve. 

Consistent  with  these  views  and  with  the  loss  of  information  attached  to 
the  axisymmetric  or  passage-averaged  approximation,  the  details  of  the 
viscous  effects  and  mechanisms  are  smeared  out  and  the  appropriate 
approximation  which  fits  into  this  framework  is  the  distributed  loss 
model .  In  this  model, the  viscous  stresses  are  replaced  by  an  external, 


distributed  friction  force  which  is  responsible  for  the  overall  energy 
losses.  The  lost  information  with  regard  to  the  viscous  shear  stresses 
will  have  to  be  replaced  by  empirical  informations  defining  the  energy 
losses  as  expressed  by  the  enthopy  production  at  the  level  of  the 
through-flow  computations.  In  a  full  quasi-three-dimensional 
computation  this  information  might  also  be  provided  by  two-dimensional 
blade-to-blade  viscous  flow  calculations  including  estimations  of  the 
three-dimensional  effects  on  the  loss  distributions. 

The  flow  models  for  multistage  turbomachines  are  therefore  defined  by 
the  combination  of  the  quasi-three-dimensional  representation  with  the 
distributed  loss  model.  Since  the  distributed  loss  model  is  essentially 
a  non-viscous  model,  coupling  with  boundary  layer  computations  can  be 
added,  in  particular  in  the  through-flow,  where  the  interaction  with  the 
end  wall  boundary  layers  is  an  important  aspect  of  the  overall 
prediction  schemes,  Mellor  and  Wood  (1971),  Hirsch  (197*0,  De  Ruyck, 
Hirsch  and  Kool  (1979),  De  Ruyck  and  Hirsch  (1981),  (1983). 

Section  2.1  will  summarize  the  system  of  equations  used  in  the 
turbomachinery  flow  models  namely  the  distributed  loss  model.  The 
following  sections  will  discuss  and  define  the  quasi-three-dimensional 
formulation  and  the  two  options,  the  streamsurface  approach  or  the 
averaging  technique  will  be  developed  in  sections  2.2  and  2.3 
respectively,  and  the  differences  and  similitudes  pointed  out.  Section 
2.H  will  describe  the  through-flow  computational  techniques.  These 
include  the  widely  used  streamline  Curvature  method  as  well  as  the 
methods  based  on  the  streamfunction  formulation  discretized  by  finite 
differences  or  by  finite  elements. 

The  interested  reader  may  also  find  some  general  reviews  of 
computational  methods,  aimed  at  turbomachinery  flow  applications,  in 
Japikse  (1976)  and  more  recently  in  McNally  and  Sockol  ( 1 98 1  ) . 


2.1.  THE  DISTRIBUTED  LOSS  MODEL 


The  distributed  loss  model  retains,  from  all  the  effects  of  heat 
conduction  and  viscous  shear  stresses,  the  contributions  to  the  entropy 
production.  That  is,  the  dissipation  or  energy  loss  resulting  from  the 
frictional  action  of  the  shear  stresses  is  considered  as  the  dominating 
contribution.  As  a  consequence,  this  action  is  expressed  by  a 

distributed  friction  force  F^,  defined  at 

every  point  in  the  flow  field.  From  a  theoretical  point  of  view,  the 
friction  force  F^,  is  equal  to  the  gradient  of  the  shear  stresses, 


PF, 


V.T 


(2.1.1) 


and  the  contributions  of  heat  conduction  and  of  the  work  of  the  shear 
stresses  to  the  overall  energy  balance  are  either  separately  neglected, 
or  considered  to  cancel  each  other  , 

7(kT7T)  ♦  V(t.w)  *  0  (2.1.2) 


The  contribution  to  the  entropy  production  reduces  to  the  work  of  the 
friction  force  as  given  by 


T 


ds 

dt 


-►  -4 


(2.1.3) 


4 

where  is  considered  to  be  in  the  direction  opposite  to  the  relative 

4 

velocity  vector  w,  and  is  defined  by  this  relation  instead  of  equation 

(2.1.1). 

With 


4 


“F_ .  1 
f  w 


(2.1 .U) 


4  4 

where  1  is  the  unit  vector  in  the  direction  of  the  velocity  w,  equation 
(2.1,3)  becomes 


T  ”  s  wF  (2.1.5) 
dt  f 

Hence,  the  entropy  equation  (2.1.5)  is  added  to  the  system  of  flow 
equations  as  an  independent  equation  and  subsequently,  one  of  the  three 
components  of  the  momentum  equation  is  redundant  and  can  be  left  out  of 
the  system. 


-  H9  - 


With  the  neglection  of  heat  conduction  effects,  of  external  forces  f 

e 

and  of  external  heat  sources  q^,  the  system  of  flow  equations  within  the 

distributed  loss  model  will  be  defined  in  a  relative  reference  system 

rotating  with  a  constant  angular  velocity  u>  around  an  axial  direction. 

■+■  ■+  •+ 

At  a  radius  r,  the  relative  system  has  a  rotation  speed  equal  to  u  *  u  x 
■* 

r.  The  system  of  equations,  defining  the  turbomachinery  flows  within 
this  approximation  level,  is  therefore 

•|£+^(pw)»0  (2.1.6) 
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p  3t 


3s 
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(2.1.7) 

(2.1.8) 

(2.1.9) 


The  rothalpy  I  is  defined  as 


I 


u.  v 


(2.1.10) 


in  function  of  the  static  enthalpy  h  or  the  stagnation  enthalpy 


H  =  h 


(2.1.11) 


and  c  is  the  vorticity  of  the  absolute  flow  v, 


x  v 


(2.1 .12) 


Equation  (2.1.7)  is  in  Crocco's  form  and  is  most  useful  in  the 
computation  of  turbomachinery  flow  because  of  the  direct  use  of  the 
energy  equation  (2.1.8)  and  of  the  entropy  equation  (2.1.9).  However 
the  standard  form  of  the  momentum  equation  can  also  be  applied  and  an 
alternative  form  of  the  above  system  is  provided  by  the  equations  in 
conservation  form. 

♦  V(pw)  =  0  (2.1.13) 
|r  (pw)  +  V(pw(^w)  »  -Vp  -  2 p ( a)  x  w)  -  pw  x  (co  x  r)  +  pF^  (2. 1.1*4) 
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(2.1 .15) 


ft  (pI)  +  ^(p”n  ’  ft 


(ps)  +  ^(pws)  -  ^  pFf  (2.1.16) 


The  fluid  constitutive  equation  has  to  be  added  to  the  above  system  and 
in  the  following,  the  perfect  gas  model  will  be  considered  as 
representing  with  sufficient  accuracy  the  fluid  state.  This  is  a  good 
approximation  for  gases  and  can  be  used  also  for  steam,  at  least  for  dry 
steam.  For  saturated  steam  the  steam  tables  have  to  be  used,  but  one 
can  still  in  this  case  refer  to  the  perfect  gas  model  approximation,  at 
least  locally. 

Hence,  with  the  perfect  gas  law  p/p=rT  the  entropy  variations  are 
connected  to  the  variations  of  stagnation  pressure  p0  and  stagnation 

temperature  T0 


ds 


r 


dp0 
P  0 


(2.1.17) 


Note  that  all  equations  are  written  here  in  the  relative  system,  that  is 
for  a  rotor  flow.  The  stator  flows  are  obtained  by  setting  u=0,  and 
from  the  velocity  composition  law, 


v  =  w  +  u 


(2.1.18) 


w  is  to  be  replaced  by  the  absolute  velocity  v. 

Equation  (2,1.7)  shows  that  the  multistage  turbomachinery  flow  models 
can  be  considered  as  basically  inviscid,  as  they  do  not  contain 

explicitly  shear  stress  gradients  and  F^.  is  considered  as  an  external 

friction  force  defined  by  equations  (2.1.*Q  and  (2.1.17)  when 
information  is  provided  with  regard  to  the  stagnation  pressure  losses. 
In  addition,  the  momentum  equation  (2.1.7)  also  shows  that  the  general 
flow  pattern  will  be  rotational  due  to  the  gradients  of  entropy,  the 
gradients  of  total  energy,  or  rothalpy  I,  and  due  to  the  presence  of  the 

distributed  friction  force  F  . 


Incompressible  flows 

For  hydraulic  turbomachines  such  as  pumps  and  hydraulic  turbines,  the 
incompressible  fluid  approximation  is  best  adapted.  In  this  case,  the 
system  of  turbomachinery  flow  equations  simplifies,  with  p-constant,  to 


V.w 


(2.1.19) 
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(2.1.20) 
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3t 


wxc 


3t  2  p 


(w.$)p* 


(2.1.21) 


where  p*  is  the  rotary  stagnation  pressure  defined  by 

*2  *2 

*  w  u 

p#-p+p— -p— 


(2.1.22) 


Note  that  (2.1.21)  is  obtained  from  (2.1.20)  by  scalar  multiplication 

with  the  relative  velocity  w  and  is  therefore  not  an  independent 
equation.  This  equation  replaces  the  entropy  equation  (2.1.9). 
Actually,  equation  (2.1.21)  is  used  to  compute  the  friction  force  F^. 

through  the  knowledge  of  the  rotary  stagnation  pressure  loss 
coefficients  which  have  to  be  provided  to  the  calculation  procedure. 

Compared  to  the  compressible  system  of  equations,  (-p*)  plays  the  role 
of  the  entropy  since  p*  decreases  in  the  downstream  direction  as  can  be 
seen  from  equation  (2.1.21)  in  steady  state  flows. 


Steady  Flow  Assumption 

It  is  also  customary,  in  the  field  of  turbomachinery ,  to  assume  that  the 
relative  flow  is  steady  for  constant  rotation  speed  of  the  machine. 
This  will  be  the  case  if  the  inlet  flow  is  uniform  in  the  tangential 
direction,  which  implies  that  the  effects  of  the  wake  of  upstream  blade 
rows  are  neglected.  Indeed  one  has  the  following  relation  between  the 
local  time-dependency  in  the  relative  and  absolute  systems 


2a 

3t 


A 


(2.1.23) 


when  6  represents  the  angular  coordinate  in  the  wheel  speed  direction 
and  A  and  R  denote  absolute  and  relative  systems  respectively. 
Therefore,  if  the  absolute  outlet  flow  from  a  preceding  stator  is 
steady,  the  flow  relative  to  the  rotor  will  be  steady  if  the  stator 
outlet  flow  is  uniform  in  the  tangential  direction.  This  will  generally 
not  be  the  case,  due  to  the  presence  of  wakes  and  other 
three-dimensional  effects,  and  hence  the  physical  flow  always  contains 
unsteady  components.  A  recent  detailed  experimental  investigation 
performed  by  Dring  et  al.(1982)  clearly  shows  the  importance  of  the 
unsteady  flow  and  pressure  variations,  although  the  time  average  flow 
was  observed  to  be  in  good  agreement  with  the  steady  state  predictions. 
Another  source  of  unsteadiness  is  due  to  the  inviscid  upstream 
propagation  of  pressure  waves  from  a  moving  rotor  to  an  upstream  located 
stator.  This  effect,  of  subsonic  nature,  can  be  locally  important  on 
the  loading  of  the  stator  trailing  edge  or  on  the  heat  transfer 
coefficients  and  exists  even  for  uniform  tangential  flows.  Therefore, 
the  detailed  interaction  between  two  consecutive  blade  rows  is  always 
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unsteady.  However,  if  one  neglects  this  interaction  or  if  time  averaged 
flows  are  considered,  the  steady  state  assumption  becomes  a  valuable 
approximation. 


In  this  case,  the  distributed  loss  model  simplifies  to 

$(p3)  =  o 

-w  X  c  =  T^s  -  $1  +  F 
(w.$)I  =  0 

.4  LJ 

(w.V)s  =  j 


(2.1.24) 

(2.1.25) 

(2.1.26) 

(2.1.27) 


for  a  compressible  fluid. 

In  particular  the  energy  equation  takes  the  simple  form  of  constancy  of 
the  rothalpy  I  along  a  streamline. 

This  is  nothing  else  than  the  Euler  equation  for  the  energy  exchange  in 
a  rotating  blade  row.  If  subscripts  1  and  2  denote  respectively  the 
inlet  and  outlet  planes  of  a  rotor  then  from  Inconstant  one  has  indeed, 

AH  =  H2  -  Hi  *  (uv^)2  ~  (uv  )|  ■  A(uv^)  (2.1.28) 


■+ 


where  v  is  the  projection  of  the  absolute  velocity  v  on 

4  4  4 

of  the  rotor  velocity  u  *  u>  x  r. 


the 


direction 


As  is  well  known,  the  stagnation  enthalpy  variation  AH  is  the  power  per 
unit  of  mass  flow  exchanged  between  the  fluid  and  the  rotor  blades. 
This  most  widely  applied  equation,  which  allows  the  determination  of  the 
power  and  of  efficiencies,  from  the  knowledge  of  the  flow  velocity 
components,  is  therefore  only  strictly  valid  for  steady  relative  flows. 
The  complete  energy  equation  in  the  relative  system,  contains  also  the 
effects  of  heat  conduction  and  of  viscous  energy  exchange  between 
streamtubes  described  by  the  two  terms  of  equation  (2.1.2). 

Hence,  it  is  essential  to  keep  in  mind  that  the  Euler  equation  for 
turbomachinery  rotors  is  an  approximation,  valid  for  steady,  inviscid 
behavior.  The  viscous  interactions  and  unsteadiness  which  are  always 
present  to  some  degree  in  turbomachinery  flows,  introduce  variations  in 
rothalpy  along  flow  paths.  This  is  of  importance  when  assessing  the 
results  of  computations  based  on  the  simplified  flow  model  described 
here.  It  is  known  from  experimental  observations  and  from  comparisons 
with  fully  viscous  computations,  Moore  and  Moore  (1980),  that  the 
overall  incidence  of  these  effects  can  reach  a  few  percent  on  global 
energy  exchange  estimations,  but  might  still  be  more  important  locally 
in  regions  of  high  viscous  effects,  for  instance  in  end  wall-blade 
corners  and  tip  regions.  For  an  incompressible  flow,  one  obtains 
similarly  the  steady  flow  model  defined  by  the  equations 
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I  Q. 


$.w  =  0 


(2.1.29) 


■»  ■»  1 


(w.$)p*/p 


$p*  +  F  (2.1.30) 

-  -wFf  (2.1.3D 
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2.2.  QUASI-THREE-DIMENSIONAL  APPROXIMATION  FOR  TURBOMACHINERY  FLOWS  : 

THE  5TREAMSURFACE  MODEL 


According  to  the  considerations  developed  in  chapter  3,  two  ways  are 
open  in  order  to  introduce  a  quasi-three  dimensional  approximation  of 
the  turbomachinery  flows. 

Either  one  can  follow  the  flow  along  families  of  streamsurf aces ,  defined 
by  the  positions  occupied  by  the  fluid  particles  lying  initially  on  a 
given  line  in  the  inlet  section  of  flow  passage,  or  one  can  consider  the 
averaged  flow  over  a  selected  direction  and  calculate  this  averaged  flow 
in  function  of  the  remaining  (two  dimensional)  variables. 

Both  methods  will  be  developed  here,  since  they  form  two  equivalent 
formulations  for  through-flow  and  quasi-three-dimensional  computations. 
The  first  approach,  introduced  by  Wu  (1952),  has  been  applied  by  a  large 
number  of  authors,  among  which  KatsSnis  (196*0,  (1966),  Novak  (1967), 
March  (1968),  Frost  (1970),  Adler  and  Krimerman  (197*0,  Biniaris  (1975), 
Novak  and  Hearsey  (1977),  Bosman  and  El-ShaSrawi  (1977).  A  more 
complete  list  of  references,  including  program  development  descriptions 
which  have  been  published  outside  the  journal  literature  can  be  found  in 
the  review  article  by  G.  Serovy  in  the  AGARD  report,  edited  by  Hirsch 
and  Denton  ( 1 98 1 ) .  The  second  approach  is  applied  essentially  by  Smith 
(1966),  Hirsch  and  Warzee  (1976),  (1979),  Jennions  and  Stow  (198*1). 

Two  families  of  surfaces  are  generally  considered  of  the  form 
illustrated  in  figure  1.1.2.  The  SI  family  is  of  the  blade-to-blade 
type  and  can  be  considered  as  generated  by  the  particles  situated 
initially  on  a  line  at  a  constant  radius  r.  The  S2  family  is  generated 
by  particles  located  initially  on  a  radial  line  and  is  designated  as  a 
hub-to-shroud  surfaces.  Due  to  the  axisymmetry  of  the  geometry  the 
cylindrical  coordinate  system  is  most  appropriate  as  reference  system 
with  the  z-axis  aligned  with  the  axis  of  the  machine  and  the  angular 
coordinate  0  defining  the  direction  of  the  blade  rotation.  In  the 
following,  the  blade  rotation  direction  will  be  selected  as  the  positive 
6-direction  and  this  will  also  serve  as  the  positive  direction  for  the 
blade  and  flow  angles. 

The  general  form  of  the  streamsurface  equations  have  been  derived 
earlier  and  are  summarized  in  Table  1.  For  a  streamsheet  thickness  B, 
the  conservation  equations  obtained  by  following  the  flow  along  the 
considered  streamsurface  can  be  written  in  a  two-dimensional  space 
defined  by  this  surface. 

In  accordance  with  the  distributed  loss  model,  heat  conduction  and  shear 

stresses  are  neglected,  but  a  friction  force  F^  is  added  to  the  external 

forces.  In  addition,  we  consider  the  flow  relative  to  a  rotating  frame 
of  reference.  As  a  consequence,  centrifugal  and  Coriolis  forces  have  to 
be  considered  in  the  streamsurface  momentum  equations  and  the  stagnation 
enthalpy  H  is  replaced  by  the  rothalpy  I  in  the  energy  equation. 

This  leads  to  the  following  system  of  equations  valid  on  a 
two-dimensional  streamsurface.  Overbars  on  the  gradient  operator 
indicate  that  the  derivatives  are  defined  by  following  the  surface  of 
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normal  n. 


For  a  streamsurface  S  defined  by  ,£2)  in  a  £  -coordinate  system, 
the  derivatives  along  the  surface,  indicated  by  an  overbar,  are  given 
by,  with  3,  representing  3/3£* 


Vg 


^g  ~  n  3sg 


(2.2.1) 


and  we  refer  to  Part  1  for  a  detailed  analysis.  The  normal  n  is  defined 
by  its  components  n^  with 


n  - 

a 


1) 


a-1 ,2,3 


(2.2.2) 


Written  in  the  relative  system  one  obtains  the  streamsurface  formulation 
of  equations  (2.1.13)  to  (2.1.16),  see  Part  1  for  a  full  derivation  and 
Table  1  for  a  summary. 

Continuity  equation 


(pB)  ♦  V  ( pwB)  -  0 


(2.2.3) 


Momentum  equation 

-► 

|^-  (pBw)  +  V[(pt$pw  +  p)B]  -  [(Ff  ♦  f  ~  2uxw  -  u  x  (uxr)]pB 


or 


3_ 

3t 
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wxc  -  TVs  -  VI  +  F  +  I'B 


(2.2.5) 


Energy  equation 

It  <>bi>  * 5  ■  ^fr 


(2.2.6) 


Entropy  equation 

|jT  ( pBs )  ♦  V  ( pBws )  -  ^  pBFf 


(2.2.7) 
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The  additional  body  force  ?  is  defined  by  equation  (1.1.55) 


3sp 


(2.2.8a) 


if  the  metric  coefficients  do  not  depend  on  see  equations  (1.1.57), 
and 


-► 

Pf 


P?B  + 


(2.2.8b) 


The  streamsheet  thickness  B  is  related  to  the  three-dimensional  flow 
behavior  by  equation  (1.1.35)  or  equation  (1.1.40)  for  b*B/h3>  where  h3 

is  the  metric  coefficient  corresponding  to  £3  in  the  curvilinear, 

a 

orthogonal  system  £  . 


VB 
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(2.2.9) 
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where  the  summation  convention  on  the 


repeated  index  a  is  assumed. 


Note  that 


the  absolute  vorticity  c  is  defined  by 


C 


V  x  v 


(2.2. 1 1 ) 


The  hub-to-shroud  equations  on  the  S2  surface  are  obtained  by  applying 
the  above  equations  for  surfaces  £3  6  0(r,z),  while  the  SI 
blade-to-blade  flow  equations  will  be  obtained  by  setting  for  instance 
C  »  r  »  r(0, z) . 


2.2.1.  Equations  for  Hub-to-Shroud  S2  Surfaces  -  Through-Flow  Equations 


In  the  selected  cylindrical  coordinate  system  0  »  0(r,z),  the  absolute 

■+ 

velocity  vector  v  and  the  relative  velocity  w  are  defined  by  their 
projections  along  the  radial,  axial  and  circumferential  (or  tangential) 
directions , 


v  1  +  v  1  +  V-1 

r  r  z  z  0  b 


u  =  Hi  x  r  »  ur  1 , 


(2.2.12) 
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Hence,  the  relative  velocity  components  are  given  by 


w 


v0  +  air 


(2.2.13) 


The  components 
defined  in  the 
example  E. 1 . 4 , 


of  the  normal  n  to  the  streamsurface  6 
natural  coordinate  system  c  ,  that  is  (1^, 
by 


-  6(r, z) 
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and  the  streamsurface  derivatives  are 
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Geometrical  Definition  of  S2-Surfaces  (Hub-to-Shroud  surfaces) 

It  is  common  practice  to  define  the  surface  S2  by  angles  0'  and  e'.  The 
intersection  of  S2  with  the  (r,6)  or  z  «  constant  plane,  makes  an  angle 
e'  with  the  local  radial  direction  and  the  intersection  of  the  surface 
with  a  cylindrical,  r  =  constant  surface,  is  at  an  angle  0'  with  the 
axial  direction,  figure  2.2.1. 


These  angles  are  related  to  the  components  of  the  normal  vector  n,  since 
they  define  uniquely  the  local  orientation  of  the  streamsurface.  They 
are  defined  by,  referring  to  figure  2.2.1, 


36 

r  r-  ■  tan  e '  -  -rn 
9r  r 


(2.2.16) 


9b 

r  jz  *  tan  0  -™z 


(2.2.17) 


and  the  n  vector  is  obtained,  following  equation  (2.2.14),  as 


4  ]  -►  1  -►  1-4 

n  -  -  -  tan  e'  .  1  -  -  tan  0'  .  1  +  -  10  -  -  n 

r  rr  z  r  0  —  r 


(2.2.18) 
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We  will  assume  in  the  following  that  the  S2  surface  is  time  independent. 
The  condition  for  this  surface  to  be  a  streamsurface  is 


With  equations  (2.2.16)  and  (2.2.17),  the  relation  between  the  velocity 
components  valid  on  S2,  is  obtained 

wn  »  w  tan  8'  +  w  tan  e'  (2.2.20) 

0  z  r 


If  an  axial  flow  angle  8  is  introduced  on  a  cylindrical  surface  of 
constant  radius  r,  figure  2.2.2,  by 

tan  8  =  w  /w  (2.2. 21 ) 

U  Zt 

and  a  meridional  angle  o  in  the  meridional  plane  (r-z)  by 

tan  o  =  w  /w  (2.2.22) 

r  z 


one  obtains  the  relation 


tan  8  =  tan  8'  +  tan  o  .  tan  e' 


(2.2.23) 


Generally,  a  flow  angle  8  is  defined  with  respect  to  the  meridional 
velocity  component  w^,  making  an  angle  o  with  the  axial  direction, 


w  =  /w  2  +  w 
m  r  z 

w0 

tan  8  =  — 
w 

m 

Between  8  and  8  one  has  the  relation 

tan  8  =  tan  8  .  cos 

and  from  equation  (2.2.11) 

tan  8  =  cos  o  .  tan  8'  ♦  sin 


(2.2.24) 

(2.2.25) 

(2.2.26) 
o 

o  .  tan  e'  (2.2.27) 


Since  w  is  in  the  surface  S2,  this  relation  can  also  be  interpreted  as 
defining  the  angle  8  of  the  intersection  of  the  S2  surface  with  a 
conical  surface  forming  an  angle  o  with  the  axial  direction. 

Blockage  coefficients,  friction  and  body  forces 
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The  dimensionless  streamtube  thickness  b,  also  called  the  tangential 
blockage  factor,  is  related  to  the  circumferential  velocity  variations 
of  the  flow  by  equation  (2.2.10),  with  b  -  B/r.  In  order  to  distinguish 
between  the  two  families  of  surfaces,  S2  and  SI,  we  will  indicate  the 
thickness  of  the  S2  streamsheet  by  a  subscript  2.  Hence  we  have  here  b2 

-  B2/r  and  equation  (2.2.10)  becomes 


^  db2 
b2  dt 


-  (-tan  6' 
r 


tan  e' 


3w  3w. 

_ L  +  _£ 

•  30  30 


1  3b2  3b2 

r—  (w  r —  +  w  - — 
b2  r  3r  z  3z 


(2.2.28) 


With  the  knowledge  of  the  normal  vector  to  the  streamsurface  and  the 
vector  relation  (2.2.9),  one  has  here 


-► 


(2.2.29) 


For  an  axisymmetric  flow,  the  tangential  blockage  factor  b2  will  be 

constant.  But  inside  blade  rows,  the  blade-to-blade  flow  creates  a 
tangential  variation  of  the  flow  field  and  hence  variations  of  b2.  As 

shown  by  Novak  and  Hearsey  (1977),  the  streamwise  variations  of  the 
tangential  blockage  factor  can  be  quite  strong  in  particular  in  the 
vicinity  of  the  leading  edge  region  of  the  blades.  In  absence  of  any 
other  information,  the  blockage  factor  b2  may  be  approximated  by  the 

tangential  spacing  variation  between  two  adjacent  blades.  A  better 
approximation  can  be  obtained  from  the  knowledge  of  the  flow  in  the 
0-direction,  that  is  from  the  flow  along  the  SI -blade-to-blade  surfaces. 
In  this  case  b2  can  be  represented  by  the  distance  between  two  adjacent 

streamlines . 


The  body  force  fg,  appearing  in  the  streamsurface  momentum  equation,  and 
defined  by  equation  (2.2.8)  is  given  here  by 


— 

-n  3p  n  1  3p 

p  30  p  r  30 


(2.2.30) 


This  force,  which  results  from  the  confinement  of  the  flow  description 
to  a  two-dimensional  surface,  is  directed  along  the  normal,  and  hence  is 
orthogonal  to  the  local  velocity  vector  following  equation  (2.2.19), 
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(2.2.31  ) 


fB  '  “  '  0 


This  relation  actually  allows  the  determination  of  two  force  components 

f_  and  f_  in  function  of  the  third  one  fnft.  One  has  from  equation 
Br  Bz  B0  M 

(2.2.18) 


fBr/fBe  ■  '  tan  £'  •  nr/ne 

(2.2.32a) 

fBz/fBe  -  '  tan  e'  ■  Vn9 

(2.2.32b) 

f  -i  i£ 

Be  p r  ee 

(2.2.33) 

The  direction  of  the  friction  force  F^.  is  defined 

by  equation  (2.1. 4)  as 

■> 

being  opposed  to  the  relative  velocity  vector  w. 

Hence ,  two  components 

of  F  are  determined  by  the  flow  angles  in  function  of  the  third  one,  or 

alternatively  in  function  of  the  magnitude  of  F^, 

we 

F  =  -F  — — 

(2.2.  34a ) 

fe  f  *  w 

w 

F,  -  -F  — 

(2.2.34b) 

fr  f  w 

R 

1 

(2.2.34c) 

f  z  f  w 

Through-Flow  Equations  in  Cylindrical  and  Axisymmetric  Coordinates 


Finally,  the  through-flow  equations  take  the  following  form  within  the 
streamsurface  approach  in  cylindrical  coordinates  when  Crocco's  form, 
equation  (2.2.5)  is  selected  for  the  momentum  equation,  see  Problem 
2.2.1. 


r—  W  -  W 

3t  r  z 


rr-  w  +  w 
3t  z  r 


3  .  .  , 

at  (pb’r)  + 

3_ 

3r 

(pb2rw  ) 
r 

‘fe 

(pb2rw  ) 
z 

-  0 

(2.2.35) 

,a  a 

3r  Wz  3z 

V 

a 

(rve) 

.  t  i? . 

31  +  F  +  f 

(2.2.36) 
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3r 

T  3r 

3r  fr  Br 

.3  a 

3r  Wz  3z 

V 

. 

3_ 

(rve) 

,Tk. 

31  +  f  +  f 

(2.2.37) 
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(2.2.38) 
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3  T  a  31 

TTT  I  +  W  r— 
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1  i-  n 
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The  operator  w.V  is  defined  here  by 


•*  5  a  3  3 

W.V  =  W  t—  «  W  r—  +  W  t — 

m  3m  r  3r  z  3z 


(2.2.39) 

(2.2. 40) 


(2.2. 41 ) 


where  m  indicates  the  distance  along  the  meridional  streamline,  see 
figure  2.2.3. 


Another  interesting  coordinate  system  is  provided  by  the  axisymmetric 
(0fm,n)  system  .  With  the  assumption  that  the  coordinate  n  is  a 
meridional  streamline,  that  is  obtained  by  the  projection  of  the 
velocity  on  a  meridional  plane,  figure  2.2.3.,  then  w^  “  vn  “  ® 

everywhere.  This  leads  to  the  following  system  of  streamsurface 
equations 


3  w6  3  ,  . 

3t  Wm  r  3m  rV6 


fm 


♦  f 
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(2.2.42) 
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♦  f 


Be 


(2.2.43) 


(2.2.44) 


4  4 

The  projection  of  Ff  along  the  direction  n  is  zero  since  Ff  is  directed 
opposite  to  wt  equation  (2.1.4).  The  radius  of  curvature  1 /R  is 
connected  to  the  angle  o  by 


9o 

9m 


-1/R 


m 


(2.2.45) 


The  continuity  equation  becomes. 


9  9 

( pb2rh2 )  +  ^  (prwmb2h2)  -  0 


(2.2.46) 


or,  for  a  steady  flow 
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prwmb2h2  -  constant  along  the  meridional  streamline  m  (2.2.47) 

Equation  (2.2.47)  expresses  the  conservation  of  mass  in  a  streamtube 
around  the  meridional  streamline  m,  see  figure  2.3.1.  Indeed,  h2  (h2 

being  the  metric  coefficient  associated  to  the  n-coordinate)  is 
proportional  to  the  thickness  of  the  streamtube  centered  on  m.  Calling 
this  meridional  streamtube  thickness  B,,  measured  perpendicularly  to  the 

m-line,  one  has  from  equation  (2.2.47) 


pw  b2rB,  -  constant 
m 


(2.2.48) 


Note  that  B2«=b2r  is  the  thickness  of  the  S2  streamsheet,  while  Bi  can  be 
considered  as  the  thickness  of  the  SI  streamsheet. 

As  mentioned  earlier,  it  is  fully  consistent  with  the  already  assumed 
approximations  to  consider  that,  for  steady  inlet  flow  conditions  to  the 
machine  and  steady  rotations,  the  relative  flow  remains  steady.  In  this 
case,  the  set  of  streamsurface-S2-equations  simplify  to  the  following 
form,  in  cylindrical  coordinates 


h  (ppbjWr)  +  h  (prb^w2)  -  0 


,3  9  ,  _  9s  91  W0  9  ,  , 

z'  9r  z  9z  r  9r  9r  r  9z  v  0;  fr  Br 

.9  9  ,  _  9s  91  W0  9  ,  . 

rv9r  z  9z  r ‘  9z  9z  r  9z  u  0;  fz  Bz 


-  w  7—  (rv.)  *  F,.  +  f-. 
r  m  9m  0  f0  B0 


9s  w 
w  —  =s  -  F 
m  9m  T  f 


w  r—  I  -  0 
m  9m 


(2.2.49) 

(2.2.50) 

(2.2.51 ) 

(2.2.52) 

(2.2.53) 

(2.2.54) 


Note  that  one  of  the  above  equations  is  not  independent  from  the  others. 

The  last  equation  (2.2.54)  states  that  the  total  energy  or  rothalpy  I  is 
constant  along  a  meridional  projection  of  the  streamline  along  the 
S2-surface.  This  implies,  the  following  relations  : 

I  -  H  -  uv  -  constant  along  m,  for  a  rotating  blade  row,  or 

(2.2.55) 

H  »  constant  along  the  meridional  streamline  in  a  non-rotating  system 
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Equation  (2.2.52)  defines  the  meridional  variations  of  the  angular 
momentum  (rvQ)  as  determined  by  the  component  of  the  friction  force  and 

the  body  force.  In  practical  computations,  equation  (2.2.53)  determines 
the  magnitude  of  the  friction  force  Ff,  while  f00  will  be  deduced  from 

equation  (2.2.52)  when  the  necessary  information  is  provided  for  the 
variations  of  rv0  as  will  be  discussed  in  the  next  section.  The  entropy 

variations  are  obtained  from  the  total  pressure  loss  coefficients  as 
discussed  next. 

The  remaining  equations  (2.2.50)  and  (2.2.51)  form  a  system  for  the 
velocity  components  w  ,w  .  Since  the  entropy  equation  has  been  used 

P  z , 

explicitly,  the  two  momentum  equations  are  not  independent  and  one  of 
them  can  be  applied ,' while  the  other  is  discarded. 

The  radial  component  of  the  momentum  equation  (2.2.50)  is  known  in  the 
turbomachinery  literature  as  the  radial  equilibrium  equation. 


2.2.2.  Principles  of  Through-Flow  Computations 


Three  sets  of  external  information  have  to  be  provided  to  the  system  of 
equations  (2. 2. ^9)  to  (2.2.5^),  expressing  the  missing  information  with 

regard  to  viscous  effects  (F^)  as  well  as  to  the  properties  of  the  flow 

in  the  omitted  direction,  that  is  the  9  or  tangential  direction 
(t>2  »rv0)  • 

i )  The  tangential  blockage  or  streamsheet  thickness  b2 


This  information  has  to  be  defined  by  the  flow  along  the  other 
family  of  surfaces,  namely  the  SI  or  blade-to-blade  surfaces 
according  to  equation  (2.2.28).  In  absence  of  this  information, 
the  following  approximation  is  usually  applied 

b  _  geometrical  width  of  channel 
2  inlet  width  of  channel 

whereby  b2  is  made  flow  independent.  For  the  S2-surface  this 
implies 

b2  -  1  -  ;  (2.2.56) 

s 

where  d  is  the  blade  thickness  in  the  tangential  direction  and  s 
the  blade  spacing  or  pitch,  see  figure  2.3.1.  In  duct  regions, 
outside  blade  rows,  b2  can  be  taken  equal  to  1. 

ii)  The  angular  momentum  rv. 

. O 

Since  equations  (2.2.50),  or  (2.2.51),  solve  for  the  velocity 
components  in  the  meridional  plane  w^  and  w^,  the  tangential 


row.  From  the  knowledge  of  the  loss  coefficient  w,  one  obtains  the 
entropy  variation 


s  -  s. 


-1  n  ( 1 


(is) 


7> 


(2.2.66) 


from  which  the  friction  force  can  be  defined. 


The  determination  of  the  stagnation  pressure  losses  are  strictly 
outside  the  present  level  of  approximation  and  have  to  be  provided 
by  external  sources  such  as  empirical  informations.  Note  however, 
that  within  the  blade-to-blade  flow  computations,  viscous 
calculations  could  be  performed,  for  instance  with  boundary  layer 
methods  in  order  to  estimate  the  losses  occuring  along  the  blade 
sections.  These  two-dimensional  profile  losses  form  only  a  part  of 
the  total  loss  and  they  would  have  to  be  supplemented  by  additional 
contributions  from  full  three-dimensional  effects  such  as,  for 
instance,  the  secondary  and  tip  leakage  losses. 

Therefore,  as  with  the  flow  angles,  reference  is  made  to  empirical 
correlations  relating  loss  coefficients  to  a  limited  number  of 
geometrical  and  flow  parameters  considered  to  have  a  dominating 
influence.  These  loss  correlations  complete  the  missing 
information  of  the  S2-surface  flow  as  defined  by  the  equations 
(2. 2. ^9)  to  (2. 2. 5*0  above.  A  large  body  of  literature  is 
available  in  this  field  and  a  recent  overview  of  the  problem  and 
limitations  of  loss  and  turning  correlations  can  be  found  in  Hirsch 
and  Denton  ( 1 981 ) . 


iv)  The  calculation  of  the  density 


The  density  is  computed  from  the  knowledge  of  the  flow  velocities 
and  the  local  stagnation  conditions,  for  instance  from  equation 
(2. 1.29),  with 


H' 

o 


C  T' 


P  o 


2  2 

w2  v1/(Y-1)  ,  jnJ_Vl/(Y-1) 

2  H'  ;  u  2  H'  ; 

o  o 


With  the  relation 


(2.2.67) 


(2.2.68) 


1 /(Y-1 ) 


-(s-sA)/r 

e 


(2.2.69) 


one  can  relate  the  local  density  to  stagnation  values  at  another 
point  A  on  the  same  streamline. 
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(2.2.70) 
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2.2.3.  Alternative  forms  of  through-flow  equations 


As  mentioned  earlier,  one  of  the  momentum  equations  is  redundant  and  the 
choice  between  equations  (2.2.50)  and  (2.2.51)  depends  on  the  general 
flow  conf iguraiton.  Equation  (2.2.50),  the  redial  component  of  the 
momentum  equation,  has  a  coefficient  wz  in  the  left  hand  side  and  will 

not  be  suitable  for  a  machine  with  a  predominant  radial  flow,  that  is 
for  radial  or  centrifugal  machines.  Similarly,  the  axial  component  of 
the  momentum  equation  ,  equation  (2.2.51),  cannot  be  applied  for  axial 
machines  where  w^  might  vanish  in  all  or  part  of  the  flow  domain.  For 

mixed  flow  machines  one  would  have  to  use  one  or  the  other  equation 
according  to  the  flow  region.  An  elegant  way  to  avoid  this  has  been 
introduced  by  Bosman  and  Marsh  (197*0  and  consists  in  using  the 

projection  of  the  momentum  equation  on  a  direction  N  situated  in  the  S2 

surface  and  perpendicular  to  the  velocity  vector  w,  see  figure  2.2.1. 

The  N-direction  is  actually  perpendicular  to  the  body  force  f  ,  which  is 
normal  to  S2,  as  well  as  to  the  friction  force  F^,  which  is  parallel  to 
w.  Hence,  one  has 


N  .  F„  -  0  and  N  .  f  «  0  (2.2.71) 

I  b 


4 

Defining  N  by  the  vector  product 

N  -  (fg  x  w) 

the  projection  of  the  streamsurface  momentum  equations  can  be 
from  equation  (2.2.5)  for  steady  state  flows,  as 

-4  -4  -4 

-N.(w  x  i)  «  N. (TVs  -  VI) 


(2.2.72) 


obtained 


(2.2.73) 


or 
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(2. 2. 7*0 


(?b.U.5‘  -  T(5.V)s  -  (5.5)1  -  T(H  f  *  N  If)  -  (Nr  |i  .  Hz  fi) 

With  equations  (2.2.32)  and  (2.2.33)  defining  the  direction  of  the  body 
force  fg,  equation  (2.2.74)  becomes 
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(2.2.75) 


,  nr  ,  ,  3s  31.  w2rnr  3  ,  ,  nz  3  , 

-  (w  -  —  w  )  (T  r  ■  5—)  ♦  ~  [—  x-  (rv  )  -  —  5-  (rv. )] 
r  n  0  dz  3z  r  n.  3z  0  nQ  3r  0 

b  b  b 


The  coefficient  in  the  left  hand  side  of  the  above  equation  is  now  w2, 
independent  of  the  flow  direction.  This  equation  has  another  advantage 
in  that  the  force  terms  have  disappeared  from  the  right  hand  side, 
although  they  are  not  neglected.  Hence,  equation  (2.2.75)  does  not 
require  the  explicit  calculation  of  the  force  components,  but  does 
however  require  the  explicit  knowledge  of  the  normal  to  the 
streamsurface ,  that  is  the  angles  c'  and  S'.  Introducing  these  angles 
explicitly,  equation  (2.2.75)  becomes 


->2,3  9  .  ,  _  o.w^  9s  31, 

w  %  wr  "  3F  wz)  “  (w~  +  w“*tan  B  )(T  ~ 
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(2.2.76) 


..L, 
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r  0  3z  3z  r  3r  0  3z  0 


In  a  duct  region,  where  no  blades  are  present  there  is  no  body  force 

because  of  the  axisymmetric  assumption  and  the  direction  N  cannot  be 
defined.  In  this  case,  because  of  the  axisymmetry,  one  can  select  the 
direction  perpendicular  to  the  meridional  streamline  m  in  the  (r-z) 
space  as  the  projection  direction,  see  figure  2.2.2.  With 


x  w 


1  - 
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(2.2.77) 


one  obtains  instead  of  equation  (2. 2. 7*0, 
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(2.2.78) 
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Note  that,  in  a  duct  region  one  can  consider  the  flow  from  the  absolute 

reference  system  and  set  I=H,  w-v  since  u-0.  Also  the  overbars  on  the 
derivatives  can  be  suppressed  because  the  flow  is  considered  as 
axisymmetric. 


2.2.4.  Properties  of  the  through-flow  equations 


It  is  important  to  analyze  the  mathematical  properties  of  the  system 
formed  by  the  continuity  equation  and  the  momentum  equation,  for 
instance  the  radial  component  for  an  axial  machine,  with  regard  to  the 
characteristics  and  their  elliptic  or  hyperbolic  properties.  This  is 
best  obtained  after  setting  the  equations  into  the  quasi-linear  form. 
The  radial  momentum  equation  is  best  written,  for  this  analysis,  in  the 
form  of  equation  (2.2.4),  after  introduction  of  the  continuity  equation 
(2.2.3),  for  steady  flows 
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F  ♦  f  +  v2/r 
fr  Br  0 


where  the  right  hand  side  is  unimportant  for  the  properties  of  the 
system.  The  continuity  equation  is  transformed  to  the  quasi-linear 
form,  developing  equation  (2.2.49) 
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(2.2.80) 


The  derivatives  of  density  and  pressure  can  be  expressed  in  function  of 
the  velocity  derivatives.  One  has,  from  the  definition,  equation 
(2.2.68),  for  any  variation  6 

<5p  -  -  fr  ( w  6w  +  w  6w  +  w  6w  )  (2.2.81) 

a  r  i  Z  Z  0  o 

and 

6p  -  (|^)  6p  -  a2  6p  (2.2.82) 

dp 
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where  a 2  is  the  square  of  the  speed  of  sound. 

The  dependent  variables  are  w^  and  w and  as  discussed  above,  w^  is 
supposed  to  be  known  from  the  blade-to-blade  information. 

If  w„  is  defined  by  the  flow  angle  &,  equation  (2.2.21),  then  wQ  will  be 
function  of  and  contribute  to  the  definition  of  the  properties  of  the 
equation.  In  the  other  case,  if  wQ  is  given  as  a  function  of 
coordinates,  as  in  the  design  option,  wQ  will  be  considered  independent 
of  the  velocity  component  w^.  We  will  combine  these  two  cases  by 
writing  formally 


6w  -  v.tan  B.6w 

V  Z 


(2.2.83) 


where  v»1  when  the  flow  angle  6  is  given  and  v»0  when  w  is  imposed. 

b 

Hence,  equation  (2.2.81)  can  be  written  as 

P 


SP 


„  [w  6w  +  w  (1  +  v  tan2B)6w  ] 
a  r  r  z  z 


(2.2.811) 


The  continuity  and  radial  momentum  equations  can  be  written  under  the 
system  form  as  follows,  writing  K  -  1  +  v  tan28 


(1-M2)  f-w  -  M  M  K  £-  w  -MM  |-  w  ♦  (1-KM2)  |-  w 
r  3r  r  r  z  3r  z  z  r  3z  r  z  3z  z 


w.  VB  2 
B  2 


(2.2.85) 


9  3  0 

■wKt—  w  +  w  — —  w  -  F  +  f  +  — 
z  3r  z  z  3z  r  fr  Br  r 


(2.2.86) 


M  and  M  are  the  Mach  numbers  corresponding  to  w  and  w  ,  that  is 
r*  z  r  z 


M  -  w  /a 
r  r 


and 


M  «  w  /a 
z  z 


(2.2.87) 


In  matrix  form,  the  above  system  can  be  written  as 
3 


A  . 


3r 


w 

- 

w 

q, 

r 

+  B  •  h 

r 

w 

w 

Qj 

z 

z 

(2.2.88) 


or  explicitly 
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X 

t. 

1 

w  U 

X 

1 

3 

wr 

-M  M 
r  z 

1  -  KM2 
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3_ 
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r 

Qi 

0  -  w  K 

z 

3r 
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z 

w 

z 

0 

3z 

w 

z 

q2 
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where  the  right  hand  sides  represent  the  terms  independent  of  the 
velocity  derivatives.  The  system  (2.2.89)  is  elliptic  or  hyperbolic  if 
the  characteristics  are  imaginary  or  real.  The  characteristic 
directions  A  are  determined  by  the  solutions  of  the  determinant 
|A+AB|«0.  That  is 


1 


M2  -  AM  M 
r  r  z 


Aw 

z 


-M  M  K  ♦ 
r  z 


A(1  - 


KM2) 

z 


(2.2.90) 


leading  to  the  two  solutions 


A  »  [M  M  K  ±  /K(M2  ♦  KM2  -  1)]/(1  -  KM2) 
±  r  z  r  z  z 


(2.2.91) 


One  obtains  the  following  important  properties,  valid  for  all 
quasi-three  dimensional  descriptions  along  the  S2  type  of  surfaces  : 

-  The  system  of  equations  (2.2.119)  to  (2. 2. 5*0  is  elliptic  if  the  expression 
under  the  square  root  is  negative,  that  is  for 

M2  +  vM2  <  1  (2.2.92) 

m  0 


Hence,  the  through-flow  model  is  elliptic  for  subsonic  meridional  Mach 

numbers  when  w„  is  imposed  (v=0),  while  the  system  is  elliptic  for 

....  —  u  A 

subsonic  relative  Mach  numbers  M2«M2  +  M2  <1,  when  the  flow  angle  6  is 

given  ( v* 1 ) •  Any  solution  method  will  have  to  take  these  properties 
into  account. 

The  first  case,  namely  the  design  option,  whereby  w^  is  considered  as  a 

known  function,  is  of  particular  interest  and  deserves  some  comments. 
The  condition  (2.2.92)  indicates  that  the  through-flow  problem  can  be 
treated  in  an  elliptic  way  for  subsonic  meridional  Mach  numbers, 
although  the  relative  Mach  number  can  be  supersonic.  Hence,  for  given 
w  distributions,  one  can  handle  relative  supersonic  flows  inside  the 

blade  rows  in  the  same  way  as  the  subsonic  relative  flows  and  compute 
the  meridional  velocity  components  within  the  blade  rows  by  implementing 
in  the  appropriate  way,  the  information  about  w0  to  be  transferred  from 

the  blade-to-blade  surfaces.  However,  if  this  information  is 

transmitted  under  the  form  of  the  flow  angles  fS  or  3,  no  calculation 
method  of  "elliptic"  type  (for  instance  the  methods  described  in  section 
2.4.)  will  be  able  to  handle  correctly  supersonic  relative  flows,  at 
calculation  points  situated  inside  the  blade  rows,  unless  the 
computational  procedure  is  adapted  for  supersonic  flows,  by  allowing  a 
correct  treatment  of  hyperbolic  regions.  For  instance,  by  applying 
techniques  similar  to  those  developed  for  the  computation  of  transonic 
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potential  flows,  or  time  dependent  formulations  of  the  type  applied  to 
the  resolution  of  the  system  of  Euler  equations. 

Note  that  outside  blade  rows,  that  is  in  the  duct  regions,  the  above 
conditions  automatically  apply  since  vQ  is  always  given,  as  a 

consequence  of  equation  (2.2.58)  expressing  the  constancy  of  angular 
momentum.  Hence  in  duct  regions  the  problem  will  always  be  elliptic  for 
subsonic  meridional  Mach  numbers. 

This  situation  is  closely  connected  to  the  analysis  of  the  time-like 
direction  in  three  dimensional  supersonic  flows.  An  arbitrary  direction 
is  to  be  considered  as  time-like  if  the  projection  of  the  velocity  along 
this  direction  is  supersonic.  When  this  projection  is  subsonic,  the 
problem  is  elliptic  with  respect  to  the  corresponding  direction,  since 
an  upstream  propagation  of  information  in  this  direction  is  possible  . 
The  same  properties  appear  here,  since  upstream  propagation  of 
information  in  the  meridional  direction  is  possible,  from  outlet  to 
inlet  of  a  blade  row  for  instance,  when  the  meridional  Mach  number  is 
subsonic  even  for  relative  supersonic  flows.  This  apparently 
contradictory  situation  can  be  understood,  even  in  the  presence  of 
choked  blade-to-blade  passages,  because  of  the  way  the  information  is 
transmitted  from  the  blade-to-blade  flow  to  the  S2-surface.  Indeed,  the 
supersonic  flow  occurs  in  the  blade-to-blade-SI -surface ,  which  is  not 
resolved  by  equations  (2.2.49)  to  (2. 2. 5*0  and,  as  discussed  above,  the 
information  on  the  properties  of  this  flow  have  to  be  transferred 
through  the  entropy  (or  stagnation  pressure)  and  w.  variations.  If  the 

flow  passage  is  choked,  implying  that  no  upstream  influence  can  be 
transmitted  in  the  relative  flow  direction,  and  with  the  mass  flow 
defined  by  the  inlet  conditions  independently  of  the  outlet  flow 
conditions,  there  is  still  an  upstream  propagation  of  pressure  waves  in 
the  meridional  direction,  assuming  subsonic  meridional  flow.  This  is 
illustrated  on  figures  2.2.4  and  2.2.5  in  the  case  of  choked  compressor 
and  turbine  passages,  respectively.  In  these  cases,  the  entropy 
variation  s=s(r,z)  or  s=s(r,m)  will  have  to  take  into  account  the 
entropy  discontinuity  over  the  normal  shock  in  figure  2.2.4  and  the 
transmitted  tangential  flow  components  will  have  to  be  consistent  with 
the  choked  mass  flow.  Similarly,  the  value  of  w  transmitted  at  blade 

outlet  of  the  choked  turbine  of  figure  2.2.5,  will  have  to  be  defined  by 
the  local  value  of  the  mass  flow  in  such  a  way  as  to  respect  mass 
conservation  between  the  throat  area  and  the  outlet.  This  has  of  course 
implications  on  the  way  the  computation  of  the  through-flow  is  performed 
since,  for  a  fully  choked  passage,  the  mass  flow  cannot  be  imposed  but 
will  result  from  the  calculations  for  an  imposed  pressure  ratio. 

A  similar  situation  occurs  for  the  supersonic  compressor  inlet  flow 
where  the  phenomena  of  unique  incidence,  at  the  leading  edge  of  the 
blades,  for  supersonic  incident  velocities,  has  to  be  taken  into  account 
and  the  corresponding  tangential  velocity  at  inlet  has  to  be  transferred 
to  the  through-flow  computation.  Similarly,  if  the  flow  is  choked  the 
outlet  tangential  velocity  should  be  consistent  with  the  critical  mass 
flow. 

Therefore,  when  information  on  wQ  is  transmitted,  this  decouples  the 

meridional  velocity  field  (w  ,w  )  from  the  tangential  velocity  field 

r  l* 
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and,  provided  the  wQ  distribution  is  consistent  with  the  physical 

u 

phenomena  in  the  blade-to-blade  surfaces,  one  is  left  with  an  elliptic 
or  hyperbolic  problem  in  the  (r,z)  space  according  to  the  meridional 
Mach  number.  However,  when  the  flow  angle  is  transmitted  this  couples 
the  tangential  flow  field  to  the  meridional  one  through  equation 
(2.2.21)  or  (2.2.25)  and  the  full  three-dimensional  conditions  for 
ellipticity  are  maintained,  namely  subsonic  relative  velocities. 


2.2.5.  Flow  Equations  along  Blade-to-Blade  SI  Surfaces 


The  equations  for  the  blade-to-blade  surfaces  of  the  SI  type  in  figure 
1.1.2,  are  obtained  by  following  the  same  procedure  as  in  the  previous 
section  starting  from  equations  (2.2.1)  to  (2.2.11). 

If  the  SI  surfaces  are  defined  by  r=r(0,z)  families,  the  continuity 
equation  becomes,  indicating  by  bi  the  thickness  of  the  SI  streamsheet, 

see  also  equation  (E. 1.1. 16), 


ft  (B‘p)  +  r  fe  (pB‘V  +  r  h  (pB‘wzr)  *  0 


(2.2.93) 


The  streamtube  thickness  B,  is  defined  by 


1  dB|  1  ,3B>  We  3B>  36 ‘  ,9_  9_  9,  .  . (2.2.94a) 

B,  dt  —  B,  9t  r  90  Wz  9z  9rWr  r‘z  9rWz  +  n0  9r  Vr 


or 


VB, 

~b7 


9n 

9r 


(2.2.94b) 


with 


-  (— ) 

dZ  ^ 

~  (— ) 

^90; 


n  -  1 


(2.2.95) 


as  the  components  of  the  normal  vector  to  the  surface  r=r(0,z).  This 
form  of  the  continuity  equation  is  not  suitable  for  radial  blade  rows, 
where  the  blade-to-blade  sections  are  essentially  in  the  (0-r)  planes 
(when  the  angle  o  comes  close  to  90  degrees)  and  the  SI  surfaces  would 
have  to  be  defined  by  z-z(0,r). 
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Therefore  a  more  general  form  valid  for  any  blade  geometry  is  obtained 
by  defining  surfaces  in  the  (m,6)  coordinate  system  where  ra  is  the 
meridional  projection  of  the  streamlines,  see  figure  2.2.3. 

Hence,  in  the  axisymmetric  orthogonal  coordinate  system  (0,m,n)  shown  on 
figure  2.2.3,  the  SI  surfaces  will  be  defined  by  surfaces  n-n(0,m). 

The  continuity  equation  takes  on  the  general  form  given  by  equation 
(2.2.3) ,  leading  to 


ft  (pB‘>  *  r  fi  <prB'"n>>  *  ?  f?  <pB,V  ‘  0 


(2.2.96) 


where  Bi«bjhj,  h3  being  the  metric  coefficient  associated  to  the 
coordinate  n,  and  with  the  definition  of  the  bi~coeff icient ,  following 
equation  (2.2.10), 


1  9b!  wQ  9b! 

-  ( -  +  -  -  +  u 

bi  9t  r  90 


B  b !  «  .  ri 

- )  =  (_  w  +  n  —  (w  /r)  +  n  —  w  ) 

m  9m  ’  '  9n  n  0  9n  v  Q  m  9n  rn 


(2.2.97) 


with 


9n 

9o 

9n 


(0,m) 


(2.2.98) 


nm  =  "  9S  (°’m) 


The  streamsurface  derivatives  are  defined  by  equation  (2.2,1) 
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90 

90 

n0  9n 

9_  _ 

9_  _ 

9_ 

9m 

9m 

nm  9n 

(2.2.99) 


Some  confusion  might  be  possible  here  between  the  normal  n  and  the  n 
coordinate  of  the  axisymmetric  coordinate  system  nfmf0.  However,  it 
should  be  clear  from  the  text  that  equation  (2.2.98)  for  instance, 

relates  the  0  component  of  the  vector  n  to  the  0  derivative  of  the 

surface  defined  by  relation  n=n(0,m)  in  the  (n,m,0)  coordinate  system. 

For  a  permanent  relative  flow,  the  steady  state  continuity  equation 
becomes 


S  <pB'rV  *  r  h  !pB‘r“e)  ■  0 


(2.2.100) 
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Note  that  B,  is  measured  here  in  a  direction  perpendicular  to  m. 


In  practical  calculations  the  function  Bi=Bi(6,m)  has  to  be  determined 

from  the  hub-shroud  S2-flow  according  to  equations  (2.2.96).  An 
interesting  and  widely  used  approximation  is  obtained  from  equation 
(2.2.48)  which  states  the  mass  conservation  along  a  streamsheet  centered 
on  m. 

Hence  one  can  define  the  function  Bi=Bi(m)  as  follows 


B,(m) 

B  i  ( m  o  ) 


(prw  b2 )  o 
_ m 

prw  b2 
m 


(2.2.101 ) 


where  (pw  )  is  obtained  from  the  hub-shroud  S2  streamsurface  solutions 
m 

to  equations  (2.2.49)  to  (2.2.54)  and  where  the  subscript  0  indicates  a 
reference  state,  generally  the  inlet  of  the  considered  blade  row.  The 
inverse  of  equation  (2.2.101),  namely  the  quantity  B^(m0)/B^(m)  is  often 

called  the  Axial  Velocity-Density  ratio  or  AVDR,  in  the  literature. 


The  momentum  equation  (2.2.5)  becomes  in  the  (0,m,n)  coordinate  system, 

with  w  =0  and  F„  =0  for  a  streamsurface  n=n(m,0) 
n  in 


9  Wm  ,3  3  . 

at  0  r  k3m  0  30  m' 

a  W0  , a  a  , 

r-  w  +  —  (;rr  w  _T“rvn) 
3t  m  r  80  m  3m  0 


r  30  3  ”  r  30  1  +  Ff0  +  fB0 
T3m  3  "  3m  1  +  Ffm  +  fBm 


-w 


Bn 


while  the  energy  and  entropy  equations  can  be  written  as 


ai  ai  w0  ai  i  a  ,  D  . 

at  wm  3m  r  80  pB,  8t  (pBl) 


(2.2.102) 

(2.2.103) 


(2.2.104) 


(2.2.105) 


3s  3s  W0  8s  w  „  (2.2.106) 

-  +  u  -  +  -  -  a  -  K 

at  m  8m  r  80  T  f 


where  R  is  the  radius  of  curvature  of  the  meridional  streamline.  The 
body  force  fg  is  defined  here  by 


pfB 


-► 

n 


3p 


*  3n 


(2.2.107) 
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In  particular,  if  the  SI -streamsurface  is  assumed  to  be  close  to  an 
axisymmetric  streamsurface,  then  SI  is  a  surface  n-constant  and  the  body 

-► 

force  fg  is  purely  directed  along  the  n-axis,  which  is  the  normal  to  the 

m-0  axisymmetric  surface.  This  follows  from  the  same  considerations  as 
those  which  lead  to  equation  (2.2.30).  Under  these  conditions,  one 
would  have 

rue  ■  fBm  ■  0  (2.2.108) 

Equation  (2.2.104)  defines  the  n-component  of  the  body  force  as  balanced 
by  the  centrifugal  force  in  this  direction  and  can  be  considered  to  be 
responsible  for  the  curvature  of  the  meridional  streamline.  The  radius 
of  curvature  Rm  will  be  positive  if  the  pressure  increases  with  n,  that 

is  from  hub  to  shroud. 

In  analogy  with  the  discussion  of  the  S2-streamsurface  computation,  the 
blade-to-blade  SI  computations  require  the  knowledge  of  the  streamsheet 

*4  4 

thickness  Bj,  and  of  the  forces  F.  and  fg. 

The  thickness  B, ,  defined  by  equation  (2.2.97)  can  be  approximated  by 

4 

equation  (2.2.101),  while  the  force  fQ  requires  the  knowledge  of  the 

D 

pressure  gradients  in  the  n  direction  at  every  point.  Both  quantities 
are  to  be  obtained  from  the  S2-streamsurface  solutions.  As  discussed 
previously,  the  friction  force  on  the  other  hand  has  to  be  provided  by 
empirical  data  or  can  be  calculated  at  least  partly  in  the 
blade-to-blade  surfaces  (for  instance  by  a  boundary  layer  approximation) 
and  averaged  out  over  the  blade-to-blade  spacing  for  transmission  to  the 
S2-hub-to-shroud  surfaces. 


The  Axisymmetric  Blade-to-Blade  Surface  Approximation 

A  considerable  simplification  of  the  whole  quasi-three  dimensional 
interaction  occurs  when  the  blade-to-blade  SI  streams urf aces  can  be 
considered  as  surfaces  of  revolution.  This  assumption  does  not  imply, 
however,  that  the  flow  is  to  be  considered  as  axisymmetric.  The 
blade-to-blade  equations  can  be  simplified  in  the  following  way,  as  a 
consequence  of  the  axisymmetry  and  of  additional  assumptions  with  regard 
to  the  tangential  uniformity  of  total  energy  I  and  entropy  s. 

i)  the  surface  n»n(0,m)  is  the  surface  of  revolution,  n«constant, 
generated  by  the  rotation  of  the  meridional  streamline  m,  see 
figure  2.2.3.  The  blade  force  reduces  to  a  force  in  the  n 

*4 

direction  according  to  equation  (2.2.108)  and  the  normal  n  to  the 
surface  is  in  the  direction  of  the  n-axis. 

Hence 


n 

m 


(2.2.109) 


ii)  In  a  steady  flow  the  total  energy  I  is  constant  along  a  streamline, 
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equation  (2.2.105)  and  hence,  the  total  energy  will  remain  constant 
in  the  whole  SI  flow  field  if  it  is  assumed  to  be  uniform  in  the 
,8-tangential  direction  in  the  inlet  section  of  the  blade-to-blade 
calculation  domain.  Thus,  if 

31  .......  .  (2.2.110) 

tt  “  0  in  the  inlet  surface 


one  has 


I  =  constant  in  the  whole  SI  surface 


(2.2.111) 


ill)  In  accordance  with  the  distributed  loss  model,  the  friction  force 
and  hence  the  entropy,  are  considered  as  an  average  representation 
of  the  losses  generated  by  the  shear  stresses.  Therefore  it  is 
consistent  to  assume  that  the  entropy  and  F^.  do  not  vary  in  the 

tangential  direction  0. 

This  corresponds  to  the  assumption  that  the  entropy  is  a  unique 
function  of  the  meridional  distance  m,  that  is 

s  -  s(m)  or  ||  -  0  (2.2.112) 


With  these  simplifications,  the  blade-to-blade  equations  for  an 
axisymmetric  streamsurface  and  tangentially  uniform  energy  and  entropy 
distributions  reduce  to  the  following  equations,  where  the  overbars  have 
been  removed  form  the  derivatives  as  a  consequence  of  (2.2.109) 

|j  (PB.)  .  |j  <pB,rWm)  *  |j  (pB,»e)  0  (2.2.113) 
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(2.2.115) 


31  +  31 

3t  +  Wm  3m 
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(2.2.116) 
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(2.2.117) 


A  further  simplification  is  introduced  for  steady  flows,  leading  to  the 
system 
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(2.2.118) 
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(2.2.119) 

(2.2.120) 

(2.2.121) 

(2.2.122) 


In  particular  for  steady  flows,  equation  (2.2.119)  shows  that  the 
n-component  of  the  absolute  vorticity  Cn  is  given  by 


„  ,  We  Ff  r  tan  B 

C  “  -r  F  /w  =■  r  —  —  *  -  .  F 

n  f0  m  w  w  w  f 

m 


(2.2.123) 


and  hence,  if  the  effects  of  the  friction  force  are  neglected  in  the 
blade-to-blade  surface,  the  absolute  flow  will  be  lrrotatlonal .  In  this 
case  the  absolute  velocity  will  depend  on  a  potential  function  <J>,  which 
satisfies  equation  (2.2.118),  in  the  steady  state  case. 

-  (2.2.124) 
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However,  it  is  to  be  noticed  that  this  assumption  is  not  consistent  with 
the  S2-hub-shroud  through-flow  since,  neglecting  the  friction  forces 
consists  in  assuming  constant  stagnation  pressure  and  hence  constant 
entropy  in  the  blade-to-blade  surface,  while  this  is  not  the  case  at  the 
corresponding  points  of  the  S2-surface.  The  introduction  of  the  losses 
leads  to  a  stagnation  pressure  variation  along  the  meridional  m-line  in 
the  through-flow  calculation,  while  the  above  assumptions  neglects  the 
same  variation  in  the  blade-to-blade  surface. 

In  addition,  the  application  of  the  isentropic  potential  equation 
(2.2.125)  to  the  internal  blade-to-blade  flows  is  faced  with  the 
non-uniqueness  problems  associated  with  transonic  potential  flows. 
Therefore,  equation  (2.2.125)  can  be  considered  as  a  good  approximation 
for  subsonic  and  mildly  supersonic  flows,  but  is  not  recommended  for 
choked  transonic  blade-to-blade  flows.  In  this  latter  case,  it  is 
advised  to  apply  the  equations  (2.2.113)  and  (2.2.114)  supplied  by 
equation  (2.2.116)  and  (2.2.117). 
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Another  option  consists  in  introducing  a  streamfunction  iKm,0),  defined 
by 


pB,w0  - 

the  vorticity  Cn  of  the 

1-  (J _  (_L_  li)  = 

30  pBir  30  3m  pBi  3m 
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3it 

30 


3 
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absolute  flow  can  be  written 


cn  *  fi(pBl“r!)  ■  r  r  BFr  * 
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as 

|^( PB,ujr 2)  (2. 2. 1 28) 


The  application  of  the  streamfunction  to  the  2D  blade-to-blade  flow, 
allows  to  take  into  account,  the  effect  of  the  friction  force  on  the 
inviscid  blade-to-blade  flow.  However  the  non-unique  relation  between 
density  and  streamfunction  makes  the  application  of  this  formulation 
difficult  for  transonic  flows. 


2.2.6.  The  Streamsurface  Quasi-Three-Dimensional  Interaction 


The  complete  Quasi-three  dimensional  interaction  can  now  be  described  as 
follows,  for  an  analysis  problem  whereby  the  geometry  of  the  machines  is 
given.  Detailed  considerations  with  regard  to  the  design  options  can  be 
found  in  Wu  (1 952) . 

i)  Define  a' series  of  S2-streamsurfaces  as  a  reasonable  approximation 
to  the  passage  flow  configurations.  A  valid  starting  option  is  to 
consider  only  one  S2-surface  assuming  an  initial  axisymmetric 
approximation . 

ii)  For  each  streamsurface  S2,  initial  approximations  have  to  be 
selected  for  the  tangential  blockage  or  sheet  thickness  coefficient 
B2 .  The  through-flow  is  computed  following  sections  2.2.1  to 

2.2.3. 

iii)  From  the  computed  S2-flows,  the  following  information  is  extracted 
in  order  to  define  the  various  SI  blade-to-blade  flows  : 

the  geometrical  coordinates  n-n(0,m)  of  the  selected 

blade-to-blade  surfaces, 

the  streamsheet  thickness  B!=Bi(0,m)  according  to  equation 

(2.2.101 ) 

-  the  body  force 

-  the  boundary  conditions  upstream  and  downstream  of  each  blade  to 
blade  section. 

iv)  The  various  SI  blade-to-blade  flows  are  computed  following  section 
2.2.5. 
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v)  From  the  computed  blade-to-blade  flows,  improved  S2-streamsurface 
coordinates  are  defined  as  well  as  blockage  coefficients  B2«B2(r,z) 

and  body  force. 

vi)  The  whole  process  is  repeated  until  a  satisfactory  convergence  is 
obtained. 

This  complete  iterative  procedure  has  been  applied  by  Krimerman  and 
Adler  (1978)  and  is  schematically  illustrated  in  figure  2.2.6 

The  iterative  process  can  be  quite  complex  and  lengthy,  especially  in 
configurations  where  the  blade-to-blade  surfaces  become  highly  twisted, 
but  it  allows  to  recover  the  full  three-dimensional  flow. 


Simplified  Quasi-Three-Dlmenslonal  Interaction  -  Mean  Streamsheet 

An  important  simplification  of  the  complete  process  consists  in  assuming 
that  the  SI -blade-to-blade  surfaces  are  surfaces  of  revolution  generated 
by  the  rotation  of  the  meridional  projection  of  the  streamlines.  This 
will  be  valid  when  the  blade-to-blade  surface  twist,  induced  by 
secondary  flows  remains  small.  In  addition,  it  is  considered  that  the 
through-flow  can  be  represented  by  a  single,  mean  streamsurface  S2, 
which  is  considered  to  be  representative  of  all  of  the  S2  surfaces. 
Hence  the  quasi-three  dimensional  flow  model  will  consist  of  one 
S2-through-flow  surface  and  several,  axisymmetric,  blade-to-blade 
surfaces . 

The  main  problem  in  this  simplified  approach  is  to  define  the  mean 
streamsheet.  Intuitively  it  is  to  be  considered  as  a  kind  of  average  or 
representative  surface  of  all  the  S2  surfaces.  However,  the  present 
approach  does  not  give  any  indication  as  how  to  select  this  surface. 

Various  choices  have  been  attempted  in  order  to  generate  the  mean  S2 
surface:  the  mean  camber  line,  Senoo  and  Nakase  (1972);  the 
streamlines  of  the  blade-to-blade  flows  separating  equally  the  mass  flow 
between  pressure  side  and  friction  side,  Novak  and  Hearsey  (1977);  the 
mass  averaged  blade-to-blade  streamlines,  these  ’'streamlines”  being 
different  from  any  actual  streamline,  Bosman  and  El-Shaarawi  (1977). 

No  indications  are  available  as  to  which  choice  is  to  be  recommended. 
As  discussed  by  Novak  and  Hearsey  (1977),  the  local  flow  values  can  be 
quite  sensitive  to  the  choice  of  the  distributions  of  the  6'  and  e' 
angles  defining  the  mean  streamsurface.  In  particular,  figure  2.2.7 
shows  a  comparison  of  the  meridional  variation  of  the  streamsurface 
angle  6’  in  function  of  meridional  distance  for  the  mean  camber  line  and 
the  equally  dividing  streamline  (4>=0.5  on  figure  2.2.8),  for  a  turbine 
nozzle  such  as  illustrated  in  figure  2.2.8,. 

Clearly,  the  strong  gradients  of  8'  along  the  camber  line  are  difficult 
to  accept  as  representative  of  the  mean  hub-shroud  flow.  In  addition, 
the  angle  e’  (lean  angle)  can  have  a  considerable  influence  on  the 
resulting  flow  distributions.  This  is  illustrated  in  figure  2.2.9  from 
Novak  and  Hearsey  (1977),  showing  a  comparison  of  two  calculations  of 
the  through-flow  along  an  S2-surface.  A  first  calculation  (Run  10) 
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corresponds  to  e'-O,  while  the  second  case  (Run  20)  had  a  mean 
streamsheet  lean  angle  e'  of  10  degrees  at  the  inner  radius  and  varying 
such  that  the  surface  was  straight.  The  computed  configuration  is  a 
turbine  nozzle  with  radius  ratio  of  0.88,  outlet  angle  of  72  degrees  and 
an  outlet  Mach  number  close  to  0.6.  In  the  second  run,  the  blades  were 
given  the  same  lean  as  the  streamsurface,  with  the  pressure  surface 
facing  towards  the  inner  radius.  As  can  be  seen  from  figure  2.2.9,  the 
radial  static  pressure  gradient  is  completely  altered  at  the  inner 
section  AA,  although  minor  differences  occur  at  the  trail  ng  edge. 

The  third  parameter  necessary  to  define  the  mean  streamsurface 
through-flow  is  the  thickness  parameter  B2.  As  already  observed  by  Wu 

(1952)  and  illustrated  by  Novak  and  Hearsey  (1977),  the  application  of 
equation  (2.2.28)  leads  to  variations  of  B2  with  very  strong  gradients 

as  shown  on  figure  2.2.9  for  the  same  computation  as  figures  2.2.7  and 
2.2.8.  In  addition,  for  subsonic  flows,  the  influence  of  the  blades 
extend  upstream  and  hence  modifies  the  value  of  B2  from  the  value  of 

unity  valid  in  bladeless  regions.  More  generally,  Horlock  and  Marsh 
(1971)  have  shown  that  no  choice  of  the  mean  streamsurface  can  be  made 
such  that  the  calculated  flow  along  this  mean  surface  would  represent 
exactly  the  actual  average  flow  over  the  passage  width,  see  also  section 
2.3.2. 

However,  the  overall  change,  from  leading  edge  to  trailing  edge  can  be 
correctly  modeled  if  the  streamsheet  angles  are  adapted  to  the  flow 
angles  at  these  locations.  But  there  is  no  guaranty  that  the  local 
values  of  the  flow  variables,  calculated  within  the  blade  passage,  are 
representative  of  the  actual  averaged  flow. 

Therefore,  although  the  simplification  provided  by  the  mean 
streamsurface  model  is  an  important  one,  its  precise  choice  can  lead  to 
uncertainties  as  to  the  representativity  of  the  calculated  local  flow 
behavior  inside  the  blade  rows  with  regard  to  the  real  passage  averaged 
flow. 

An  alternative  approach,  which  avoids  the  arbitrariness  of  the  choice  of 
the  mean  S2  surface,  can  be  defined  by  the  averaging  procedure  to  be 
discussed  next.  As  will  be  seen,  this  approach  which  solves  the  exact 
through-flow  equations  for  the  passage  averaged  flow  in  a  meridional 
plane  does  not  have  the  same  uncertainties  on  the  angular  parameters, 
8',  e'.  The  third  variable  B2  on  the  other  hand,  will  be  seen  to  be 

only  function  of  geometry,  see  equation  (2.3.5).  However,  the 
through-flow  equations  contain  additional  terms  depending  on  the 
blade-to-blade  flow  distributions,  which  have  to  be  estimated.  Since 
the  estimation  of  these  terms  can  be  derived,  if  necessary  from  the 
blade-to-blade  solutions,  one  has  a  general,  rigorous  approach  for  the 
simplified  quasi-three-dimensional  flow  model  consisting  of  the  passage 
averaged  flow  and  axisymmetric  blade-to-blade  surfaces.  The  averaging 
procedure  on  the  other  hand  cannot  be  applied  to  the  estimation  of  the 
full  three-dimensional  flow,  as  can  be  done  with  the  streamsurface 
approach. 
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2.3.  THE  PASSAGE-AVERAGED  QUASI-THREE  DIMENSIONAL  APPROXIMATION 


The  averaging  procedure  for  obtaining  a  quasi-three  dimensional 
approximation  to  the  turbomachinery  flow  consists  in  defining  a 
two-dimensional  flow  in  the  meridional  cross-section  of  the  machine, 
representative  of  an  average  flow  over  the  passage  width  of  the  blade 
row. 


Referring  to  figure  2.3.1.  the  averaged  flow  at  a  given  point  P  of  the 
meridional  cross-section  is  obtained  by  performing  an  integration  of  the 
flow  properties  from  the  pressure  surface  to  the  suction  surface  of  the 
next  blade  along  the  axisymmetric  cross  section  generated  by  the 
rotation  of  the  meridional  streamline  m.  In  order  to  apply  the  general 
theory  of  section  1.2,  a  cylindrical  coordinate  system  (r,0,z)  is  chosen 
and  the  integration  direction  The  limits  of  integration  are 
defined  by  the  blade  surfaces  and  the  width  b  of  the  integration  region 
is  given  by  b-e^-e^  where  0^  and  0g  are  the  angular  position  of  tne 

pressure  and  suction  surfaces  of  two  consecutive  blades.  If  N  is  the 
number  of  blades,  d  the  tangential  thickness  of  an  individual  blade,  and 
s  the  pitch  or  blade  spacing  at  radius  r,  one  has 
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(2.3.1) 


defining  a  streamtube  thickness  parameter  b2,  function  of  the  geometry, 

b2  =  1  -  d/s  (2.3.2) 


The  averaged  equations  have  the  following  form,  obtained  from  section 
1.2  (Table  3.1)  by  removing  the  shear  stress  and  heat  conduction  terms. 
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(2.3.5) 


(pE*b2)  ♦  $  ( pwl  b2 )  =  0 


A  body  force  f  appears  in  the  momentum  equation  as  a  result  of  the 

averaging  process  as  seen  in  section  1.2.  This  additional  force  f  is 
defined  by  equation  (1.2.37)  and  represents  here  the  resultant  of  the 
pressure  forces  on  the  blade  sections 
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where  n  and  n  are  the  normals  to  the  pressure  and  suction 
surfaces,  respectively. 

The  energy  equation  is  written  in  the  relative  system  with 


E* 


(2.3.7) 


The  above  equations  are  fully  rigorous. 


Various  averaged  quantities  can  be  defined  and  the  most  frequently 
occuring  choices  are  the  geometrical,  density  weighted  or  mass  flux 
weighted  averages.  A  consistent  set  of  equations  can  be  obtainedfor  the 
density  weighted  average  which  has  the  property  ofsatisfying  exactly 
mass  conservation,  Hirsch  and  Warzee  (1979). 

The  density  weighted  average  of  a  quantity  A  is  defined  by 
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where  A"  are  the  deviations  of  the  quantity  A  from  the  axisymmetric 
value  A. 

Applying  this  definition  with  the  general  procedure  of  section  1.2, 
leads  to  the  following  averaged  equations  for  the  through-flow  in  a 
meridional  cross  section  of  the  machine. 


2.3.1.  Passage  Averaged  Continuity  Equation 


The  averaged  continuity  equation  (2.3.3).  in  the  relative  system, 
becomes 
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where  and  are  the  density  weighted  averaged  relative  velocity 
components . 
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This  is  an  exact  equation  to  be  compared  with  equation  (2.2.35)  where  b2 

is  given  by  equation  (2.2.28).  Formally,  both  equations  are  identical, 
but  in  the  streamsurface  approach  b2  is  flow  dependent  and  is  only 

approximately  equal  to  the  value  given  by  equation  (2.3.2),  see  equation 
(2.2.56).  On  the  other  hand,  equation  (2.2.35)  is  the  mass  conservation 
law  for*  the  local  velocities  along  an  S2  streamsurface  which  has  to  be 
defined,  while  equation  (2.3.10)  is  the  mass  conservation  law  for  the 
averaged  meridional  flow  considered  in  a  true  meridional  cross  section 
of  the  machine. 

For  an  axisymmetric  flow,  however,  both  approaches  lead  to  exactly  the 
same  equations. 


2.3*2.  Passage  Averaged  Momentum  Equation 


Applying  the  results  summarized  in  equation  (1.2.39)  with  the  body  force 

f  defined  by  equation  (2.3.6),  the  momentum  equations  in  the  relative 
system  become 
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where  V  is  the  gradient  operating  on  the  r  and  z  coordinate  space.  That 
is,  the  explicit  components  of  the  gradient  operator  are  obtained  from 
the  three-dimensional  projections,  where  all  derivatives  with  respect  to 
0  are  put  to  zero. 

■+ 

The  normals  n  to  the  blade  surfaces,  figure  2.3.2,  appearing  in  the  body 
■+ 

force  term  f,  are  defined  by  the  angles  B'  and  e'  is  in  equations 
(2.2.16)  arid  (2.2.17)  and  are  given  respectively,  for  the  blade  pressure 
and  suction  surfaces,  by 
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Hence,  the  components  of  the  body  force  f  are 
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for  the  radial  component  and  for  the  axial  component 
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while  the  tangential  component  f  is 
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This  last  relation  shows  that  the  body  force  originates  from  the 
presence  of  the  blades  in  the  flow  field  as  solid  surfaces  able  to 
sustain  pressure  forces.  Actually  f.  is  nothing  else  than  the 

tangential  projection  of  the  local  lift  force  on  the  blade.  Hence, 

outside  the  blade  row  f  -  0  and  tne  flow  is  considered  as  axisymmetric. 


It  is  customary  to  separate  in  the  body  force  the  effect  of  the  blade 
camber  from  the  blade  thickness.  If  6^  is  the  angular  coordinate  of  the 

mean  camber  line,  one  has 
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Introducing  these  relations  in  the  definition  of  the  blade  force, 
equation  (2.3.6),  leads  to 
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The  first  term  pf  corresponds  to  the  lift  force  generated  by  the  mean 
camber  line  and  is  orthogonal  to  the  mean  camber  surface.  The  second 
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term  pf  is  a  correction  depending  on  the  radial  and  axial  variations 
of  the  blaie  thickness  blockage  factor  b2.  Explicitly,  one  has 
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The  last  term  of  equation  (2.3.12),  has  the  same  structure  as  the 
Reynolds  stres  terms  in  turbulent  averaged  Navier-Stokes  equations  and 
describes  the  influertces  of  the  full  three-dimensional  flow  on  its 
averaged  components.  We  will  define  a  secondary  stress  tensor  by 
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since  the  velocity  fluctuations  w"  are  generated  by  the  deviations  from 
axisymmetry,  in  particular  the  secondary  flow  contributions,  or  the 
blade-to-blade  flow  variations.  We  will  designate  the  gradients  of  the 
secondary  stresses  as  the  interaction  terms,  since  their  influence 
arises  essentially  from  the  interaction  betwen  the  averaged  flow  and  the 
secondary  flows.  These  terms  can  be  considered  as  the  sole 
contributions  from  the  three-dimensional  flow  distribution  to  the 
averaged  momentum  equation. 

An  alternative  form  to  equation  (2.3.12)  is  obtained  by  taking  into 
account  the  continuity  equation  (2.3*9),  leading  to 
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where  the  pressure  fluctuations  have  been  introduced  in  the  second  term. 

The  projections  of  the  momentum  equation  (2.3*21)  can  be  obtained  from 
the  general  transformation  laws  to  cylindrical  coordinates,  as  follows  : 


Radial  momentum  equation 
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where  the  secondary  stresses  are  represented  by  the  following  three 
contributions 
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Observe  that  the  non-derivative  velocity  terms  of  equation  (2.3.23)  can 
be  grouped  as 


V  (W  ♦  u>r)2  W  * 

— -  =  — — -  .  ——  +  oj2r  +  2ojW 

r  r  r  6 


(2.3.2b) 


Similarly  the  axial  component  of  the  momentum  equation  becomes 


It  “z  ‘  Wr  5F  “z  *  “z  h  \  :  £  '  Ffz  *  rBz  *  r-  t?(^<S)b.n2  (2.3.26) 

P  pt>2 


where  the  last  term  represents 


■(S) 


1  3 


1  3 


-[V(t  'b2)]z  *  -  (b2rpwJ)w|J)  ♦  -  —  (b2rpw|^w^) 


r  z  r  3z 


z  z 


(2.3.27) 


Tangential  component 


|r  Wq  ♦  W  f-  W  ♦  W  WQ  ♦  r  W  WQ 
3t  0  r  3r  6  z  3z  6  r  r  0 


’  -2“Ur  *  Ff»  *  fB6  * 


pb2 


T^b2)l 


(2.3.28) 


where 
‘(S), 


-[$( TvS,b2)j0  -  ;*  I-  (b2rpw"w")  +  ■■  |-  (b2rpw"w")  +  ]■  pw"w"b 


1  3 


(2.3.29) 


Here  also  the  non-derivative  velocity  terms  can  be  written  as 


l  WrlW0  ♦  2wr )  =  l  Wr(V0  ♦  wr) 


(2.3.30) 


A  detailed  investigation  of  the  influence  and  order  of  magnitude  of  the 
interaction  terms,  originating  from  the  blade-to-blade  flow  variations 
inside  a  compressor  blade  row,  has  b6en  presented  by  Smith  (1966).  The 
main  conclusion  is  that  the  interaction  terms  are  proportional  to  the 
blade  loading  or  lift  coefficient  but  that  their  effect  is  rather  small, 
although  their  influence  might  be  non-negligible  in  certain 
configurations  with  high  loadings  and  secondary  flows,  for  instance  in 
turbine  passages. 

A  similar  analysis,  performed  by  Hirsch  (1975),  for  the  contributions 
arising  from  the  two-dimensional  wakes  of  the  blades  shows  that  the 
interaction  terms  are  proportional,  outside  the  blade  passage  in  the 
vicinity  of  the  trailing  edge,  to  the  loss  coefficients  of  the 
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corresponding  blade  section.  The  influence  of  these  terms  has  been 
found  to  be  generally  negligible  although,  in  regions  of  high  losses 
such  as  the  sections  close  to  the  end-walls,  their  effect  might  need  to 
be  accounted  for. 

In  addition  to  the  two-dimensional  effects,  three-dimensional  effects 
due  to  secondary  flows  and  local  vorticity  production  can  strongly 
influence  the  averaged  flow,  through  significant  contributions  to  the 
interaction  terms. 

This  is  also  confirmed  by  an  experimental  investigation  in  a  transonic 
compressor  rotor,  Sehra  (1979).  It  was  shown  that  the  magnitude  of 
these  terms  has  a  non-negligible  effect  on  the  averaged  flow' in  the  end 
wall  regions  and  that  the  omission  of  their  contribution  led  to  an 
inacurate  prediction  of  the  radial  flow  distribution,  see  also  Sehra  and 
Kerrebrock  (1981).  A  recent  evaluation  of  the  relative  importance  of 
these  interaction  terms  has  been  presented  by  Jennions  and  Stow  (1985), 
for  turbine  vanes  with  strong  three-dimensional  flow  patterns. 

Hence  in  general,  the  contributions  of  the  secondary  stresses  will  have 
some  influence  on  the  averaged  flow,  mainly  in  the  endwall  regions  and 
in  the  main  blade  part,  unless  the  blade  rows  contain  lightly  loaded 
blades  or  generate  small  secondary  flows. 


2.3.3.  Passage  averaged  energy  equation 


The  averaged  energy  equation  is  obtained  from  the  general  form  given  by 
equations  (1.2.50)  to  (1.2.52)  written  for  a  relative  system.  Writing  I 

4  4 

instead  of  H  and  w  instead  of  v,  one  obtains,  within  the  assumptions  of 
the  distributed  loss  model,  equation  (2.1.9), 


|^-  (pE*b2)  +  ^  (pWIb2) 


-  $  (pw"I"b2)  ♦  b2Q. 


(2.3.31) 


where  the  additional  source  term  is  defined  by  equation  (1.2.52),  but 

takes  on  the  particular  form  of  equation  (1.2.53)  in  the  present  case  of 
a  turbomachinery  blade  passage. 


3b2 

Qb  ■  HT  [p  n  3  ■  5f-3B 


(2.3.32) 


In  this  relation,  f  is  the  body  force  defined  by  equation  (2.3.6)  or 
(2.3.15),  and  Ug  is  the  blade  displacement  speed  in  the  considered 
reference  system.  That  is,  for  a  rotor  in  the  relative  system 


(2.3.33) 
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and 


u  »  u)  x  r 


in  the  absolute  system. 

It  is  seen  from  the  above  equations  that  the  density  weighted  averaged 
total  energy  represented  by  the  averaged  rothalpy  I 


I 


-  H 


h  + 


2 


ul 

2 


pw" . w" 
2p 


I  + 


pw” . w" 
2p 


(2.3.3^) 


cannot  be  considered  as  constant  anymore  along  a  streamline  for  the 
steady  relative  flow.  Indeed,  in  the  relative  steady  state  situation, 
equation  (2.3.31)  becomes,  after  introduction  of  the  continuity  equation 


(piU)I  -  -  $(pw"I"b2) 

D  2 


or 


where 


Is  i  '  -  k 


W  T-“W 

m  9m  r  9r  z  9z 


(2.3.35) 


(2.3.36) 


(2.3.37) 


is  the  derivative  along  the  averaged  meridional  streamline  m.  This 
important  equation  (2.3*35)  shows  that  the  averaged  total  energy  varies 
along  a  streamline,  as  a  consequence  of  the  non-axisymmetric  energy 

fluxes  pw"I". 

These  energy  fluxes  are  due  to  the  transport  of  the  energy  fluctuations 
I"  by  the  non-axisymmetric  velocity  field  w" ,  where 


I" 


h" 


+  w" .  W  + 


-►  -► 
w"w" 

2 


pwftwH 

2p 


(2.3.38) 


represents  the  deviations  from  the  axisymmetric  energy  I. 
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One  of  the  main  contributions  to  these  energy  fluxes  will  come  from  the 
secondary  flows,  transporting  low  energy  fluid  from  the  end  wall  regions 
to  the  blade  surfaces  and  the  main  stream  region  of  the  blade  passage. 
These  effects  can  be  important  and  influence  the  overall  energy  exchange 
within  a  blade  row,  see  for  instance  Sehra  and  Kerrebrock  (1981),  Adkins 
and  Smith  ( 1 981 ) . 

The  energy  conservation  law  can  also  be  written  for  the  total  energy  of 
the  averaged  flow  I  defined  by  equation  (2. 3. 3*0.  One  obtains,  instead 
of  equation  (2.3.31) 

(plb2)  +  $(pW.Ib2)  =  -  $(pw"I"b2)  -  $( pW.kb2 )  ♦  b2Qb  (2.3-39) 


where  k  represents  the  average  kinetic  energy  of  the  large  scale 
fluctuations , 


pk 


pk" 


(2.3- ^0) 


Equation  (2.3*39)  shows  that  the  total  energy  of  the  density  averaged 

flow  does  not  remain  constant  along  an  averaged  streamline  following  W, 
in  the  steady  state  limit.  This  is  an  essential  limitation  of  the 
averaged  flow  model  for  turbomachinery  computation,  since  the  general 
procedure  in  practical  calculations  is  based  on  the  constancy  of  energy, 
in  conjunction  with  Crocco's  form  of  the  momentum  equation. 

2.3.11.  Crocco's  Form  for  the  Averaged  Momentum  Equations 


In  the  averaged  formulation,  Crocco's  form  can  be  derived  either 
directly  by  averaging  the  three-dimensional  form  of  equation  (2.1.8), 
Hirsch  and  Warzee  (1979),  or  by  introducing  the  entropy  and  the  enthalpy 
in  the  averaged  momentum  equation  (2.3.21),  following  section  1.2.1*. 
This  second  procedure  requires  the  averaging  of  the  entropy  relation 

p  TVs  =  p$h  -  $p  (2.3.111) 

The  averaged  pressure  gradient  gives  rise  to  the  body  force  as  seen  from 
section  1.2.  and  the  definition  (2.3.6).  Hence  one  has,  after 
averaging 


$p  -  pfD  -  -j-  $(pb2)  -  pf  -  ptfh  -  pT$s  (2. 3- !J2) 

O  D  2 


This  could  be  worked  out  following  section  1.2.11;  but  the  second  term 
on  the  right  hand  side  can  be  simplified  if  it  is  considered,  in 


-  92  - 


accordance  with  the  distributed  loss  model  assumptions,  that  the  entropy 

s  and  the  friction  force  F^.  are  averaged,  mean  quantities.  Hence  with 

s«s,  the  entropy  term  becomes  equal  to  pT$s.  The  enthalpy  term  is 
averaged  in  the  following  way,  applying  the  rules  of  section  1.2. 

p$h  -  g-  $(phb2)  ♦  Cph.n]  g  -  h$p  (2. 3. H3) 

The  last  term  can  be  decomposed  as  follows,  introducing  h  -  h  +  h" 

h.$p  -  £-  $pb2  +  h".$p  +  ~  Cpn]  (2.3.^^) 

D  2  D 


Combining  this  relation  with  equation  (2. 3. ^3),  gives  finally 


pVh  -  p'Jh  +  [ph".n]  g  -  h"$p 


(2- 3. H5) 


where  the  second  term  on  the  right  hand  side  is  the  new  expression  of 
the  body  force 


pb2f(h)  *  -[ph".n] 


2ir/N 


(2. 3.  ^6) 


Equation  (2.3*  *12)  becomes 


$p  -  pf_  -  —  $pb2  -  pf  -  p^h  -  pf'h'  -  h"$p  -  pT$s 

D  D  2 


(2. 3. 47) 


An  additional  interaction  term  is  introduced  here,  namely  h"Vp  which  can 
be  approximated  in  practical  computations  by 


h"$p  -  h"$p  +  h"$p ’  «  h"$p 


(2-3- ^8) 


since  the  second  term  can  be  expected  to  be  of  a  lower  order  of 
magnitude. 

Introducing  equation  (2. 3* **7)  in  the  averaged  momentum  equations  leads 
to  the  following  form  of  the  averaged  Crocco's  momentum  equation 

f-  V$  -  W  x  c  «  T$s  -  $1  +  +  ?(h) - —  $(pw"  x  w"b2 )  ♦  (2. 3.  ^9) 

b2p  p 
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where  the  averaged  absolute  voriticity  is  defined  by 


C  -  V  x  V 


(2.3.50) 


The  total  energy  gradient  term  appearing  in  this  equation  contains  the 

A 

total  energy  of  the  averaged  flow  I.  In  practical  computations,  where 

~  ■+ 

the  averaged  quantities  h,  W  are  calculated,  this  quantity  is  the 
obvious  one  to  consider,  since  it  is  only  dependent  on  the  averaged  flow 

variables.  This  is  not  the  case  for  the  averaged  total  energy  I  which 

depends  also  on  the  large  scale  flow  deviations  from  axisymmetry  w". 
However,  neither  of  these  quantities  are  conserved  along  a  streamline  of 
the  density  averaged  flow  in  steady  state  conditions.  The  implications 
of  this  situation  have  been  discussed  by  Sehra  (1979)  with  respect  to 
the  interpretations  of  overall  radial  equilibrium  as  described  by 
equation  (2.3.49)  and  more  recently  by  Hirsch  and  Dring  (1986). 

Compared  to  the  corresponding  momemtum  equation  (2.2.5)  of  the 
streamsurface  approach,  similar  terms  are  obtained  here,  with  the 
addition  of  the  already  discussed  interaction  terms.  The  projections  of 
equation  (2.3.49)  follow  closely  the  formulas  given  in  section  2.2.1, 
with  the  projections  of  the  interaction  terms  calculated  in  section 
2.3.2.  The  particular  projection  defined  by  equation  (2.2.73)  is 
obtained  in  the  same  way  but  the  projection  of  the  secondary  stress 

+ 

gradients  on  the  direction  N,  has  to  be  added. 


Alternative  formulation  of  the  averaged  energy  equation 

By  introducing  the  total  energy  fluctuations  (2.3-38)  into  the  energy 
equation  (2.3.39),  one  obtains 


|jr  (plb2)  ♦  $(pWIb2)  =  |jr  (pb2 )  -  $(pw"h"b2)  -  $(pwk"b2) 


(2.3.51 ) 


+  $(T(3).Wb2)  ♦  b2U. 


One  can  notice  the  appearance  of  a  contribution  from  the  work  of  the 
secondary  stresses  against  the  averaged  flow,  in  addition  to  terms 
describing  a  diffusive  effect  of  the  fluctuating  energies  h"  and  k" . 

Observe  also  that  the  term  V(pwk"b2)  contains  the  full  velocity  vector  w 
-  W  +  w". 

In  steady  state  conditions,  the  above  equation  becomes 
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(2.3.52) 


$( pWIb2 )  -  -$( pw"h"b2 )  -  $(pwk"b2)  +  $(x3Wb2)  ♦  b2Qb 


The  secondary  stress  term  can  be  decomposed,  in  analogy  with  the  viscous 
shear  stress  contributions  in  a  pseudo  dissipation  term  (t  .v)W  and  an 

4  4  «g 

energy  producing  term  (W.v)t  .  Sehra  and  Keerebrock  (1981)  postulated 

that  the  former  term  plus  the  kinetic  energy  term  $(pwk"b2)  were  lost  to 

the  main  flow  and  hence  dissipated,  leading  to  an  irreversible  entropy 
production . 

However,  there  is  presently  too  little  experimental  evidence  to  confirm 
or  infirm  this  hypothesis. 


2.3.5.  Blade-to-Blade  Flow  Equations  in  the  Averaged  Flow  Model 


Since  the  through-flow  is  represented  by  an  average  flow,  the  only 
consistent  assumption  with  regard  to  the  blade-to-blade  flows  consists 
in  assuming  the  corresponding  streamsurfaces  to  be  surfaces  of 
revolution . 

In  (0,m,n)  coordinates,  the  blade-to-blade  flow  equations  are  obtained 
by  considering  the  tangential  component  of  the  momentum  equations  for  a 
streamsheet  of  thickness  B,-  Bi(m).  Therefore  one  obtains  exactly  the 

same  equations  as  described  in  section  2.2.5. 


2.3.6.  Quasi-Three  Dimensional  Interaction  for  the  Averaged  Flow  Model 


The  only  quasi-three  dimensional  interaction  possible  here  is  the 
simplified  interaction  whereby  the  role  of  the  mean  streamsheet  is 
played  by  the  averaged  flow  considered  in  a  true  meridional  plane. 

The  blade-to-blade  streamsheet  thickness  B,  as  well  as  the  location  of 

the  meridional  streamlines  m  are  an  output  of  the  through-flow 
computation  of  the  averaged  flow.  The  averaged  flow  is  first  computed 
with  an  initial  approximation  whereby  all  the  secondary  stresses  are  set 
to  zero,  that  is,  an  axisymmetric  flow  approximation.  This  delivers 
also  the  boundary  conditions  on  the  different  blade-to-blade  sections. 
After  the  computation  of  the  blade-to-blade  flows  an  estimation  of  the 
interaction  terms  can  be  obtained  from  the  blade-to-blade  flows 
distribution,  but  also  from  additional,  direct  estimations  of  secondary 
flows  and  streamwise  vorticity,  Hawthorne  (1955).  A  new  iteration  cycle 
can  then  be  initiated  for  the  through-flow  including  the  interaction 
terms . 

This  procedure  has  been  applied,  and  developed  into  a  prediction  and 
design  system,  by  Jennions  and  Stow  (1985),  who  performed  also  a 
detailed  investigation  of  the  comparative  influence  of  the  different 


-  95  - 


contributions  to  the  radial  equilibrium  equation,  including  the 
interaction  terms. 


Relation  between  streamsurface  and  averaged  formulations 

The  question  arises  actually  in  the  description  of  turbomachinery  flows 
as  to  the  equivalence  of  the  averaged  flow  representation  and  the  meaan 
streamsheet  representation  used  in  the  simplified  quasi-three 
dimensional  flow  description  discussed  in  the  previous  section.  This 
topic  has  been  discussed  in  particular  by  Horlock  and  Marsh  (1971). 

If  a  mean  streamsheet  is  defined  by  angles  S'  and  e',  is  there  a  choice 
of  values  of  8'  and  e'  which  will  render  the  two  representations 
identical  ? 


Comparing  the  equations  of  the  two  models,  a  first  conclusion  arises 
from  the  energy  equations  (2. 2. 5*0  and  (2.3.35).  It  is  clear  that 

A 

neither  the  total  energy  of  the  averaged  flow  I,  nor  the  averaged  total 

energy  I,  are  able  to  represent  the  total  energy  along  the  mean  S2 
streamsurface.  For  a  steady  state  flow,  equation  (2.2.5^)  gives  with 
the  definitions  (2. 2.1^),  (2.2.15)  and  the  orthogonality  relation 

(2.2.19) 


31  /  31 x 

m  3m  '  e  36; 


(2.3.53) 


while  equation  (2.3*36)  leads  to 


31 

p  m  3m 


-  ClF 


(2.3.510 


Unless  the  right  hand  sides  of  both  equations  are  zero,  which  is  the 

case  in  an  axisymmetric  flow,  no  general  choice  of  I  would  render  the 
energy  fluxes  of  the  interaction  terms  zero.  With  regard  to  the  body 
forces  in  the  two  models,  a  comparison  has  to  take  into  account  the 
contributions  represented  by  equation  (2.3.22). 

If  the  pressure  variation  between  pressure  and  suction  surface,  at  a 
constant  radial  and  axial  position,  is  linear,  equation  (2.3.22)  shows 

that  the  force  pf^  becomes  identical  to  the  corresponding  term  in 

equation  (2.2.8)  if  the  mean  S2  streamsurface  is  selected  to  be 
identical  with  the  camber  line. 

Indeed,  with  this  assumption  the  body  force  equation  (2.2.30)  becomes 

'  '  "<H>  '  (2.3.55) 


-  96  - 


and  is  identical  to  if  &'  -  and  e'  ■  e'^. 

However,  for  a  more  general  pressure  variation  law,  the  two  body  forces 
will  not  be  equal,  unless  the  blade  thickness  is  constant  or  zero.  This 
corresponds  to  an  infinitely  thin  blade  and  a  mean  streamsheet  following 
the  blade  surface. 

In  addition,  as  discussed  previously  the  interaction  terms  in  the 
momentum  equations  are  non-negl igible  except  for  lightly  loaded  bl  ides 
generating  negligible  secondary  flow  effects.  This  would  be  the  'ase 
for  instance  for  highly  spaced,  low  cambered  cascade  blades  or  with 
thin,  many  bladed  cascades. 

Some  direct  estimations  in  simplified  cases  of  the  different 
contributions  to  various  flow  representations  can  be  found  in  Horlock 
and  Marsh  (1971). 

It  is  to  be  concluded  from  these  considerations,  that  in  general,  it  is 
not  possible  to  define  a  mean  S2  streamsurface  on  which  the  flow 
variations  are  identical  to  the  averaged  flow  variations.  Considering 
that  the  averaged  flow  model  describes  correctly  the  behavior  of  the 
real  physical  average  flow  over  the  blade  spacing,  the  flow  on  a  mean 
streamsurface  will  not  be  able  to  represent  correctly  the  averaged  flow. 

Only  in  the  hypothetical  case  of  an  axlsymmetrlc  flow,  will  both  models 
be  identical .  Therefore,  the  interpretat ion  of  the  calculated  local 
behavior  of  the  flow  variables  along  a  mean  S2  streamsurface  should  be 
taken  with  caution  when  considering  them  as  representat ive  of  the 
averaged  actual  flow. 
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2.4.  THROUGH-FLOW  COMPUTATIONAL  METHODS 


Several  methods  have  been  developed  for  the  numerical  resolution  of  the 
through-flow  equations  in  an  arbitrary,  multistage,  turbomachine  and 
they  can  be  classified  in  two  broad  families.  The  first  family  is  based 
on  the  solution  of  the  momentum  and  continuity  equations  for  the 
physical  variables,  while  the  second  family  is  based  on  the  introduction 
of  a  streamf unction  in  the  meridional  coordinates  (r,z)  leading  to  a 
second  order  partial  differential  equation  for  the  streamf unction. 

Within  the  first  family,  the  streamline  curvature  method  is  based  on  an 
equation  for  the  meridional  velocity  vPf-v^y,  along  a  given  calculation 

station,  obtained  from  the  stationary  momentum  equation.  The  dependence 
of  the  flow  on  the  Other  coordinate  is  expressed  through  the  curvature 
of  the  meridional  projection  of  the  streamlines;  hence  the  name  of  the 
method.  The  streamline  curvature  method  was  developed  for  through-flow 
computations  in  its  contemporary  form  in  the  early  sixties,  Katsanis 
(1964),  (1966),  Smith  (1966),  Novak  (1967),  Jansen  and  Moffatt  (1967) 
and  is  still  widely  used. 

The  streamf unction  approach  was  introduced  by  Wu  (1952)  in  a  finite 
difference  discretization  and  renamed  matrix  through-flow  method  by 
Marsh  (1968).  The  same  approach  has  been  followed  by  Davis  and  Millar 
(1972),  (1976),  Davis  (1975),  Biniaris  (1975),  Bosman  and  Marsh  (1974), 
Bosman  and  El  Shaarawi  (1977). 

The  discretization  of  the  same  streamfunction  equations  by  finite 
elements  has  been  introduced  independently  by  Adler  and  Krimerman  (1974) 
and  Hirsch  and  Warzee  (1974),  (1976). 

The  various  methods  will  be  discussed  in  this  section  from  the  point  of 
view  of  the  numerical  approach,  leaving  aside  any  detailed 
considerations  with  regard  to  the  necessary  empirical  input  under  the 
form  of  loss  coefficients  and  turning  angles  as  discussed  in  section 
2.2.1 . 


2.4.1.  Axlsymmetrlc  Through-Flow  Equations 


The  current  approximations  in  through-flow  computations  assume  relative 
steady  flow  situations  at  constant  wheel  speed  u>r  and  some  form  of 
axisymmetry  of  the  flow  configuration.  Actually,  within  the 
streamsurface  approach,  the  equations  along  the  S2  streamsheet  are 
formally  identical  to  the  axisymmetric  flow  equations  since  the 
replacement  of  the  streamsurface  derivatives  by  ordinary  derivatives 
leads  to  the  axisymmetric  equations.  Hence,  within  the  streamsurface 
approach,  the  distinction  between  the  flow  along  an  arbitrary  S2 
streamsheet  and  an  axisymmetric  flow  is  a  matter  of  interpretation,  the 
numerical  resolution  techniques  being  identical  in  both  cases. 

This  is  not  the  case  within  the  passage  averaged  model,  since  the  extra 

i  *S 

terms  represented  by  the  secondary  stresses  v.t  in  the  momentum 
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equations  and  by  the  additional  energy  exchange  and  dissipation  terms  in 
the  energy  equations,  vanish  only  in  the  axisymmetric  assumption. 
Within  a  complete  quasi-three-dimensional  computation,  as  discussed  in 
section  2.3.6,  the  first  step  of  the  calculation  will  be  based  on  an 
axisymmetric  through-flow,  while  the  next  iterations  will  introduce  some 
estimations  of  the  interaction  terms  between  the  averaged  flow  and  the 
non-axisymmetric  contributions,  and  take  into  account  their 
contributions  to  the  momentum  equations  as  well  as  on  the  energy  and 
entropy  equations.  Since  the  basic  approximation  is  based  on  the 
axisymmetric  assumption,  we  will  consider  in  the  following  all  methods 
as  applied  to  the  axisymmetric  through-flow  models. 

As  shown  in  the  previous  sections,  the  through-flow  equations  consist  of 
the  steady  state  form  of  the  mass  conservation  equation,  one  of  the 
meridional  (axial,  radial  or  any  combination  such  as  the  N-component) 
projections  of  the  momentum  conservation  law,  the  tangential  momentum 
equations,  the  energy  donservation  and  the  entropy  equation  connecting 
the  stagnation  pressure  variation  to  the  empirical  loss  coefficients. 

The  through-flow  equations  can  be  summarized,  within  the  axisymmetric 
assumption  as  follows,  in  cylindrical  coordinates 


Continuity  equation 


(2.4.1) 


Radial  equilibrium  equation 


w  -r—  w 
m  3m  r 


1  £P 
p  dr 


(2.4.2) 


Alternate  form 


(2.4.3) 


Axial  momentum  component 


m  dm  z 


(2.4.4) 
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Alternate  form 


,9  9  n 

W  (t—  w  -  —  w  ) 
r  9r  z  9z  r 


9s 
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9z 
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’  9z  U  9  fz  Bz 


(2.4.5) 


Tangential  momentum  component 


l 

r 


“m  3m  (rv9) 


fe 


+  f 


B9 


(2.4.6) 


Energy  equation 


pw 


91 
m  9m 


0 


(2.4.7) 


Entropy  equation 


The  friction  force  F_=-F„1 

I  I  w 

_  9s 
Tw  t —  m 
m  9m 


is  defined  by 


(2.4.8) 


and  the  meridional  convection  operator  is  defined  by 

9  9  9 

w  r—  -  w  —  +  w 
m  9m  r  9r  z  9z 


(2.4.9) 


The  body  force  is  normal  to  the  streamsurface  S2,  in  the  streamsurface 
approach,  or  to  the  mean  camber  line  in  the  averaged  formulation,  under 
the  axisymmetric  assumption,  that  is 

?B  =  rfgii  (2.4.10) 


according  to  equations  (2.2.30)  to  (2.2.32)  or  equations  (2.3-22)  in  the 
appropriate  interpretation  with 

n  n 

—  -  -  tan  ef  —  -  -  tan  6f  (2.i4.11) 

n^  n^ 


and 
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(2.4.12) 


w  =  0 


For  the  averaged  flow  model  it  is  assumed  that  the  flow  is  tangent  to 
the  mean  camber  line,  according  to  equations  (2.3*18)  and  (2.3*22). 

The  above  equations  are  valid  outside  the  boundary  layer  regions,  since 

the  friction  force  Fj.  is  not  resolved  and  is  related  to  overall 

empirical  loss  coefficients.  However  the  complete  through  flow  is 
influenced  by  the  end  wall  boundary  layers,  through  the  blockage  effect 
introduced  by  their  displacement  thickness  and  through  their  effect  on 
the  overall  energy  exchange  process,  Smith  (1969). 

The  passage  averaged  equations  allow  in  addition  to  define  a  coherent 
interaction  between  the  meridional  through-flow  and  the  end-wall 
boundary  layers  developing  along  the  hub  and  shroud  walls  of  a 
turbomachine,  figure  2.4.1. 

A  consistant  theory  describing  the  influence  and  interaction  of  the  end 
wall  regions  and  the  main  flow  has  been  presented  by  Mellor  and  Wood 
(1971),  and  further  extended  by  Hirsch  (1974),  (1976),  Horlock  and 
Perkins  (1974),  De  Ruyck,  Hirsch  and  Kool  (1979),  De  Ruyck  and  Hirsch 
(1981),  (1983). 

Within  the  passage-averaged  flow  representation,  the  mainstream  through 
flow  will  be  described,  for  instance  by  the  radial  component  of  the 
momentum  equation  in  an  axial  flow  machine,  in  addition  to  the 
continuity  equation  and  the  energy  and  entropy  laws.  Within  the  same 
representation,  the  end  wall  boundary  layer  equations  will  make  use  of 
the  other  two  components  of  the  momentum  equations. 


2.4.2.  The  Streamline  Curvature  Method 


The  streamline  curvature  method  is  based  on  the  transformation  of  the 
"meridional"  momentum  equation,  that  is  either  the  fadial  or  axial 
component  or  the  projection  of  the  momentum  equation  along  any  direction 
in  the  (r,z)  surface,  into  an  ordinary  differential  equation  for  the 
meridional  velocity  w^  in  function  of  the  distance  along  selected 

calculation  stations  joining  corresponding  points  on  the  hub  and  shroud 
surfaces . 

An  arbitrary  calculation  station  is  displayed  on  figure  2.4.2,  joining 
two  points  A  and  B  on  the  hub  and  shroud  surfaces.  The  station  may  be 
curved  and  is  defined  by  the  local  angle  Y  with  respect  to  the  local 
radial  direction.  The  angle  Y  is  defined  as  positive  from  the  radial 
axis  towards  the  station  AB.  Hence  Y  is  positive  on  figure  2.4.2. 

The  distance  along  the  calculation  station  AB  is  denoted  by  1  and  an 
equation  is  sought  for  the  variation  dw^/dl  along  that  station  AB. 

The  whole  computational  procedure  consists  in  sweeping  through  the 
meridional  Section  of  the  machine,  from  a  calculation  station  at  inlet 
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of  the  machine  towards  the  outlet  station,  solving  for  the  radial 
distribution  of  flow  properties  along  each  station.  The  meridional 
velocity  distribution  along  a  calculation  station  is  adapted  in  order  to 
satisfy  mass  flow  conservation  along  each  station  and  the  whole 
procedure  is  repeated  until  convergence  is  reached.  The  radial  flow 
distributions  at  different  stations  have  to  be  consistent  with  each 
other  and  the  normal  acceleration  which  a  fluid  particle  undergoes  has 
to  be  taken  into  account  as  it  influences  the  interaction  of  flow 
properties  between  successive  stations.  This  is  obtained  through 
assumptions  for  the  meridional  streamline  shape  and  its  slope  and 
curvature  as  it  interconnects  the  different  stations.  The  details  of 
the  procedure  are  given  in  the  following. 

The  basic  equation  for  dw^/dl  can  be  obtained  in  a  variety  of  ways  and 

written  in  many  different,  but  equivalent  forms,  Katsanis  (196^), 
(1966),  Smith '( 1 966) ,  Novak  (1967),  Wennerstrom  (197H). 


+T 

Defining  a  total  force  vector  F 


ft 


F  +  f 
f  B 


(2.i».13) 


the  radial  equilibrium  equation  becomes 
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(2.4.15) 


while  the  axial  component  takes  the  form 
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d_ 
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or 


,d  d  , 

w  (t—  W  -  —  W  ) 
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T  3s  di  we  d  ,  . 

T  —  -  —  +  —  —  (rv  ) 
3z  3z  r  dz  d 


+  F 


(2.11.17) 


If  the  complete  form  of  the  passage  averaged  equations  are  solved,  F 
will  be  defined  by  adding  the  secondary  stress  term  to  equation 
(2.11.13).  see  equations  (2.3.21)  to  (2.3.30). 


Momentum  equations  along  calculation  stations 

Referring  to  figure  2.11.2,  the  derivatives  with  respect  to  (r,z)  can  be 
replaced  by  derivatives  with  respect  to  m  and  1  along  the  meridional 
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streamline  and  the  calculation  station  AB. 


One  has,  with  Y  positive  when  1  is  located  left  of  the  radial 
on  figure  2.4.2, 
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The  radial  and  axial  momentum  equations  are  combined  by 
equation  (2.4.14)  by  cosY  and  equation  (2.4.16)  by  (-sin  Y) 
up,  applying  equation  (2.4.18) 


9  3 

w  (cos  Y  r—  w  -  sin  Y  7—  w  ) 
m  3m  r  3m  z 


- 1  If  *  r  <=os  y  . 


where 


T  T  T 

F.  -  cos  Y  F  -  sin  Y  F 
1  r  z 


is  the  projection  of  the  total  force  F  on  the  1  direction. 
With 


w  »  w  sin  o 
r  m 


one  has 


.  -  W  .  cos  0 

d  .  d  m 
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3m  r  3m  m  R 

m 

where  the  radius  of  curvature  of  the  meridional  streamline  R 

i 

introduced  according  to  equation  (A. 6. 40) 

2l  ,  3o 

R  3m 
m 


axis,  as 

(2.4.18) 

(2.4.19) 

(2.4.20) 

(2.4.21 ) 

multiplying 
and  adding 

(2.4.22) 

(2.4.23) 

(2.4.24) 

(2.4.25) 

has  been 

i 

(2.4.26) 
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Note  that  with  this  definition  Rm  is  negative  for  a  concave  streamline, 
such  as  shown  on  figure  2.4.2. 


Similarly 
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leading  to 
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(2.4.27) 


(2.4.28) 


This  equation  has  to  be  coupled  to  the  energy  equation  in  order  to 
obtain  an  equation  for  This  is  achieved  explicitely  by  using 

Crocco's  form  of  the  momentum  equations,  equations  (2.4.15)  and 
(2.4.16). 


The  same  procedure  leads  to  the  equation  for  wm=wm(D  along  the 
calculation  station  AB  : 


m 


3_ 

31 


w- 


m 


sin(o-Y)w  — 
m  3m 


*  -  COS(O-Y)^  -  T  If 

m  R  31 

m 


31 

31 


w 


r  |r(rve)”FI(2*1**29) 


The  estimation  of  the  gradients  of  s,  rv.  and  I  is  performed  according 

O 

to  the  discussion  of  section  2,2,1.  Note  that  the  particular  case  T-0 
corresponds  to  a  radial  station,  while  T-±tt/2  is  an  axially  oriented 
station.  For  a  radial  station,  one  obtains 
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(2.4.30) 


Outside  bladed  regions,  the  force  components  vanish  as  well  as  the 

meridional  entropy  and  total  energy  gradients  because  of  the  assumption 
Of  axi3ymmetry.  However,  if  wake  effects  are  taken  into  account  these 
terms  might  not  be  zero  and  have  to  be  estimated  through  the 
introduction  of  appropriate  wake  flow  models,  Hirsch  (1975),  Hearsey 
(1975). 

When  the  computing  station  AB  is  orthogonal  to  the  meridional 

streamline,  that  is,  Y»o,  the  first  term  on  the  right  hand  side  of 

equation  (2.4.29),  describing  the  streamwise  acceleration  due  to  the 

meridional  variation  of  the  meridional  velocity  component  w  ,  vanishes. 

m 

However,  for  arbitrary  calculation  stations,  this  will  not  be  the  case 


and  this  term  has  to  be  evaluated.  In  early  through-flow  calculations 

in  axial  flow  machines,  with  calculation  stations  located  only  outside 

the  bladed  regions  this  term  was  often  neglected.  In  more  general  cases 

and  when  calculation  stations  are  located  inside  blade  rows,  this  will 

not  be  Justified  and  the  term  3  w  has  to  be  taken  into  account.  So 

m  m 

more  that  its  variation  might  become  important  at  transonic  Mach  number 
ranges.  Indeed,  equation  (2.2.84)  shows  that  the  variation  of  mass  flux 
pw^  is  related  to  the  variation  of  wm  by  the  following  equation  valid 

for  lsentroplc,  lsoenergetic  conditions,  H-constant,  in  an  absolute 
reference  system 
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m  mm 
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Hence  an 

alternative  form  of  the 

meridional 

velocity  equation 

is 

obtained 

from  an  explicit  evaluation 

Of  jp  w  . 

3m  m 

Evaluation  of  meridional  streamwlse  acceleration  -  Design  problem 


3w 

In  order  to  evaluate  the  term  the  continuity  equation  is  introduced 

with  equations  (2.4.20)  and  (2.4.21),  leading  to 
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0  (2.4.32) 


Since,  in  general  the  turbomachinery  flow  is  not  isentropic  because  of 
the  friction  force  F^.  representing  the  stagnation  pressure  losses, 

equation  (2.2.82)  has  to  be  replaced  by  the  more  general  relation, 


6p  6h  6s 

p  a"5"  R 


(2.4.33) 


where  6h  is  obtained  from  the  energy  equation  (2.4.7) 

Here  R  has  been  used  for  the  gas  constant  in  order  to  avoid  confusion 
with  the  radius  r. 


Equation  (2.4.31)  is  to  be  replaced  by  a  more  general  expression 
obtained  from  (2.4.33)  where  the  definition  of  I  has  been  used.  Since 
the  right  hand  side  of  equation  (2.4.31)  corresponds  to  the  first  term 
of  (2.4.32),  the  density  variation  in  the  m-direction  becomes 
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where  the  relative  Mach  number  component  MQ  is  introduced 
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(2.^35) 


with  the  definition  sin  o  -  3r/3m. 

Introducing  equation  (2.44. 3*4)  into  equation  (2.44.32)  leads  to  the 
following  expression  for  the  acceleration  term  3wm/3m 

u2.  sin  o  1  3o  tan(o-Y)1 

a?  r  cos(o-Y)  31  R  J 

m  (2. H.36) 


The  meridional  derivative  of  the  tangential  relative  velocity  can  be 
eliminated  through  the  use  of  the  tangential  momentum  equation  (2.44.6), 
written  as 
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T 

where  F„  will  generally  be 
o 

the  design  problem,  rv 
analysis  case  v  has  to  be 
and  the  specified  flow 
becomes 


taken  to  be  zero  outside  the  blade  rows.  In 

is  specified  within  a  blade  row,  while  in  the 

expressed  in  function  of  meridional  velocity 

angles.  The  3w  /3m  term  in  equation  (2.H.36) 
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leading  to  the  following  form  for  the  meridional  variation  of  w 
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Introduced  in  equation  (2.44.29),  the  basic  equation  written  as  an 
ordinary  differential  equation  for  wm“wm(D  along  the  station  AB,  is 

obtained 
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The  equations  (2.4.40)  to  (2.4.42)  are  to  be  used  outside  blade  rows 
with 
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outside  blade  rows 


(2.H.1J3) 


and  also  inside  blade  rows  with  a  design  option  whereby  rv.  is  specified 

as  the  input  information  from  the  blade-to-blade  flow.  Note  also  the 
appearance  of  the  (1-M*)  factor  which  is  an  expression  of  the  range  of 

ellipticity  of  the  through-flow  equations  as  discussed  in  section  8.2, 
when  rvQ  is  specified. 

When  flow  angles  are  specified,  -  the  analysis  case  -  the  terms 
containing  derivatives  of  tangential  velocity  components  have  to  be 
explicitely  calculated  in  function  of  meridional  velocity  and  flow 
angle . 


Analysis  Problem  -  Flow  angles  specified 

With 
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Within  blade  rows,  the  full  dependence  of  3wm/3m  on  the  other  flow 

variables,  in  particular  the  flow  angle  6,  becomes  instead  of  equation 
(2.4.36) 
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Introducing  equation  (2.4.47)  in  equation  (2.4.29)  and  replacing  the 

g 

term  w  /r  rr-(rv  )  with  equation  (2.4.46),  leads  to  the  following  form 
0  ol  0 

for  the  momentum  equation 


A, 


(1 


+  tan  B  tan  cOw^  rr  w„ 
m  dl  m 


A,  w2 
1  m 


+  B, 


(2.4.48) 


( r ) 2 

1-M^  cos2(o-Y) 

(1-M(r)  )R  cos(o~7) 
m 


(i 


S>* 


1  Mi 

R  3mJ 


sln(o-y)  j-l_  3b  + _ ] _  3o  +  sin  o 

,  M(r)2  Lb  3m  cos(o-Y)  31  r 
.  (2.4.49) 


B, 


31  3s 

31  31 


(2.4.50) 


(r)2 

Note  the  appearance  of  the  factor  (1-M  )  in  the  denominator.  This 
form  is  limited  to  problems  with  fully  subsonic  relative  Mach  numbers 
(or  fully  supersonic). 

In  practical  computations  inside  blade  rows,  it  might  be  more 
appropriate  to  calculate  Sw^/Sm  directly  by  a  finite  difference 

approximation,  using  values  from  the  previous  iteration  and  apply 
equation  (2.4.29)  with  equation  (2.4.46)  introduced  for  the  last  term 
when  dealing  with  the  analysis  case. 


For  incompressible  flows,  equation  (2.4.32)  becomes 
sin  o 


jn  3j) 
b  3m 


m 


1  r  3°  ^  sin(o-Y)  ,  A  3 

cos  ( o-7)  L  di  31  R  mJ  3m  m 

m 


(2.4.51 ) 


while  equation  (2.4.28)  can  be  transformed  to  the  following  form  after 
the  introduction  of  the  rotary  stagnation  pressure 


P* 


v2 

p  +  p  ~2  -  puvQ 


P  + 


(2.4.52) 
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w  *  w 

J  -  ir  cos(o-Y)  -  sin(o-Y)w  Ir  W  ♦  w  L  w  ♦  -5  |r(rv  )  ♦  fJ  (2.U.53) 
p  3  1  R  m  9m  m  m  31  m  r  31  9  1 

K  m 

Combining  both  equations,  one  obtains  the  ordinary  differential  equation 

for  w  -w  (1)  valid  for  incompressible  flows 
m  m 


W  -rrr  W  -  A  W2  +  B 

m  31  m  m 


(2.H.5H) 


with 

A  - 


R  coS(o-Y) 
m 


sin(Y+o)  [~  + 


1 


b  3m  cos(o-Y)  31 


3°  ♦  gAL.0-] 


(2.14.55) 


n  1  3p*  63,  . 
B  “  p  31 - r  3l  rv9 


(2. 14 . 54) 


The  above  equations  are  the  basic  formulation  for  streamline  curvature 
methods,  coupled  to  the  mass  flow  conservation  law  in  integral  form, 

stating  that  the  mass  flow  rate  m  is  constant  throughout  the  machine. 
One  has,  for  the  mass  flow  through  station  AB 


2n 


pw  cos(o-Y)b.B^e^ 

'ab  m 


rdl 


(2. 14 . 57 ) 


(e) 

where  B  is  the  blockage  factor  introduced  by  the  presence  of  the  end 
wall  boundary  layers.  If  6^  and  6*  are  respectively  the  hub  and  tip 

displacement  thickness  for  the  meridional  velocity  profiles  along  the 

(  6  ) 

station  AB,  then  the  factor  B  is  defined  by 


1-6* 
B  r 


mAB 


2tt 


r 

ji.+ 


6* 
A  °H 


pw  cos  (o-Y)brdl 
m 


(2.K.58) 


For  uniform  flow  conditions,  one  has  equating  (2.14.57)  and  (2.14.58) 

B(e)  .  -  SAfi  -  2urA6»  -  2irrB6$  (2.14.59) 

where  SAfi  is  the  geometrical  section  through  AB,  given  by 
r  r  ♦  r 

-  2*  j  rdl  -  2(-A  2-  B)  ArB  -  rft)2  ♦  (zQ  -  z ft)2  (2.14.60) 

AB 
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for  a  straight,  non  curved,  station  AB. 

(e) 

Equation  (2.4.59)  can  be  used  to  define  the  blockage  factor  B  when 
the  displacement  thicknesses  are  known  or  estimated. 


Determination  of  force  terms 

T 

The  total  force  component  can  be  obtained  from  the  tangential 

momentum  equation  (2.4.37)  and  the  definitions  (2.4.8)  to  (7.4.12). 


The  blade  force  component  fD1  is  defined  by 

Dl 

fBl  '  fBr  003  Y  -  fB*  Sln  1 

-  f00(tan  g'  .  sin  7  -  tan  e'  . 


cos  7) 


(2.4.61 ) 


If  a  blade  angle  is  defined  in  the  plane  (AB,0)  that  is  by  a  section 
along  the  calculation  station,  then 

,  ,9G, 

tan  e'  -  r(^) 


,30.  3r  ^  ,3e.  3z 

r(5F)  51  *  r<5i)  51 


(2.4.62) 


tan  e'  .  cos  7  -  tan  0'  sin  7 


Hence,  the  blade  force  component  fD1  can  be  written  as 

D  1 


fBl  "  "fB0  tan  Gb 


(2.4.63) 


The  friction  force  is  obtained  from  the  entropy  variation,  as 


'  f  0 


-Ff  r 

m 


w  w„  . 
m  0  3s 

w2  3m 


3s 

-  -T  sin  0  cos  0  t— 
dm 


(2.4.64) 


The  tangential  component  of  the  total  force  is  given  by 


T 

F  -  F  ♦  f 
0  rf0  B0 


m3-  . 

-  35  (rV 


(2.4.65) 
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and  from  equation  (2.14.63) 


fBi  •  'tan  ci  r  Is  <rV* F  f8tan  E'b 


w 


-tan  e'  .  —  (rv  )  -  T  1^-  sin0  cos0  tan  e' 
b  r  dm  0  3m  b 


(2. H.66) 


The  1-component  of  the  friction  force  is  obtained  from 
F. 


-  F, 

COS 

y  -  f,  sin  y 

fr 

fz 

w 

w 

-  -F„ 

(— 

cos  y  -  —  sin 

f 

w 

w 

’  -Fr 

cos 

0  sin  (o-T) 

(2. H.67) 


ds 

-T  .  cos20  sin  (o-T) 
dm 


Finally,  the  1-component  of  the  total  force  can  be  calculated  from 


F  -  F  +  f 
1  fl  B1 


(2.14.68) 


w 

-tan  e '  —  t—  (rv  )  -  T  [cos  0  sin(o-T)  ♦  sin  0  tan  e']cos  0 
b  r  3m  6  dm  b 


In  duct  regions,  outside  blade  rows,  fD1«  0  and  eventually  only  the 

Dl 

friction  force  might  be  accounted  for  due  to  wake  effects.  Then  F^  is 
given  by  equation  (2.14.67). 


Ill 


2.4.3;  Numerical  Techniques  for  Streamline  Curvature  Methods 


The  complete  numerical  procedure  for  streamline  curvature  methods  can  be 
described  by  the  following  steps  and  considerations. 


a.  Initial  mesh  point  distribution 

Given  the  meridional  cross-section  of  the  machine,  calculation  stations 
are  distributed  between  an  inlet  and  an  Outlet  section  and  the  current 
station  is  indicated  by  I  (figure  2.4.3). 

In  addition,  the  through-flow  is  divided  into  a  specified  number  of  M 

streamtubes  separated  by  (M  +  1 )  meridional  streamlines,  the  hub  wall 

being  streamline  j-1 ,  the  shroud  wall  is  streamline  j-M+1 .  Since  the 

streamline  slopes  o  and  curvatures  R  have  to  be  determined  at  every 

m 

point  (I,j),  the  positions  of  the  streamlines  separating  the  different 
streamtubes  has  to  be  known.  Therefore,  from  one  iteration  to  the 
other,  the  streamline  position  will  be  shifted  in  order  to  be  in 
accordance  with  the  newly  obtained  meridional  velocity  distribution  w^. 

This  distinguishes  the  streamline  curvature  method  from  classical  finite 
difference  or  finite  element  methods,  in  that  the  mesh  points  (I,j)  are 
not  fixed  in  the  flow  domain  but  are  displaced  along  the  I-stations  such 
as  to  follow  the  streamlines.  Hence,  at  any  interation,  the  j-lines 
will  be  considered  as  formed  by  streamlines  limiting  streamtubes  of 

constant  mass  flow  (ro/M)  and  displaced  at  the  end  of  the  iteration  after 

the  w  distributions  have  been  obtained  at  all  stations, 
m 

An  initial  location  of  the  streamlines  can  be  obtained  by  assuming 

uniform  w  distribution  along  each  calculation  station, 
m 


b.  Estimation  of  streamline  slopes  and  curvatures 

This  is  a  central  point  in  the  calculation  procedure,  since  the 
curvature  corresponds  to  a  second  derivative  of  the  function  r  -  r(z) 
along  a  streamline  and  its  numerical  estimation  from  a  limited  number  of 
points  generates  high  frequency  errors. 

One  has 


1  d2r  /  dz2 

R  ”  T7T 

m  (1  +  tan2o) 


(2.4.69) 


A  detailed  investigation  of  various  formulas  for  dr/dz  and  d2r/dz2  along 
a  streamline  has  been  presented  by  Wilkinson  (1970).  From  this  analysis 
results  that  the  best  choice  between 

i)  spline  curve  fitting  on  streamline  positions  and  (or)  on  its 
derivatives , 

ii)  standard  finite  difference  formulas  of  second  order  accuracy, 
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iii)  various  polynomial  fittings, 


is  undoubtly  obtained  for  the  standard  finite  difference  formulas. 
Hence  the  best  recommendation  is  to  use  the  following  formulas 


dr 

dz 


I.J 


rl+1,J 

*I+i.J 


rl.J  +  rU  " 

ZI,J  ZI,J  '  ZI-1,J 


tan  o 


I.J 


- 1 - (  rJ'A  -  ji-J  illihA  ) 

ZI+1.J  -  ZI-1,J  ZI*1.J'ZI,J  ZI,J"ZI-1,J 
2 


(2.4.70) 


(2.4.71  ) 


Another  variant  is  used  by  Hearsey  (1975),  to  the  same  order  of  accuracy 


o.  -  (  tan’1  ♦  tan'1  )  1 

I,J  zi*i,j  ZI,J  ZI,J  "  ZI-1,J  2 


(2.4.72) 


°I.J  *  °I.J 


I_.  -  °I.J  ~  °I.J 

Rm  ml+1,J  ~  mi-1,J 
2 


(2.4.73) 


The  theoretical  analysis  of  Wilkinson  has  been  confirmed  by  Novak 
(1976). 

Another  way  of  investigating  the  influence  of  the  estimation  of 
streamline  slopes  and  curvatures  on  the  solution  is  reported  by  J. 
Denton  in  Hirsch  and  Denton  (1981). 

A  simple  test  case  was  chosen,  consisting  of  a  single  blade  row  in  an 
axial  channel  with  a  hub  wall  at  constant  radius  and  a  casing  wall 
formed  by  two  circular  arcs  and  straight  lines  such  that  discontinuities 
in  wall  curvature  at  the  leading  and  trailing  edge,  from  circular  to 
straight  lines,  and  at  mid-chord  between  the  two  circular  arcs,  provide 
severe  tests  of  the  approximations. 

The  number  of  stations  was  varied  between  a  coarse  grid  with  no  internal 
stations  except  at  leading  and  trailing  edge  and  a  finer  grid  with  1,  3 
and  5  additional  internal  stations.  Using  the  same  second  order 
difference  formulas,  it  is  shown  that  with  no  internal  stations  the 
curvature  at  leading  edge  is  underestimated  leading  to  an  innacurate, 
overpredicted  estimation  of  axial  velocity  and  hence  of  flow  angles. 
The  inadequacy  of  computations  with  no  internal  stations  is  also 
confirmed  by  Novak  in  the  same  reference,  Hirsch  and  Denton  (1981),  page 
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245.  A  single  internal  station  covers  most  of  the  effects  of  curvature 
although  the  accuracy  is  moderably  improved  by  additional  stations. 
Similar  conclusions  on  the  necessity  of  including  internal  stations  have 
also  been  reached  by  Frost  (1970).  At  the  trailing  edge,  the  curvature 
is  shown  to  be  insensitive  to  the  number  of  internal  stations  as  a 
consequence  of  the  imposed  values  of  flow  angles. 

However,  the  estimation  of  the  streamline  curvature  at  trailing  edge  has 
a  strong  effect  on  the  prediction  of  the  downstream  flow  field.  This 
can  be  shown  in  various  ways,  through  equation  (2.4.29).  Assuming  that 
s,I  and  rvQ  as  well  as  wm  are  uniform  at  trailing  edge,  the  downstream 

variation  of  meridional  profiles  is  uniquely  influenced  by  the  radius  of 
curvature,  since  equation  (2.4.39)  reduces  to 


w 


tanU  -  K)  55  "m  -  r 

m 


It  can  also  be  shown  from  actuator  disk  theories,  Horlock  (1978)  that  a 

shift  in  streamline  of  (Ar)  will  generate  a  curvature  of  the  order  of 

Ar/(Az)2  where  Az  is  the  axial  distance  between  two  points  and 

innacuracies  on  1 /R  greatly  affect  the  meridional  velocity 

m 

distributions.  This  is  fully  confirmed  by  Denton's  model  calculation 
(see  figure  2.4.4)  where  the  effect  of  the  number  of  stations  used  to 
calculate  the  curvature  is  clearly  seen.  These  effects  will  be 
amplified  for  machines  with  strong  curvatures  such  as  centrifugal  flow 
machines  for  which  the  computation  could  be  excessively  sensitive  to  the 
numerical  estimation  of  the  curvature. 


c.  Determination  of  flow  properties  at  station  I  +  1 

The  flow  variables  are  supposed  to  be  known  at  station  I.  Within  the 
axisymmetric  approximation  the  stagnation  values  at  point  A  at  station  I 
+  1  are  obtained  from  the  known  values  at  station  I  in  the  following 
way. 

The  variations  of  angular  momentum  and  entropy  are  determined  as 
described  in  section  2.2.2. 

Within  a  blade  row  : 

v  “  vQ(m)  for  the  design  problem 

0  0 

B  -  B  (m)  for  the  analysis  problem 

s  -  s  (m) 

Outside  a  blade  row  : 


0 
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when  no  wake  or  mixing  effects  are  considered,  otherwise  equation 
(2.4.8)  has  to  be  applied.  The  stagnation  temperature  in  the  relative 
system  is  obtained  from  the  energy  equation  . 

From 


I 


A 


I 


C 


and  a  perfect  gas  assumption,  one  obtains  the  relative  stagnation 
temperature  TJ  through 


Evaluation  of  the  density 

The  density  is  evaluated  from  the  velocity  field  through 


P  -  P  •  ( 1 


w 


2 


2c  .T 
P 


1/Y-1 

7> 


(2.4.74) 


where  the  velocities  are  taken  from  the  previous  iteration  and  pj  is  the 
stagnation  density  in  the  relative  system. 

From  the  knowledge  of  the  density,  stagnation  temperature  TJ  and 
stagnation  pressure  pj  (obtained  from  the  entropy  equation,  see  equation 
(2.2.60))  one  can  estimate  the  static  temperature  and  speed  of  sound. 

The  previous  steps  allow  the  determination  of  all  quantities  in  the 
right  hand  side  of  equation  (2.4.40),  and  the  new  distribution  of 
meridional  velocities  along  station  1+1  can  be  obtained  from  the 
integration  of  the  ordinary  differential  equation  (2.4.40). 

d .  Integration  of  meridional  velocity  equation 

Various  methods  can  be  applied  to  integrate  equation  (2.4.40),  such  as  a 
Runge-Kutta  method  or,  when  linearising  the  coefficients,  a  first  order 
linear  equation  is  obtained,  which  can  be  exactly  solved. 

Whatever  method  is  applied,  an  initial  value  at  a  certain  point  is 
needed  to  initiate  the  numerical  integration.  This  reference  value  of 
w  ,  w  can  be  taken  at  hub,  shroud  or  any  point  along  station  (1+1). 
m  mo  m+1 

Generally  the  value  at  mid  station  is  chosen  that  is  the  value  at  J— 5— 

M+1  ^ 

and  the  integration  proceeds  in  both  directions  from  j— =—  to  (m+1 )  and 

from  this  mid-point  to  j-1 .  The  initial  value  of  w  is  arbitrary  and 

mo 
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will  be  corrected  iteratively  until  the  w  -w  (1)  distribution  satisfies 

m  m 

the  conservation  of  mass  flow  at  that  station  up  to  a  prescribed 
accuracy,  see  equation  (2. ^.57). 


This  is  not  a  trivial  step  since  the  sign  of  the  variation  of  mass  flow 

with  w _  depends  on  the  Mach  number  and  is  different  in  subsonic  and 

mo 

supersonic  regions,  following  equation  (2.M.31).  Following  Hearsey 

(1975),  Novak  (1976),  the  sign  of  the  correction  on  w  is  obtained  from 

mo 

the  variation  of  mass  flow  with  respect  to  wmQ.  Taking  into  account 

equation  (2.H.31)  and  (2.4.57)  •  one  can  write  for  fixed  stagnation 
conditions  and  tangential  velocity  v„ 


9m 


6w 


9w_ 


P(1 


m  r 


'  Mm} 

m 


6  w 


d  S 


(2.4. 75 ) 


where  dS  is  the  area 

dS 


element  normal  to  the  station 
-  2nrbB^e^dl  cos(o-Y) 


(2.4.76) 


If  the  flow  angle  is  fixed  instead  of  v  ,  one  has  from  w  =w  tan&  and 

u  o  m 

equation  (2.4.31 ) 

-  [  p(1  -  M(r)2)  --  dS  (2.4.77) 

m0  m0 


( p )  r2 

where  M  is  the  relative  mach  number  M  =w2/a2. 

The  mass  flow  rate  passes  through  a  maximum  at  sonic  meridional  Mach 

number  or  sonic  relative  Mach  numbers  depending  on  the  formulation  of 

the  problem  that  is,  v.  or  &  fixed.  This  is  an  alternative  formulation 

0 

for  the  properties  of  the  through-flow  equations  as  discussed  in  2.2.1. 

The  variation  ratio  6w  /6w  can  be  evaluated  either  in  an  approximate 

m  m0 

way,  through 


6  w  w 
m  m 


or  through  a  more  rigorous  way  by  a  direct  computation. 
With  the  first  assumption,  equation  (2. 4.77)  becomes 


(2.11.78) 
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(|e1,  -  Jp  j  p(,  .  M;,wmds 

m  m  * 


m  o  in 


vQ  fixed 
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(2.4.79) 


or 


3nT_ 

3wm 
m  a 


H 

m  * 


p(i 


M(r)2)w  dS 
m 


B  fixed 


(2.4.80) 


The  sign  of  the  variations  of  mass  flow  with  respect  to  w  is  positive 

m  © 

for  subsonic  Mach  numbers  and  negative  for  supersonic  velocities.  Also, 
the  local  mass  flow  in  a  streamtube  dS  is  maximum  at  sonic  conditions 
and  when  this  maximum  value  is  reached  the  streamtube  is  said  to  be 
choked.  The  large  sensitivity  of  mass  flow  with  respect  to  w^,  makes 

the  estimation  of  w  delicate  particularly  in  the  vicinity  of  sonic 

m  q 

velocities  or  choking  conditions. 

(r) 

If  the  Mach  number  M^  or  M  ,  according  to  the  problem  treated,  is 

subsonic,  the  computed  correction  on  wm  can  be  determined  from  a 

m  o 

relation  of  the  following  type,  for  subsonic  flows 


(w  ) 

m  o 


new 


(w  )  +  (m  -  m  )/( , 
m0  c  3w 


3m 


mc 


(2.4.81 ) 


where  m  is  the  calculated  mass  flow  with  the  current  distribution  of  w 
c  m 

and  the  given  w  .  The  same  relation  is  used  for  supersonic  velocities 
m  o 

when  the  supersonic  solution  is  desired. 


However,  when  a  subsonic  solution  for  instance,  is  looked  for  and  the 

calculated  derivative  8m/3w  is  negative,  which  indicates  that  the 

m0 

initial  chosen  value  of  w  is  on  the  supersonic  side  of  the  m-w 

m0  m 

relation  then  is  is  necessary  to  force  (w  )  to  move  towards  the 

m  new 

subsonic  branch  and  impose  a  relation  such  as 


) 

new 


(2.4.82) 


with  K-0.9  for  instance. In  the  opposite  case,  K  will  be  taken  greater 
than  1 . 
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Additional  recommendations  to  control  the  variations  of  w  are  given  by 

m0 

Hearsey  (1975)  and  Novak  (1976).  For  instance,  for  a  subsonic  solution, 

(w  )  should  not  be  lower  than  a  previous  estimation  w  which  led  to 
m0  _  m„ 

0  new  0 

a  supersonic  solution  and  inversely  when  a  supersonic  solution  is 

searched. 

With  the  new  value  of  (w  )  the  numerical  integration  of  the  momentum 

m0  . 
new 

equation  (2.4.40)  is  started  again  and  the  (w  )  values  recalculated 

m  q 

until  conservation  of  mass  flow  is  achieved  to  an  imposed  accuracy. 


Instead  of  the  approximation  (2.4.78)  a  better  estimation  of  (9w  /9w  ) 

m  m  0 

can  be  obtained  by  the  Taylor  expansion,  Hearsey  (1975) 


w  (1 )  -  w  ♦  (1-10)  aT  Wm 

m  m0  91  m 


CS.4.83) 


from  which  one  obtains 


9w 


m 
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9w 


m0  1 


(1”lo)  9w  (91  wm) 
m 


(2.4.84) 


The  derivative  9/9w  (9,w  )  can  be  calculated  from  the  momentum 

m  1  m 

equations  (2.4.40)  -  (2.4.42)  or  similar  forms,  (2.4.48)  to  (2.4.50). 


The  above  procedure  can  allow  the  computation  of  relative  supersonic 
flows  in  the  analysis  case,  when  the  type  of  solution  which  is  sought, 

subsonic  corresponding  to  9m/9wm  >  0,  or  supersonic  9m/9wm  <  0  is 

specified  by  the  user  at  every  station  inside  a  blade  row.  Otherwise, 
the  present  formulation  does  not  allow  to  select  the  correct  solution 
due  to  the  non-unique  relation  between  mass  flow  and  Mach  number;  see 
also  Marsh  (1971)  for  a  discussion  on  this  point. 

e .  Displacement  of  streamlines 


Once  the  new  distributions  of  meridional  velocities,  consistent  with 
mass  flow,  are  obtained  at  all  the  stations,  new  estimations  of  the  the 
streamline  positions  can  be  calculated. 

Referring  to  figure  2.4.3,  the  calculated  position  of  streamline  j  at 
station  I  can  be  obtained  from 


118  - 


m 


-  w 

m 


(A1'l,J  -  1I.J  calculated  "  1I,J|old 


hi in-e-w _ LA\2±*  (2.U.85) 

(—  w  ) 

V31  m'l.J 


A  detailed  theoretical  analysis  of  the  stability  of  the  streamline 
relocation  process  has  been  developed  by  Wilkinson  (1970).  From  this 
analysis  follows  that  the  repositioning  of  the  streamlines  has  to  be 
strongly  underrelaxed  in  order  to  achieve  convergent  iteration  cycles. 
Based  on  a  simple  model  of  a  swirling  flow  in  a  cylindrical  annulus,  an 
estimation  of  the  optimal  relaxation  factor  can  be  found.  A  simplified 
version  of  Wilkinson's  analysis  is  presented  by  Hearsey  (1975)  and  Novak 
(1976). 


According  to  the 
procedure 

problem 

solved. 

one  has  the  following 

relaxation 

1I,J 

new 

*  “ai.j 

1 

-lul  ) 

calc.  ,J|old 

(2.44.86) 

p  - 

1  ♦ 

K(1  -  M*) 
m 

(2.44.87) 

where  Ar  »  1  .  .  -  1.  .  is  the  length  of  the  station  and  Az  is  the 

shroud  hub 

meridional  distance  to  the  nearest  station. 

Analysis  problem 


p 


K(1 


,(r): 


(2.14.88) 


The  constant  K  is  in  both  cases  of  the  order  of  1/8  to  1/6.  These 

values  of  the  optimum  relaxation  factors  tend  to  zero  when  the  spacing 

between  stations  is  reduced,  that  is  when  the  mesh  is  made  finer. 

Therefore,  it  is  necessary  to  limit  the  spacing  between  stations  on  the 

lower  side  in  order  to  avoid  exessive  slow  convergence  rates.  The  new 

values  of  1  , I  are  then  rised  in  order  to  displace  the  streamlines 

I.J |new  i 

and  start  a  new  iteration  cyclus  with  updated  values  of  streamline 

slopes  and  curvature.  Another  variant  is  applied  by  Denton  (1978), 

where  the  curvature  is  underrelaxed  from  one  iteration  to  the  other. 


Summary 

The  streamline  curvature  method  solves  for  the  meridional  velocity 
distribution  along  arbitrary  calculation  stations  coupled  to  the  mass 
conservation  equation. 

The  streamline  positions  are  modified  from  one  iteration  to  the  next  and 
this  step,  coupled  to  the  estimation  of  the  slope  and  curvature  of  the 
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streamlines,  is  the  most  delicate  operation  in  the  whole  procedure. 
Practical  experience,  see  for  instance  Davis  and  Millar  (1975)  for  an 
interesting  account  of  some  available  user's  experience,  shows  that  with 
the  necessary  care  in  the  estimation  of  slopes  and  curvatures,  accurate 
results  can  be  obtained  enven  for  radial  and  mixed  flow 
machines , including  return  ducts.  A  recent  report  on  the  application  of 
the  streamline  curvature  method  to  radial  compressors  has  been  presented 
by  Casey  and  Roth  (1984). 

With  regard  to  Mach  number  limitations,  the  general  properties  discussed 
in  section  2.2.1  are  confirmed  by  the  properties  of  the  mass  flow 

derivatives  Om/Sw^)  through  the  sign  of  this  function.  As  pointed  out 

by  Marsh  (1971),  when  this  derivative  has  a  fixed  sign,  a  unique 
solution  is  obtained  for  the  iteration  on  the  meridional  distribution. 
This  is  the  case  for  subsonic  meridional  velocities  in  the  design 
problem  and  for  fully  subsonic  or  fully  supersonic  relative  Mach  numbers 
in  the  analysis  problems.  However  transonic  relative  velocities  can  be 
treated  if  the  user  specifies  which,  from  the  subsonic  or  supersonic 
solution  of  the  relation  between  mass  flow  and  Mach  number  is  to  be 
selected . 
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2.4.4.  Streamf unction  Methods 


An  alternative  to  the  streamline  curvature  method  is  provided  by  the 
introduction  of  a  streamfunction  ty,  in  order  to  satisfy  the  continuity 
equation.  This  leads  to  a  second  order,  quasi-linear ,  partial 
differential  equation  for  the  streamfunction,  which  can  be  discretized 
on  a  fixed  mesh,  either  by  finite  difference  or  by  finite  element 
methods . 


The  continuity  equation,  (2.4.1)  can  be  satisfied  for  any  continuous 
function  i|>,  by  setting 


1  34- 

w  *  -  — — 

z  pbr  3r 

_  _  3^ 

r  pbr  3z 
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Inserting  these  definitions  in  equation  (2.4.15)  for  instance,  leads  to 
the  following  streamfunction  equation 
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This  equation  can  be  applied  as  long  as  wz  does  not  vanish,  that  is  for 
axial  flow  configurations.  For  radial  flow  configurations,  the  axial 


momentum  equation  (2.H.17) 
streamfunction  equation 

can 

be 

used , 

leading 

to  an 
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(2.4.91 ) 

This  equation  in  turn  ceases 

to 

be 

valid 

when  the 

radial 

velocity 

component  vanishes. 

A  general  form,  valid  for  any  flow  configuration  and  which  can  therefore 
be  applied  in  mixed  flow  machines  with  axial-radial  transitions  is 
obtained  by  introducing  the  N-component  of  the  momentum  equation,  see 
section  2.3*4. 

Within  the  notation  conventions  of  this  section,  where  no  distinction  is 
made  between  the  strearnsheet  and  the  averaging  approaches  the 
streamfunction  equation  becomes 

Blade  region 
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(2.14.92) 
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In  the  axisymmetric  approximation  B',  e'  are  the  streamsurface  or  mean 
Camber  surface  angles  with  respect  to  the  axial  and  radial  directions 
res pectively . 


Duct  region 
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For  Incompressible  flows,  the  same  equations  remain  valid  with 


s  -  0 
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In  the  following  we  will  consider  either  of  these  three  forms,  recalling 
that  equation  (2.14.92)  is  valid  for  all  flow  configurations. 

Hence,  the  streamfunction  equation  will  be  written  as 


3  /  _J 3j(\  +  3_  /_J 3  V  t 
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(2.14.95) 


An  alternative  form  is  obtained  by  working  out  the  left-hand  side  into  a 
Laplace  operator 
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or 


3  3 

Ay  -  q. pbr  -  wp  —  (pbr)  +  wz  —  (pbr) 


(2.14.97) 


Equation  (2.H.95)  is  elliptic  for  subsonic  relative  Mach  numbers  in  the 
analysis  case  and  for  subsonic  meridional  Mach  numbers  in  the  design 
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option  when  rv  is  specified.  This  results  from  the  general  analysis  of 

section  2.2.1  but  can  also  be  shown  directly  on  equation  (2. 4. 95)  using 
equation  (2.2.81). 


In  the  analysis  case,  the  (rv0)  derivatives  of  the  right-hand-side 

contribute  to  the  type  of  the  streamf unction  equation  through  equation 
(2.2.83).  For  the  design  problem,  applying  equation  (2.2.81)  gives 
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and  one  can  write  the  streamf unction  equation  in  the  following  form, 
with  double  subscripts  on  ip  indicating  second  derivatives 
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where  D  contains  the  terms  of  lower  order  in  the  ^-derivatives, 
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This  equation  is  elliptic  for 
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In  the  analysis  case,  when  the  flow  angle  is  imposed,  the  streamf unction 
equation  becomes,  Wu  (1952),  with 


0 

tans  -  — 
w 

z 


(1  +  tan2B)  (1  -  M2)  ^ 
r  rr 


2(1+  tan2e)  M  H  i|i 

r  z  rz 


+  (1 


M2)  4. 


zz 


D 


which  is  elliptic  for 
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that  is  for  relative  Mach  numbers  M  lower  than  one,  and  hyperbolic 
for  supersonic  relative  Mach  numbers. 


The  determination  of  density 

One  of  the  central  problems  with  the  application  of  streamfunction 
methods  for  compressible  flows  is  connected  to  the  non-unique  relation 
between  the  streamfunction  and  the  density. 


From  the  isentropic  relation  between  static  and  relative  stagnation 
variables,  one  has 
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where 
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Dtp 

1 3r 

Dz 
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and  where  the  stagnation  quantities  are  considered  In  the  relative 
system. 


For  fixed  values  of  the  streamfunction  gradient  |Vtp|2,  relative 
stagnation  enthalpy  Hq  and  stagnation  density  pQ,  equation  (2.4.101)  has 

generally  two  solutions,  one  corresponding  to  a  supersonic  flow,  the 
other  to  a  subsonic  relative  flow. 

Indeed,  equation  (2.4.101)  can  be  written  under  the  form,  valid  for  the 
design  option,  v.  or  wQ  fixed, 
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(2.4.104) 


For  the  analysis  option,  (5  fixed,  one  has 
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The  relations  p/p^”F(p/p^)  are  represented  on  figure  2.^.5  and  the 

solutions  can  be  found  as  the  intersections  SI  and  S2  between  the 
function  and  the  first  quadrant  bissector. 


The  function  F  has  an  asymptotic  limit  for  large  values  of  p  given  by 
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The  number  of  solutions  for  the  density  depends  on  the  value  of  the 
constant  A.  Generally  one  will  have  two  values  of  the  density,  given  by 
the  points  SI  and  S2  of  figure  2.^.5,  for  a  given  value  of  the  constant 
A. 


Point  SI  is  the  subsonic  solution  while  point  S2  corresponds  to  a 
supersonic  velocity.  Points  C  are  defined  by  the  condition  that  the 
tangents  to  F  equal  one.  If  point  C  lies  below  the  line  F-p/p^  for 

instance  curve  2,  there  will  be  no  solutions  for  the  density 
corresponding  to  the  current  values  of  the  constants  A.  Indeed,  the 
condition  to  be  imposed  in  order  to  have  two  solutions 
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where  (p/p')r  is  defined  by  the  condition 
0  c 
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the  notation  F'  indicating  a  derivative  with  respect  to  (p/p^). 
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The  limiting  case  being  reached  for  curve  3  where  F(C3)  *  (p/p')r,. 

O  U  j 

From  the  definitions  above,  one  has 


(2.4.110) 


Considering  a  perfect  gas  assumption  and  the  particular  limiting 
position  of  curve  3,  for  which  F  -  p/pQ  in  point  C3,  one  obtains  the 

conditions 
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This  shows  that  point  C3  of  curve  3  corresponds  to  the  critical  Mach 
number  namely 
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(2.4.112) 


and  that  the  critical  density  ratio  is  reached 

(p/Po)*C3  -  (Z£f1/(1M)  (1 
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The  lowest  value  of  (p/pj)  is  obtained  by  the  intersection  of  the  F 
curve  with  the  horizontal  axis,  that  is  by  F-0.  Hence 
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The  critical  value  of  A,  A*  corresponds  to  curve  C3,  with  the  condition, 
deduced  from  equation  (2.4.110),  that 
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with  (p/pj)*  given  by  equation  (2.4.113). 

The  condition  for  the  density  relation,  equation  (2.4.103).  to  have 
solutions,  can  be  expressed  by  the  condition  that  the  values  of  A 
corresponding  to  curve  Cl  should  be  larger  than  A*. 

This  is  actually  a  condition  on  the  mass  flow  per  unit  area  pwn,  since 

the  critical  mass  flow  is  attained  for  w*-c*  the  speed  of  sound 
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Hence,  the  condition  A  <  A*  becomes  in  all  cases,  with  either  w.  or  tan 
B  fixed, 
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In  practical  computations,  when  the  computed  streamfunction  is  far  from 
the  converged  value,  a  careful  check  has  to  be  performed  in  order  to 

insure  that  the  current  approximation  of  | Vip |  2  is  not  too  high  such  as 

to  violate  the  condition  (2.4.117). 

When  (2.4.117)  is  satisfied,  two  solutions  occur,  corresponding  to  the 
Subsonic  branch  point  SI ,  or  to  the  supersonic  branch  of  the  mass  flow 
relation  versus  density  or  Mach  number.  In  order  to  solve  the  density 
equation  and  obtain  a  unique  solution,  one  can  apply  a  Newton  method  or 
a  fixed  point  method.  The  former  will  allow  to  obtain  either  of  the  two 
solutions,  if  the  user  knows  which  one  to  select,  while  the  latter  will 
always  lead  to  the  subsonic  solution. 
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Indeed,  writing  equations  (2.11.103)  or  (2.4.105)  under  the  form 


(2.4.118) 


where 
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the  Newton  iteration  reads 
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The  function  Gn  is  the  value  of  G  for  p-p°  and  the  Jacobian  G'  is 
defined  by 


(2.4.121) 


As  expected,  the  zero  of  the  Jacobian  G'  is  obtained  at  the  critical 
mass  flow  condition  (2.4.115)  but  the  jacoblan  will  remain  of  constant 
sign,  positive  or  negative  for  subsonic  or  supersonic  flows 
respectively. 

Graphically,  it  is  easy  to  see  that  the  Newton  iteration  will  converge 
to  the  subsonic  or  the  supersonic  solution  depending  on  the  value  chosen 
to  initiate  the  iteration,  figure  2.4.6. 

If  the  iteration  starts  at  A1  at  a  subsonic  velocity,  the  Newton 
procedure  will  converge  to  SI,  while  starting  with  A2  at  a  supersonic 
point  will  lead  to  the  supersonic  solution  S2. 

On  the  other  hand,  the  fixed  point  iteration 


(p/pj)n+1  “  F(pn/pJ) 


(2.4.122) 


converges  always  to  the  subsonic  solution  SI  as  illustrated  on  figure 
2.4.6,  excepted  when  the  initial  value  is  left  of  the  supersonic 
solution.  In  this  case  the  iteration  diverges. 


Boundary  conditions  for  the  streamf unction  equation 

The  boundary  conditions  for  the  solution  of  the  streamf unction  equation 
are  easily  defined  and  applied. 
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At  the  Inlet  station,  the  flow  is  assumed  to  be  known  and  therefore  it 
is  easy  to  determine  the  streamfunction  distribution. 

For  1-1,  figure  2.4.7, 


is  given  for  all  values  of  J,  from  the  inlet  velocity  distribution. 

At  the  outlet  station  it  is  nearly  always  possible  to  select  a  station 
which  is  situated  in  a  region  of  uniform  meridional  flow  direction. 
Hence,  if  the  last  station  is  set  perpendicularly  to  this  direction,  one 
can  assume  the  normal  derivative  to  be  zero.  That  is 

*  0  at  the  last  station  I-n  (2.4.124) 

an 


where  n  is  the  normal  direction. 


The  hub  and  shroud  walls  can  be  taken  as  streamlines  if  no  end  wall 
boundary  layer  is  to  be  considered.  In  this  case  one  can  select 
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However  if  an  end  wall  boundary  layer  computation  is  performed  or  if  the 
wall  layer  thicknesses  are  given,  one  can  either  displace  the  hub  and 
shroud  wall  points  by  an  amount  equal  to  the  local  displacement 
thickness  maintaining  the  above  boundary  conditions,  or  one  can 
introduce  a  fictive  mass  flux  compensating  for  the  mass  flow  defect  due 
to  the  shear  layers. 


If  and  6*  represent  the  hub  and  shroud  displacement  thicknesses  with 
respect  to  the  meridional  velocity  profiles,  one  can  impose 
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In  this  case  the  physical  end  walls  are  not  streamlines  anymore. 


2.4.5.  Finite  Difference  Solutions  of  the  Streamfunction  Equation 


The  discretization  of  the  streamfunction  with  finite  differences  has 
been  initially  proposed  by  Wu  (1952).  The  approach  of  Wu  was  closely 
followed  and  developed  to  practical  computer  programs  by  Marsh  (1968), 
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Davis  (1971),  (1975),  Katsanis  and  McNally  (1973),  (1977),  Biniaris 
(1975),  Bosnian  and  El  Shaarawi  (1977). 

All  these  authors  use  central  difference  schemes  and  hence  the  methods 
are  limited  to  subsonic  flow  regions,  either  to  meridional  subsonic 
velocities  for  the  design  option,  or  to  relative  subsonic  velocities  in 
the  analysis  case.  The  densities  are  obtained  by  the  fixed  point 
iteration  scheme  in  all  the  presently  published  methods.  Of  the  above 
mentioned  methods,  Davis  (1971)  and  Katsanis  and  McNally  (1973)  contain 
full  documented  codes.  In  particular  the  latter,  called  MERIDL,  is 
extensively  documented. 

In  most  methods  an  arbitrary  grid  is  generated,  although  Katsanis  and 
McN&lly  work  with  a  curvilinear  orthogonal  grid  and  Biniaris  uses  a 
cartesian  grid  in  the  meridional  surface  (r,z).  Second  order  central 
differences  are  applied  by  these  two  authors,  while  Marsh  (1968)  and 
Davis  (1971)  apply  up  to  third  order  central  difference  formulas  on 
arbitrary  grids. 

The  algebraic  matrix  system  is  mostly  solved  by  LU  decompositions  with 
the  exception  of  Katsanis  and  McNally  (1973)  and  Biniaris  (19J5)  who 
apply  successive  overrelaxation  techniques. 

2. *1.6.  Finite  Element  Solutions  of  the  Streamfunction  Equation 


Finite  element  solution  procedures  for  the  through-flow  problem  have 
been  introduced  by  Adler  and  Krimerman  ( 1 97 ^ ) ,  Hirsch  and  Warzee  (197*0. 
The  first  authors  define  triangular  meshes  with  linear  elements  while 
Hirsch  and  Warzee  apply  biquadratic  elements  leading  to  third  order 
accuracy  for  the  streamfunction  and  hence,  second  order  accuracy  for  the 
velocities . 


SUMMARY 


The  specific  approximations  introduced  in  order  to  calculate  multistage 
turbomachinery  flows,  are  a  combination  of  a  distributed  loss  model  and 
a  quasi-three-dimensional  decomposition.  Two  approaches,  the 
streamsurface  method  or  the  averaging  technique  can  be  applied  to  define 
the  quasi-3D  formulation.  The  latter  allows  a  through  flow  definition 
without  ambiguity  or  arbitrary  selected  streamsurface  shapes.  However 
it  does  not  permit  a  reconstruction  of  the  full  three  dimensional  flow. 

The  quasi~3D  approximation  enables  to  reduce  the  problem  to  a  succession 
of  two-dimensional  flows,  basically  of  inviscid  nature.  The  interaction 
between  the  through-flows  and  the  blade-to-blade  flows  can  be  extended 
to  the  transfer  of  losses  and  flow  angles  from  the  blade-to-blades  to 
the  meridional  flow.  However,  this  can  be  replaced  by  empirical 
information,  which  is  required  in  all  cases  if  some  three-dimensional 
effects  have  to  be  accounted  for. 

The  specific  properties  and  solution  methodology  related  to  through-flow 
computations  have  been  presented  and  different  resolution  techniques 
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have  been  described.  None  of  the  methods  described  however,  is  able  to 
treat  correctly  relative  supersonic  flows  in  the  analysis  mode.  This 
would. require  a  different  simulation  of  the  problem,  see  for  instance 
Spurr  (1981). 

Two  methods  are  currently  applied  for  through-flows  :  streamline 
curvature  and  streamsurface  methods.  The  latter  can  be  discretized 
either  by  finite  differences  or  by  finite  elements. 

Solution  techniques  for  blade-to-blade  flows  have  not  been  discussed  in 
this  chapter  since  they  do  not  require  specific  methods.  Depending  on 
the  level  of  approximation,  potential,  Euler  or  Navler-Stokes  models  can 
be  applied  to  the  blade-to-blade  flows  following  the  methods  presented 
in  the  corresponding  chapters  dealing  with  the  discretization  of  these 
flow  models.  Several  examples  of  blade-to-blade  flow  computations  have 
already  been  shown  in  chapters  2  and  7  and  other  examples  will  be 
presented  in  the  following  chapters. 
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F igure  2.1.1.  :  Complex  blow  Phenomena  in  compressor  rotor  blade  row. 


Figure  2,1,2..  :  Flow  visualization  of  horseshoe  vortex  in  the  leading  edge 

region  of  a  turbine  cascade.  Courtesy,  C.  Sieverding,  Von  * 

Karman  Institute,  Belgium. 
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Figure  2.1.3.  :  Transonic  pressure  field  in  a  compressor  cascade.  From  AGAKD 
LS  83. 
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Figure  2.1.14.  :  Transonic  flow  shock  model  associated  to  figure  2.3.  From 
AGARD  LS  83 
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Figure  2.1.5.  :  Meridional  section  of  a  two-stage  compressor. 


Fig.  4  Second  ststs  stator  sxlt  total  praaaurs  contours,  aCPT  -  0.1 


Figure  2,1,6.  :  Stagnation  pressure  distribution  at  outlet  of  the  stator  of 
the  second  stage  of  an  axial  compressor.  Courtesy  R.  Dring, 
United  Technology  Research  Center,  USA. 
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F igare  2.2.1.  :  Angles  of  S2  surface  with  cylindrical  coordinates  (r,6,z). 


Figure  2.2.2.  :  Definition  of  flow  angles  and  velocity  components. 
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Figure  2.2.3»  :  Cross  section  of  a  blade  row  by  an  axisymmetric  surface 


Figures  2.2.M.  and  2.2.5.  :  Through  flow  interaction  with  blade  to  blade 

flow  for  choked  compressor  and  turbine  configurations 
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Figure  2.2.6.  :Scheme  of  three  dimensional  flow  computation  with  the 
-  streamsurface  formulation,  from  Adler  and  Krimerman  (1978). 
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Figure  2.2.7.  rComparison  of  ip  -  0.5,  mean  streamline  angles  with  mean 
camberline  angles,  from  Novak  and  Hearsey  (1977) 


Figure  2.2.8.  :Meridional  and  blade  to  blade  cross  section  of  contoured 
turbine  nozzle,  from  Novak  and  Hearsey  (1977) 


(o)  Station  A  A  (b)  Stotlon  T.  E. 


Figure  2.2.9.  :Radial  static  pressure  variation  for  leaned  and  unleaned 
annular  turbine  nozzle,  from  Novak  and  Hearsey  (1977) 


Meridional  Distance -in. 

Figure  2.2.10.  :Meridional  distribution  of  or  spacific  mass  flow  ratio 

distribution  for  turbine  nozzle,  from  Novak  and  Hearsey  (1977) 
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Figure  2.3.1.  :  Meridional  and  blade  to  blade  sections  of 
axial  flow  turbomachine. 
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Figure  2.^.1 .  :  Meridional  cross  section  of  multibladed  turbomachine. 
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Figure  2.^.2.  :  Layout  of  meridional  section  with  arbitrary  calculation 
stations . 
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Figure  2.4.3.  :  Typical  lay-out  for  streamline  curvature  methods 

I  indicates  the  station  number,  J  the  streamline  number. 


y/h 


Figure  2.4.4  :  Influence  of  curvature  estimation  on  meridional  velocity 
distribution  in  Denton's  model  problem;  from  AGARD  AR175, 
(1981 ). 
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Figure  2.4.5.  :  Relation  between  density  and  streamf unction  gradients. 


Figure  2.4.6.  :  Iterative  solution  procedure  for  the  density  equation 
(2.4.118),  in  function  of  streamfunction  gradients. 
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Figure  2.)).7.  :  Meridional  section  of  mixed  flow  machine  and  grid  layout. 
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8.  National  Aeronautics  &  Space  Administration 

Lewis  Research  Center 

21000  Brookpark  Road 
Cleveland,  OH  44135 

Attentions  Chief,  Internal  Fluid  Mechanics  Division  3 

Dr.  J.  Adamczyk,  MS  5-9  1 

Library  1 

9.  Dr.  Bao  B.  Hwang  1 


Propulsion  and  Auxiliary  Department 
David  Taylor  Naval  Ship 
Research  and  Development  Center 
Annapolis,  MD  21402 

1 0 •  Library  1 

General  Electric  Company 
Aircraft  Engine  Technology  Division 
DTO  Mail  Drop  H43 
Cincinnati,  OH  45215 

1 1 •  Library  1 

Pratt  &  Whitney  Aircraft  Group 
Post  Office  Box  109600 
West  Palm  Beach,  FL  33410-9600 

12.  Library  1 

Pratt-Whitney  Aircraft  Group 

East  Hartford,  CT  06108 

13.  Library  1 

Curtis  Wright  Corporation 

Woodridge,  NJ  07075 

1 4 •  Library  1 

AVCO/Lycoming 
550  S.  Main  Street 
Stratford,  CT  06497 

15.  Library  1 

Teledyne  CAE,  Turbine  Engines 

1330  Laskey  Road 
Toledo,  OH  43612 

16.  Library  1 

Williams  International 

P.  O.  Box  200 
Walled  Lake,  MI  48088 

1 7 •  Library  1 

Detroit  Diesel  Allison  Division  G.M.C. 

P.  O.  Box  894 
Indianapolis,  IN  46202 
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18.  Library 

Garrett  Turbine  Engine  Company 
111  S.  34th  Street 
P.  0.  Box  5217 
Phoenix,  AZ  85010 

19.  Professor  J.  P.  Gostelow  1 

School  of  Mechanical  Engineering 

The  New  South  Wales  Institute  of  Technology 

New  South  Wales 

AUSTRALIA 

20.  Dr.  G.  J.  Walker  1 

Civil  and  Mechanical  Engineering 

Department 

The  University  of  Tasmania 
Box  25 2C 

GPO  Hobart,  Tasmania  7110 
AUSTRALIA 

21 •  Professor  F.  A.  E.  Breugelmans  1 

I ns ti tut  von  Karma n  de  la  Dynamique 
des  Fluides 

72  Chausee  de  Waterloo 
1640  Rhode-St.  Genese 
BELGIUM 

22.  Professor  Ch.  Hirsch  2 

Vrije  Universiteit  Brussel 

Pleinlaan  2 
1050  Brussels 
BELGIUM 

23.  Director  1 

Gas  Turbine  Establishment 

P.  0.  Box  305 
Jiangyou  County 
Sichuan  Province 
CHINA 

24.  Professor  C.  H.  Wu  1 

P.  0.  Box  2706 

Beijing  100080 
CHINA 

25.  Director,  Whittle  Laboratory  1 

Department  of  Engineering 

Cambridge  University 
ENGLAND 
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26.  Professor  Jacques  Chauvin 
Urdversite  d* Aix-Marseille 
1  Rue  Honnorat 
Marseille 
FRANCE 

27.  Mr.  Jean  Fabri  1 

ONERA 

29,  Ave.  de  la  Division  Leclerc 

92  Chatillon 

FRANCE 

28.  Professor  D.  Adler  1 

Technion  Israel  Institute  of  Technology 

Department  of  Mechanical  Engineering 

Haifa  32000 

ISRAEL 

29.  Dr.  P.  A.  Paranjpe  1 

Head,  Propulsion  Division 

National  Aeronautics  Laboratory 
Post  Bag  1700 
Bangalore  -  17 
INDIA 

30.  Dr.  W.  Schlachter  1 

Brown,  Boveri  Company  Ltd. 

Dept.  T-T 

P.  O.  Box  CH-5401  Baden 
SWITZERLAND 

31 •  Professor  Leonhard  Fottner  1 

Department  of  Aeronautics  and  Astronautics 
German  Armed  Forces  University 
Hochschule  des  Bundeswehr 
Werner  Heisenbergweg  39 
8014  Neubiberg  near  Munich 
WEST  GERMANY 

32.  Professor  Dr.  Ing.  Heinz  E.  Gallus  1 

Lehrstuhl  und  Institut  feur  Strahlantiebe 

und  Turbourbeitsmashinen 
Rhein. -Westf.  Techn.  Hochschule  Aachen 
Templergraben  55 
5100  Aachen 
WEST  GERMANY 

33.  Dr.  Ing.  Hans-J.  Heinemann  1 

DFVLR-AVA 

Bunsens trasse  10 
3400  Geottingen 
WEST  GERMANY 
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34.  Dr.  H.  Weyer  1 

DFVLR 

^  Linder  Hohe 

505  Porz-Wahn 
WEST  GERMANY 
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35.  Dr.  Robert  P.  Dring  1 

United  Technologies  Research  Center 

East  Hartford,  CT  06108 

36.  Chairman  1 

Aeronautics  and  Astronautics  Department 

31-265  Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

37.  Dr.  B.  Lakshminarayana  1 

Professor  of  Aerospace  Engineering 

The  Pennsylvania  State  University 

233  Hammond  Building 

University  Park,  Pennsylvania  16802 

38.  Dr.  Steven  Shamroth  1 

Scientific  Research  Associates,  Inc. 

PO  Box  498 

Glastonbury,  CT  06033 

39.  Professor  Alan  H.  Epstein  1 

Gas  Turbine  Laboratory 

Massachusetts  Institute 
of  Technology 
Cambridge,  MA  02139 

40.  Mechanical  Engineering  Department 
Virginia  Polytechnic  Institute  and 

State  University 
Blacksburg,  VA  24061 

Attn:  Professor  W.  O' Brian  1 

Professor  H.  Moses  1 


41 .  Professor  T.  H.  Okiishi  1 

Professor  of  Mechanical  Engineering 

208  Mechanical  Engineering  Building 
Iowa  State  University 
Ames,  Iowa  50011 

42.  Dr.  Fernando  Sis to  1 

Professor  and  Head  of  Mechanical 

Engineering  Department 
Stevens  Institute  of  Technology 
Castle  Point 
Hoboken,  NJ  07030 
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43.  Dr.  Leroy  H.  Smith,  Jr. 

Manager,  Compressor  and  Fan 

Technology  Operation 
General  Electric  Company 
Aircraft  Engine  Technology  Division 
DTO  Mail  Drop  H43 
Cincinnati,  OH  45215 
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44.  Dr.  W.  Tabakoff  1 

Professor,  Department  of  Aerospace 

Engineering 

University  of  Cincinnati 
Cincinnati,  OH  45221 

45.  Mr.  P.  Tramm  1 

Manager,  Research  Labs 

Detroit  Diesel  Allison  Division 
Genteral  Motors 
P.  0.  Box  894 
Indianapolis,  IN  46206 

46.  Professor  Frank  Moore  1 

Sibley  School  of  Mechanical  & 

Aerospace  Engineering 
291  Gramman  Hall 
Cornell  University 
Ithaca,  NY  14853 

47.  Mr.  P.  F.  Yaggy  1 

Director 

U.  S.  Army  Aeronautical  Research  Laboratory 
AMES  Research  Center 
Moffett  Field,  CA  94035 

48.  Library  2 

Code  1424 

Naval  Postgraduate  School 
Monterey,  CA  93943 

49.  Office  of  Research  Administration  1 

Code  01 2 

Naval  Postgraduate  School 
Monterey,  CA  93943 

50.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22314 

51 •  Naval  Postgraduate  School 

Monterey,  CA  93943 

Attn:  Professor  M.  F.  Platzer  (67PL)  1 
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